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RESUMO: A espécie florestal Mimosa 
scabrella Benth. conhecida popularmente 
como bracatinga, pertence à família Fabaceae. 
O uso de sementes de qualidade é um fator 
determinante para o sucesso do empreendimento 
florestal, e a principal característica da qualidade 
a ser analisada é a capacidade germinativa 
das sementes. Assim o presente trabalho teve 
como objetivo avaliar a qualidade fisiológica de 
sementes de dez matrizes de M. scabrella, bem 
como avaliar o coeficiente de repetibilidade e 
a dissimilaridade genética em características 
biométricas das sementes. As sementes foram 
coletadas de uma população do município de 
Curitibanos, no mês de fevereiro de 2018. Foram 
realizadas análises como: peso de mil sementes 
(PMS), teor de água (TA), condutividade elétrica 
(CE) e teste de germinação (G) onde foram 
avaliados índice de velocidade de germinação 
(IVG) e comprimento de plântulas (CMP). Todas 
as análises estatísticas foram realizadas por 
meio do software GENES. Dentre as matrizes 
avaliadas por meio das análises de peso de mil 
sementes, teor de água, germinação, índice de 
velocidade de germinação, comprimento de 
plântula e condutividade elétrica, as matrizes 3 e 9 
apresentaram menor qualidade fisiológica, sendo 
indicada a formação de um lote de sementes de 
com as matrizes 1,2,4,5,7 e 10. O coeficiente 
de repetibilidade foi alto para a característica 
diâmetro longitudinal da semente, intermediário 
para o peso da semente, e, baixo para diâmetro 
equatorial e espessura da semente. Para garantir 
90% de determinação, é necessária a medição 
de 21 sementes, para compreensão de todas as 
características. Com base nas características 
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biométricas avaliadas, indica-se o cruzamento entre as matrizes 4 e 7 por apresentarem 
maior dissimilaridade genética. Ainda, com base na distância de Mahalanobis foi possível 
formar dois grupos, demonstrando uma variabilidade restrita da população estudada. 
PALAVRAS - CHAVE: Bracatinga, divergência genética, espécie florestal, biometria. 


PHYSIOLOGICAL QUALITY, REPEATABILITY AND GENETIC DISSIMILARITY 
FOR BIOMETRIC CARACTERISTICS OF SEEDS OF Mimosa scabrella BENTH 


ABSTRACT: The forest species Mimosa scabrella Benth, popularly known as bracatinga, 
belongs to the family Fabaceae. The use of quality seeds is a determinant factor for the 
success of the forest growth, one of the main characteristics of the germinative capacity of 
the seeds. Thus the present work had the objective of evaluating the physiological quality of 
seeds of ten matrices of M. scabrella, as well as to evaluate the repeatability coefficient and 
genetic dissimilarity in biometric characteristics of the seeds of the study population. The 
seeds were collected from a population of the municipality of Curitibanos, in the month of 
February, 2018. Analyzes were performed as: thousand seed weight (PMS), water content 
(TA), electrical conductivity (CE) and germination (G), where they were evaluated germination 
speed index (IVG) and seedling length (CMP). All statistical analyzes were performed using 
the GENES software. Among the matrices evaluated by means of the tests of thousand seed 
weight, water content, germination, germination speed index, seedling length and electrical 
conductivity 3 and 9 presented lower physiological quality, indicating the formation of a batch 
of seed of qualities of the matrices 1, 2, 4, 5, 7 and 10. The repeatability coefficient was high 
for the characteristic longitudinal diameter of the seed, intermediate for the weight of seed, 
and, down to equatorial diameter and seed thickness. To ensure 90% determination, it is 
necessary to measure 21 seeds, to understand all the characteristics. Based on the biometric 
characteristics evaluated, the crosses between matrices 4 and 7 are indicated because they 
present greater genetic dissimilarity. Also, based on the Mahalanobis distance, it was possible 
to form two groups, demonstrating a restricted variability of the studied population. 
KEYWORDS: Bracatinga, genetic divergence; forest species, biometrics. 


INTRODUÇÃO 


A bracatinga é um exemplo de espécie nativa com potencial para tornar-se viável 
comercialmente. Isso será possível devido à natureza perene, e ao desenvolvimento de 
estudos genéticos de populações e de trabalhos de melhoramento genético à longo prazo. 
Nessa perspectiva, a realização de trabalhos de melhoramento genético de espécies 
nativas, no meio acadêmico, é possibilitada pela continuidade dos estudos realizados numa 
mesma população, em períodos diferentes por pesquisadores distintos. Sendo assim, 
através das sucessivas gerações melhoradas, pode-se obter genótipos mais produtivos a 
serem disponibilizado aos produtores (NASCIMENTO, 2010). 

Através de testes de germinação realizados em laboratórios é possível avaliar o 
potencial fisiológico das sementes. O uso de sementes de qualidade é um fator determinante 
para o sucesso do empreendimento florestal, e a principal característica relacionada à 
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qualidade a ser analisada é a capacidade germinativa das sementes, pois, sem esta, a 
semente não possui valor para a semeadura, e dela também dependem a qualidade das 
mudas e a implantação de um reflorestamento com bons resultados (GONÇALVES et al., 
2013). 

De forma complementar ao teste de germinação são realizados testes de vigor, 
sendo estas ferramentas importantes para determinar a qualidade fisiológica de sementes 
florestais, sendo aperfeiçoados constantemente. Dentre estes, podem ser citados os testes 
de envelhecimento acelerado, condutividade elétrica, lixiviação de potássio e tetrazólio. 

Juntamente, o conhecimento sobre produção e tecnologia de sementes florestais 
assume importância fundamental no processo de manejo, conservação e melhoramento 
genético dessas espécies. Sendo assim, é de extrema importância a escolha de plantas 
matrizes sadias, com boa capacidade de produção de sementes, associada a um 
monitoramento adequado do processo de produção e coleta de sementes (ROVERI NETO, 
2014). 

Em programas de melhoramento, os testes genéticos demandam seleção de plantas 
superiores, as quais apresentem características interessantes para recombinação, e ao 
escolher um genótipo, acredita-se que sua superioridade perdure durante toda a sua vida 
(CRUZ et al., 2012). Pode-se estimar o coeficiente de repetibilidade realizando várias 
medições em um mesmo indivíduo com variações no tempo ou no espaço. Através desse 
coeficiente é possível avaliar se a seleção estabelecida para uma característica fenotípica 
será confiável, isto é, se os genótipos selecionados manterão sua superioridade. Possibilita 
ainda, determinar o número de medições necessárias em cada indivíduo (DANNER et al., 
2010). 

Desta forma, o objetivo deste trabalho foi avaliar a qualidade fisiológica de sementes 
de M. scabrella, e estimar o coeficiente de repetibilidade para características biométricas. 


MATERIAL E MÉTODOS 


Descrição da Área de Coleta das Sementes 


Os frutos foram colhidos de 10 matrizes localizadas no entorno do Campus de 
Curitibanos da Universidade Federal de Santa Catarina (27º17'7” Sul, 50º32'3” Oeste) a 
uma altitude de 987 m, no município de Curitibanos — SC. As sementes de M. scabrella 
foram coletadas no mês de fevereiro de 2018, diretamente na árvore através do método de 
vibração da copa com o auxílio de podão, fazendo com que as sementes caíssem sobre a 
lona instalada na projeção da copa (MAZUCHOWSKI, 2014). 

Após a coleta as sementes foram beneficiadas manualmente com o auxílio 
de peneiras, sendo descartadas as sementes imaturas, deterioradas ou danificadas. 
Posteriormente, foram determinados o teor de água pelo método de estufa a 105 + 
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3ºC durante 24 horas, com 4 repetições, e o peso de mil sementes. O teor de água foi 
realizado para cada matriz, porém o peso de mil sementes realizou-se somente para sete 
matrizes, devido à falta de sementes para realização desta análise. Ambas as análises 
foram realizadas conforme recomendações das Regras de Análise de Sementes — RAS 
(BRASIL, 2009). O trabalho foi desenvolvido no Laboratório de Biotecnologia e Genética da 
Universidade Federal de Santa Catarina — Campus de Curitibanos. 


Teste de Germinação, Índice de Velocidade de Germinação e Comprimento de 
Plântulas 

O teste de germinação foi realizado de acordo com as instruções para análise de 
sementes de espécies florestais (BRASIL, 2013). Para realização do teste inicialmente 
a dormência das sementes foi superada com água a 80ºC por 24 horas. Para o teste 
foram utilizadas caixas de plástico transparente do tipo gerbox, com tampa e duas folhas 
de papel germitest como substrato, umedecidas com 2,5 vezes a massa do papel com 
água destilada, o teste foi incubado em câmara de germinação à temperatura constante 
de 25ºC. As avaliações do número de sementes germinadas foram realizadas diariamente, 
possuindo como critério de germinação a emissão da raiz primária. Ao concluir o teste, 
que teve duração de 10 dias, foram determinados: porcentagem de germinação (%G), 
considerando a relação percentual entre o número de sementes com emissão de raiz 
primária e o número de sementes colocadas para germinar e o índice de velocidade de 
germinação (IVG), resultante do somatório da razão da germinação diária pelo tempo, em 
dias, decorrido do início do teste. Juntamente com as avaliações destes testes realizou-se 
a avaliação do comprimento de plântula (CMP), a qual foi obtida com o auxílio de régua 
graduada em centímetros. As variáveis G(%), IVG e CMP, foram submetidas a Análise 
de variância e as médias comparadas pelo teste de Tukey a 5% de probabilidade. Para 
realização das análises foi utilizado o software GENES (CRUZ, 2013). 


Teste de Condutividade Elétrica 


Para realização do teste de condutividade elétrica (CE) foram utilizadas 4 repetições 
de 25 sementes, as quais foram pesadas em balança analítica e postas em recipiente 
contendo 75 mL de água destilada, colocados em câmara germinadora tipo B.O.D. com 
temperatura constante de 25 ºC durante um período de embebição de 24 horas. Após esse 
período de embebição, realizou-se a leitura da condutividade elétrica da solução na qual as 
sementes ficaram imersas, utilizando um condutivímetro digital de bancada (MS Tecnopon, 
modelo mCA 150). Os resultados de leitura foram divididos pelos respectivos valores de 
massa das amostras das sementes, sendo os valores expressos em uS.cm.g de semente. 
A variável CE foi submetida a Análise de variância e as médias comparadas pelo teste de 
Tukey a 5% de probabilidade. Para realização da análise foi utilizado o software GENES 
(CRUZ, 2013). 


Conceitos e Conhecimentos de Métodos e Técnicas de Pesquisa Científica em Engenharia E 
Florestal 2 Capítulo 1 


Repetibilidade e Dissimilaridade Genética de Características Biométricas 


Para estimar o coeficiente de repetibilidade foram avaliadas 100 sementes de cada 
matriz, para determinar as características: peso de sementes (PS), diâmetro longitudinal da 
semente (DLS), diâmetro equatorial da semente (DES) e espessura da semente (ES). O 
PS (9) foi determinado em balança analítica com 0,0001 g de precisão. Os DLS, DES e ES 
(mm) foram obtidos com o auxílio de paquímetro digital com precisão de 0,01 mm. 

A estimativa dos coeficientes de repetibilidade foi realizada por meio da análise 
de variância de fator único (ANOVA); componentes principais com base na matriz de 
correlações (CPC) e de covariânicas (CPCV); e análise estrutural (AE) com base na matriz 
de covariância. 

No método de análise de variância o coeficiente de repetibilidade é estimado por 
meio de resultados da própria variância, que é calculada de acordo com a equação: 


Yi; = u + gi + a; + Eij 


Em que: Y,: observação referente ao i-ésimo indivíduo na j-ésima medição; uy: média 
geral; g;: efeito aleatório da i-ésima subamostra sob influência do ambiente permanente 
(i=1,2,... n indivíduos); a; efeito da j-ésima medição (j = 1, 2, 3); e; erro experimental 
associado à observação Y,. 

O coeficiente de repetibilidade mede a maior ou menor capacidade das plantas de 
repetir a expressão fenotípica de uma determinada característica, e é obtido pela equação: 


2 
Ig 





2 2 
“+ 0 


Em que: o? = estimativa da variância entre genótipos e o? = estimativa da variância 
do erro experimental. 

O coeficiente de determinação, que representa a porcentagem de certeza na 
predição do valor real dos indivíduos selecionados com base em n medições é obtido pela 
equação: 


- nr 
“a+r(ng-1) 


R? 
Após estimado o coeficiente de repetibilidade (r), foi realizada a estimativa do 
número de medições (n,) necessárias para predizer o valor real dos indivíduos em diferentes 
porcentagens de determinação (80, 85, 90, 95 e 99%), obtido pela equação: 
RE) 
> q-RJA 
A dissimilaridade genética foi avaliada para as 10 matrizes a fim de obter o 
agrupamento das mesmas por similaridade. O método de agrupamento utilizado foi 
o UPGMA (Unweighted Pair Group Method Using Arithmetic Average) obtendo-se o 
dendrograma pela distância generalizada de Mahalanobis (D?). A contribuição relativa de 
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cada variável para a divergência foi avaliada pelo método de Singh (SINGH, 1981). As 
análises estatísticas foram realizadas por meio do software GENES (CRUZ, 2013). 


RESULTADOS E DISCUSSÃO 


As sementes de M. scabrella apresentaram após o beneficiamento, teor de água 
médio de 6,88 +1,08%, e peso de mil sementes de 16,46 + 1,90 g (Tabela 1). A umidade 
na semente proporciona condições diversas no armazenamento, portanto quanto maior o 
teor de água da semente armazenada, mais elevado será o número de fatores adversos à 
conservação da sua qualidade fisiológica, e no presente trabalho o teor de água encontrado é 
considerado favorável ao armazenamento em embalagens impermeáveis (MAZUCHOWSKI 
et al., 2014). Conforme as Regras para Análise de Sementes (BRASIL, 2009), as sementes 
de bracatinga (M. scabrella) são classificadas como sementes pequenas. 

Ressalta-se que para algumas matrizes não foram obtidas sementes em quantidades 
suficientes para avaliação do peso de mil sementes e outros parâmetros (Tabela 1), o que 
pode ser devido ao processo de maturação da espécie variar entre plantas, na mesma região 
e no mesmo ano, ou seja, diferentes matrizes apresentam maturação de frutos/sementes 
em períodos diferentes, o que pode ter ocasionado a redução do número de sementes em 
algumas matrizes, no momento da operação de colheita (MAZUCHOWSKI et al., 2014). 











Matriz Teor de água (%) Eee 
1 6,32 18,03 
2 7,25 15,35 
3 8,85 15,16 
4 6,16 14,23 
5 5,67 17,76 
6 7,14 + 
7 7,23 19,37 
8 6,66 -* 
9 8,15 -* 
10 5,85 15,33 
Média 6,88 16,46 
Desvio Padrão 1,08 1,90 
CV (%) 15,67 11,54 





Tabela 1 — Teor de água (%) e peso de mil sementes (9) para sementes de dez matrizes de Mimosa 
scabrella Benth., colhidas em Curitibanos — SC. 


*Matrizes com sementes insuficientes para determinação deste parâmetro. 


A análise de variância constatou diferenças significativas (Tabela 2), para todas as 
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características avaliadas, germinação (G%), índice de velocidade de germinação (IVG), 
comprimento de plântula (CMP) e condutividade elétrica (CE). 

Verifica-se que no teste de germinação se destacaram as matrizes 3 e 9 como de 
menor qualidade fisiológica, diferindo estatisticamente das demais, ou seja, as matrizes 
1,2,4,5, 7 e 10 foram agrupadas num mesmo nível de viabilidade, não diferindo entre 
si (Tabela 2). Os resultados obtidos no teste de germinação corroboram com Marcos 
Filho (2015), uma vez que, que os objetivos dos testes de vigor são identificar diferenças 
significativas na qualidade fisiológica de lotes de sementes com germinação similar. 
Com base no padrão apresentado para germinação da espécie M. scabrella é de 71% 
(WIELEWICKI et al., 2006), verifica-se que as sementes de todas as matrizes avaliadas 
estão acima do padrão. 





Matrizes G (%) IvG CMP (cm) CE (uS cm” 97) 

1 100,00 a 55,23 ab 16,26 a 27,97 ab 

2 99,25 a 55,26 ab 15,58 ab 31,65 bed 

3 83,25 b 48,46 bc 13,21 c 33,61 d 

4 98,00 a 54,37 ab 13,80 bc 34,30 d 

5 100,00 a 59,65 a 14,49 abc 27,42 a 

7 98,25 a 583,17 ab 14,98 abc 26,39 a 

8 E” E” E” =” 

9 83,00 b 40,27 c 14,79 abc 29,14 abc 
10 97,00 a 59,27 a 15,03 ab 32,46 de 


Tabela 2 — Médias de germinação (G), índice de velocidade de germinação (IVG), comprimento de 
plântula (CMP) e condutividade elétrica (CE), para sementes de dez matrizes de Mimosa scabrella 
Benth, colhidas em Curitibanos — SC. 


* Médias seguidas pela mesma letra na coluna, não diferem estatisticamente entre si, pelo teste de 
Tukey, a 5% de probabilidade. 


**Matrizes com sementes insuficientes para determinação deste parâmetro. 


Menegatti (2015), em estudo com M. scabrella, observou que o teste de germinação 
é de grande importância para explicar as variações existentes entre procedências distintas 
de ocorrência natural. Conforme Alves et al. (2005), em estudo sobre a germinação de 
sementes de M. caesalpiniifolia Benth, oriundas de diferentes procedências, estas não 
foram influenciadas pelo tamanho das sementes, demonstrando que pode não haver 
correlação entre os caracteres biométricos de sementes e as variáveis do processo 
germinativo, explicando a não adoção de classes de tamanho como indicativo de sucesso 
no estabelecimento de plântulas. 

Ao final do teste de germinação, foi avaliado o IVG, sendo, os maiores IVG 
encontrados para as matrizes 5 e 10, que não diferiram estatisticamente, das matrizes 1, 
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2,4 e 7, respectivamente. Menegatti (2015), considerou a variável IVG adequada para a 
diferenciação de matrizes. Diversos autores relacionam os dados de IVG para distinguir os 
melhores tratamentos e/ou lotes de sementes, uma vez que, quanto maior o valor de IVG, 
maior é a germinação diária, sendo estes lotes os de maior vigor (GOMES, 2013). 

A determinação do comprimento médio das plântulas normais é realizada, 
considerando que as amostras que expressam os maiores valores são mais vigorosas 
(NAKAGAWA, 1999). As médias do comprimento de plântula variaram de 13,21 a 16,26 
cm, dentre as matrizes, sendo que as matrizes 1, 2,3, 5, 7, 9 e 10 apresentaram maior 
qualidade fisiológica, diferindo-se estatisticamente das matrizes 3 e 4, demonstrando que 
este parâmetro pode ser utilizado para diferenciar matrizes quanto ao vigor. Em estudos 
com M. scabrella, Menegatti (2015), encontrou valores variando de 10,44 a 14,85 cm, 
inferindo que esta variável se mostra eficiente na detecção das diferenças e separação das 
matrizes, quanto à qualidade fisiológica das sementes. 

Com relação ao teste de condutividade elétrica, a análise dos dados mostrou 
relação direta com os testes anteriores na avaliação da qualidade das sementes das 
matrizes. Destacando-se as matrizes 5 e 7 como as de melhor qualidade e as matrizes 
3 e 4 com o valor elevado de condutividade elétrica, sendo, portanto, as que apresentam 
menor qualidade fisiológica. À medida que a semente vai envelhecendo, ocorre a sua 
deterioração, resultando na perda da integridade dos sistemas de membranas da célula, 
resultando no aumento de sua permeabilidade, o que causa a lixiviação de eletrólitos e 
maior condutividade elétrica (SANTOS e PAULA, 2009). 

Na análise de repetibilidade (Tabela 4) houve diferenças significativas entre as 
matrizes avaliadas, para todas as características de acordo com a análise de variância, 
isto também foi verificado por Santos et al. (2010), e demonstra que o componente de 
variância genético, confundido com os efeitos do ambiente é significativo nessa população. 
As estimativas dos coeficientes de repetibilidade (r) obtidos com relação as quatro 
características por meio dos quatro métodos estatísticos utilizados. Observa-se que 
o melhor método para obtenção das estimativas do coeficiente de repetibilidade é dos 
componentes principais de covariância (CPCV), isso ocorre pelo fato de este considerar o 
comportamento cíclico dos caracteres (ABEYWARDENA, 1972). 

Os resultados obtidos nos diferentes métodos seguem o padrão de resultados 
observados por Farias (2016), ao analisar frutos de Passiflora spp., no qual as estimativas 
do coeficiente de repetibilidade obtidas pelo método da análise de variância (ANOVA) 
foram sempre inferiores ou iguais as estimativas obtidas pelos demais métodos. Conforme 
Costa (2003), os valores inferiores obtidos por meio do método da ANOVA, são resultantes 
da variância genotípica, empregada para estimar a repetibilidade, já que esta pode não ser 
completamente de origem genética, uma vez que o componente de variância do ambiente 
entre indivíduos permanece somado com a variância genotípica, enquanto o método de 


componentes principais permite isolar o efeito da alternância. 
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Característica Coeficiente ANOVA CPCV CPC AE 





SE r 0,03 0,46 0,28 0,03 
R2 76,95 98,85 97,55 76,95 
r 0,56 0,61 0,59 0,56 
DLS 
R2 99,22 99,37 99,33 99,22 
r 0,21 0,29 0,29 0,21 
DES 
R2 96,47 97,71 97,63 96,47 
Eé r 0,22 0,33 0,29 0,22 
R2 96,68 98,08 97,69 96,68 





Tabela 4 — Estimativa dos coeficientes de repetibilidade (r) e dos coeficientes de determinação (R2) 
utilizando os métodos de análise de variância (ANOVA), dos componentes principais baseado na 
matriz de correlação (CPC) e de covariância (CPCV), e análise estrutural (AE) baseado na matriz de 
covariância. 


Legenda: PS: peso de sementes, DLS: diâmetro longitudinal de sementes; DLS: diâmetro equatorial de 
sementes; ES: espessura de sementes. 


Conforme a classificação de Resende (2015), a característica que apresentou 
alta repetibilidade (r = 0,60) foi o DLS, para o método CPCV, apresentando repetibilidade 
intermediária para os demais métodos, com variação de 0,56 <r < 0,59. Para a característica 
de OS, obteve-se repetibilidade intermediária, com valor de r < 0,46, para o método CPCV. 
As demais características avaliadas DES e ES, apresentaram repetibilidade baixa para 
o método CPCV, com valores de 0,21 <r < 0,29 e 0,22 <r < 0,33 respectivamente, para 
coeficiente de determinação, acima de 97%. 

Oliveira e Moura (2008), por meio dos métodos de análise do coeficiente de 
repetibilidade é possível determinar um número mínimo de medições necessárias por 
indivíduo, observando quais características apresentam maiores contribuições, para que 
se possa considerar no processo de melhoramento, para os diferentes métodos pode ser 
visualizado na tabela 5. 





Coeficiente de determinação (%) 





Caract. Método 
80 85 90 95 99 
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ANOVA 119,7(120) 169,7(170) 269,5(270)  569,0(569) 2965,0(2965) 





= CPCV 4,6(5) 6,5(7) 10,4(11) 21,9(22) 114,5(115) 
CPC 10,0(10) 14,2(15) 22,5(23) 47,6(48) 248,2(249) 
AE 16,4(17) 23,2(24) 36,9(37) 77,9(78) 406,0(406) 
ANOVA 3,1(4) 4,4(5) 7,0(7) 14,8(15) 77,3(78) 
suê CPCV 2,5(3) 3,5(4) 5,6(6) 11,8(12) 61,8(62) 
CPC 2,6(3) 3,7(4) 6,0(6) 12,7(13) 66,2(67) 
AE 2,9(3) 4,2(5) 6,6(7) 14,1(15) 73,5(74) 
ANOVA 14,6(15) 20,6(21) 32,8(33) 69,3(70) 361,3(362) 
Ee CPCV 9,3(10) 13,2(14) 21,0(21) 44,5(45) 231,9(232) 
CPC 9,7(10) 13,7(14) 21,8(22) 46,0(47) 240,1(241) 
AE 13,6(14) 19,3(20) 30,7(31) 64,9(65) 338,3(339) 
ANOVA 13,7(14) 19,4(20) 30,8(31) 65,1(66) 339,6(340) 
CPCV 7,7(8) 11,0(11) 17,5(18) 37,0(37) 192,8(1983) 
Ee CPC 9,4(10) 13,3(14) 21,2(22) 44,8(45) 233,6(234) 
AE 14,3(15) 20,3(21) 32,3(33) 68,2(69) 355,7(356) 


Tabela 5 — Número de medições necessárias para obtenção de diferentes coeficientes de determinação 
para as características peso de semente (PS), diâmetro longitudinal de semente (DLS), diâmetro 
equatorial de semente (DES), e espessura de semente (ES), para cada um dos métodos testados. 


* Valores em parênteses indicam valores arredondados de número de sementes. 


Legenda: ANOVA: Análise de variância; CPC: componentes principais baseados na matriz de 
correlação, CPCV: componentes principais baseados na matriz de covariância; AE: análise estrutural 
baseada na matriz de covariância. 


Conforme Cruz et al. (2012), o número mínimo de medições necessárias é 
inversamente proporcional ao coeficiente de repetibilidade, pois, quanto maior a estimativa 
do coeficiente, menor será o número de medições necessárias. Isto justifica os números 
mínimos de medições necessárias encontrados para CPCV serem iguais ou inferiores aos 
outros métodos para todas as características em todos os graus de liberdade. 

Farias Neto et al. (2003) considera como satisfatório para a tomada de decisão sobre 
a superioridade dos indivíduos ter como base um coeficiente de determinação de 90%. 
Assim, considerando-se o método CPCV, o número mínimo de medições para alcançar 
90% de determinação é de 11 sementes para PS, 6 sementes para DLS, 21 sementes para 
DES e 18 sementes para ES. De modo geral, indica-se a medição de 21 sementes. Desse 
modo, pode-se afirmar que um R? de 90% é suficiente para garantir, além do mínimo gasto 
dos recursos e mão de obra para a realização da caracterização, a confiabilidade mínima 
dos dados. 

Além da estimativa do coeficiente de repetibilidade para as características de 
interesse, é de extrema importância mensurar a contribuição relativa que os caracteres 
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desempenham na divergência genética entre matrizes (ROVERI NETO, 2014). Dessa forma, 
com base na metodologia proposta por Singh (1981), com relação a contribuição relativa 
dos caracteres para diversidade, pode-se observar na tabela 6 que as características de 
DLS e ES são as que apresentam maior contribuição, entre as características avaliadas, 
sendo estas responsáveis por 77,67 e 12,74%, respectivamente, de toda a variação para a 
diversidade encontrada entre as matrizes estudadas. 








Característica Contribuição Relativa (%) 
os 0,64 
DLS 77,67 
DES 8,95 
ES 12,74 





Tabela 6 — Contribuição relativa das características morfológicas em sementes de Mimosa scabrella 
Benth. para a variabilidade entre matrizes 


Legenda: PS = Peso de sementes; DLS = Diâmetro Longitudinal de Sementes; DES = Diâmetro 
Equatorial de Sementes e ES = Espessura de Sementes. 


Menegatti (2015), em estudo avaliando a dissimilaridade para sementes de M. 
scabrella, verificou que os caracteres que mais contribuíram para a divergência foram o 
comprimento e largura de sementes, porcentagem de plântulas normais, porcentagem 
de germinação, índice de velocidade de germinação e massa verde de plântulas. Dessa 
forma características relacionadas ao teste de germinação, como índice de velocidade de 
germinação e percentual de germinação, foram mais relevantes para explicar as variações 
existentes entres as diferentes procedências, porém, neste trabalho características 
fisiológicas não foram avaliadas, impedindo assim uma comparação direta. 

Flôres Júnior et al. (2018), em estudo com sementes de Acacia mearnsii, também 
observou que a característica peso de sementes apresentou menor contribuição para as 
análises de diversidade, porém, a retirada de tais características com baixas contribuições 
pode acarretar em alteração do agrupamento inicial, não sendo recomendadas, 
principalmente quando não se realizaram várias análises, em diferentes anos. 

As análises das características biométricas das sementes permitiram a diferenciação 
das dez matrizes, bem como seu agrupamento de acordo com a similaridade (Figura 1). As 
plantas desse estudo se dividem em dois grupos, sendo um representado pelas matrizes 
4,10,3,2,e 1 e ouro pelas matrizes 6, 8, 9, 5 e 7. Ainda que o método UPGMA favoreça a 
separação de indivíduos, como enfatizado por Silva (2012), nesse estudo não foi verificada 
a presença de indivíduos de outras populações, o que é indicado pela ausência de grupos 
individualizados (Figura 1). 
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MATRIZ 04 
MATRIZ 10 
MATRIZ 03 
MATRIZ 02 
MATRIZ 01 
MATRIZ 06 
MATRIZ 08 
MATRIZ 09 
MATRIZ 05 
MATRIZ 07 


Figura 1 — Dendrograma com agrupamento UPGMA, através da distância genética de Mahalanobis 
para 10 matrizes de Mimosa scabrella Benth. com base nas características biométricas das sementes 
colhidas em 2018. 


As matrizes 4 e 10 foram as que apresentaram maior similaridade entre si. De 
modo geral pode-se observar que a distância entre as matrizes foi baixa, isso se deve ao 
fato de essas matrizes pertencerem à mesma população, indicando que as mesmas são 
aparentadas. 

Conforme Cruz et al. (2012), o cruzamento entre famílias de grupos similares 
deve ser evitado para garantir a variabilidade, pois a mesma é fundamental em programa 
de melhoramento e conservação genética. Desse modo, as matrizes 4 e 7 são as mais 
divergentes e devem compor programas de intercruzamento para obtenção de genótipos 
superiores, em futuros programas de melhoramento para a espécie em estudo. 


CONCLUSÃO 


Dentre as matrizes avaliadas por meio da análise de qualidade fisiológica, para a 
formação de um lote de sementes de alta qualidade indica-se apenas as matrizes 1, 2,4,5, 
7 e 10. O coeficiente de repetibilidade foi alto para a característica DLS, intermediário para o 
PS, e, baixo para DES e ES. É necessária a medição de 21 sementes para garantir 90% de 
determinação. Verificou-se formação de dois grupos através na análise de dissimilaridade 
genética. 
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RESUMO: Objetivou-se avaliar o crescimento 
em diâmetro do colo e altura de mudas de Cordia 
trichotoma produzidas a partir de diferentes 
tipos de propágulo e submetidas a diferentes 
adubações de cobertura e tempo de avaliação. O 
experimento foi conduzido em casa de vegetação 
e de sombra com miniestacas enraizadas 
cultivadas em tubetes de 110 cm* preenchidos 
com substrato comercial e vermiculita na 
mesma proporção. Testaram-se três tipos de 
propágulo (uma gema, duas gemas e apical), 
quatro adubações de cobertura (diferentes 
concentrações de nitrogênio e potássio), em 
quatro avaliações (30, 60, 90 e 120 dias). Mudas 
de C. trichotoma produzidas por miniestaquia 
e adubadas com sulfato de amônio e cloreto 
de potássio apresentam bom crescimento em 
diâmetro do colo e altura. 

PALAVRAS - CHAVE: Espécie nativa; adubação 
de cobertura; tipo de propágulo; produção de 
mudas. 


GROWTH OF SEEDLINGS OF CORDIA 
TRICHOTOMA SUBMITTED TO 
DIFFERENT TREATMENTS IN NURSERY 
ABSTRACT: The objective was to evaluate the 
growth in diameter and height of Cordia trichotoma 
seedlings in different type of propagule, cover 


Capítulo 2 


fertilizations and time of evaluation. The experiment was carried out in a greenhouse and 
shade house with rooted mini-cuttings from 110 cm* tubes filled with the same proportions 
of commercial substrate and vermiculite. The treatments were: three types of propagule 
(one bud, two buds and apical), four cover fertilizations (different concentrations of nitrogen 
and potassium), and four evaluations (30, 60, 90 and 120 days). Plantlets of C. trichotoma 
produced by mini-cuttings and fertilized with ammonium sulfate and potassium chloride show 
good growth in diameter and height. 

KEYWORDS: Native species; cover fertilization; propagule type; plantlet production. 


11 INTRODUÇÃO 


Cordia trichotoma (Vell.) Arrab. ex Steud., conhecida popularmente como louro- 
pardo, é uma espécie de porte arbóreo nativa do Brasil pertencente à família Boraginaceae. 
Integra o grupo sucessional secundário inicial com tendência a ser pioneira, podendo atingir 
até 35 m de altura e diâmetro à altura do peito (DAP) entre 70 e 90 cm. O louro-pardo é uma 
planta caducifólia na estação de repouso vegetativo (CARVALHO, 2003; LORENZI, 2008). 

A espécie apresenta potencial para o plantio com fins econômicos, devido a sua 
grande durabilidade e qualidade da madeira, sendo empregada na confecção de móveis, 
entre outros. Entretanto, para se estabelecer plantios de qualidade é indispensável que 
se tenha mudas com excelentes caracteres silviculturais. Nesse sentido, a propagação 
vegetativa através da miniestaquia permite propagar genótipos de interesse. 

A propagação vegetativa de plantas consiste na utilização de partes de uma planta 
doadora, denominados propágulos, destinados à produção de novos indivíduos completos 
e idênticos à planta que deu origem ao material (SOUZA JUNIOR; WENDLING, 2003). A 
miniestaquia visa a formação de mudas por meio do enraizamento de um propágulo, o qual, 
por sua vez, consiste em um processo complexo que envolve o redirecionamento de células 
para formação de um novo sistema radicular adventício. Por ser um processo complexo, 
depende de diferentes fatores, como por exemplo, o tipo de propágulo, o substrato, a 
concentração de ácido indolbutírico (AIB), entre outros. 

O crescimento de mudas é também uma fase importante, uma vez que as plantas 
mais vigorosas permitem maior chance de sucesso no estabelecimento de povoamentos 
florestais (BARBOSA et al., 2003). Vários aspectos também afetam a qualidade da muda, 
entre eles, o tipo de recipiente, substrato, a adubação, irrigação e o manejo em viveiro. 
Em razão do grande número de espécies, e de suas exigências nutricionais, não há 
recomendações específicas de adubação (GONÇALVES et al., 2005). 

Gonçalves et al. (2005) recomendam repetir a adubação de cobertura em intervalos 
de 7 a 10 dias para as espécies pioneiras e secundárias iniciais, como é o caso do louro- 
pardo, cuja recomendação é de 200 g de sulfato de amônio [(NH,))SO )Je 150 g de cloreto de 
potássio (KCl), em 100 litros de água, com aplicações intercaladas de cloreto de potássio. 
Já Valeri e Corradini (2005) indicam a fertilização de mudas de Eucalyptus e Pinus, com 
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1.000 g de sulfato de amônio e 300 g de cloreto de potássio, dissolvidos em 100 litros de 
água, a cada 7 ou 10 dias, alternando a aplicação de cloreto de potássio. 

Apesar do louro-pardo ser uma espécie com aptidão madeireira, ainda são escassos 
os estudos que abordam a produção de mudas por meio de propagação vegetativa. Assim, 
objetivou-se avaliar o crescimento em diâmetro do colo e a altura de mudas de Cordia 
trichotoma produzidas a partir de diferentes tipos de propágulo e submetidas a diferentes 
adubações de cobertura e tempo de avaliação. 


21 METODOLOGIA 


Este estudo foi conduzido durante os meses de março a setembro de 2020 no Núcleo 
de Melhoramento e Propagação Vegetativa de Plantas (MPVP), da Universidade Federal 
de Santa Maria. As mudas utilizadas para estudar o crescimento em diâmetro do colo e 
altura foram originadas pelo enraizamento de miniestacas coletadas em minijardim clonal. 

O minijardim clonal foi construído em bandejas de polietileno com sistema fechado 
de cultivo e areia grossa como substrato (BISOGNIN et al., 2015). As minicepas foram 
fertirrigadas de 2 a 3 vezes ao dia, por 15 min., com o auxílio de um programador digital e 
bomba submersa de baixa vazão. A solução nutritiva foi preparada conforme descrito por 
Pimentel et al. (2020). 

As brotações das minicepas foram coletadas e seccionadas em miniestacas com 
uma gema (GU), duas gemas (2G), e a porção apical (AP) com área foliar original foi 
reduzida em 50%. As miniestacas tiveram sua porção basal imersa na solução hidroalcoólica 
de ácido indolbutírico com concentração de 3.000 mg L* por 10 s. O enraizamento das 
miniestacas foi realizado em câmara úmida automatizada com temperatura média de 32º C 
e nebulização a cada 30 min., por 10 s., entre as 7 às 19 horas. 

As miniestacas enraizadas foram transferidas para a casa de vegetação com 25º 
C de temperatura média, para as quais foi fornecida irrigação manual, uma vez ao dia por 
30 dias. Decorrido esse período, as mudas foram transferidas para a casa de sombra com 
tela de 50% de permeabilidade dos raios solares e equipada com irrigação automatizada 2 
vezes ao dia, durante 15 min. 

O experimento foi um fatorial (3 x 4 x 4) no delineamento experimental inteiramente 
casualizado sendo testado três tipos de propágulo (uma gema, duas gemas e porção 
apical), quatro adubações de cobertura (concentrações de sulfato de amônio - (NH,)SO, - e 
cloreto de potássio - KCl) em quatro avaliação (30, 60, 90 e 120 dias). As concentrações 
do (NH,)SO, e do KCl foram, respectivamente, de: A1 = 200 e 150 9; A2 = 1.000 e 300 g; 
A3 = 750 e 2259; A4 = 1.250 e 3759. Para a fertilização Al, intercalou-se a aplicação de 
cloreto de potássio, a cada 10 dias e a fertilização dos tratamentos A2, A3 e A4 ocorreu a 
cada 30 dias. 

O diâmetro do colo e a altura foram mensuradas mensalmente até 120 dias com 
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auxílio de um paquímetro eletrônico graduado em milímetros e uma régua graduada em 
centímetro, respectivamente. Os dados foram submetidos à análise de covariância através 
do procedimento GLM do pacote estatístico SAS OnDemand for Academics (SAS, 2014). 
Quando constatada a existência de diferenças significativas em ao menos um dos fatores 
analisados e/ou suas interações, estes foram comparados pelo teste de Tukey (a=0,05). 


31 RESULTADOS E DISCUSSÃO 


O conjunto de dados utilizado no presente estudo foi composto por 422 pares de 
observações de crescimento inicial das mudas de louro-pardo. O resumo estatístico das 
variáveis analisadas no presente estudo (valores médios, mínimos e máximos, bem como 
o desvio padrão da variável diâmetro do colo e altura) é apresentado na Tabela 1. 














Ea Média Desvio padrão Mínimo Máximo 
Diâmetro do colo (mm) 
30 2,32 + 0,59 1,02 4,17 
60 217 +0,53 0,85 3,74 
90 2,40 + 0,57 1,31 3,95 
120 2,40 +0,61 1,08 3,91 
Altura (cm) 
30 1,95 + 1,66 0,20 8,00 
60 2,21 + 1,64 0,40 9,00 
90 2,53 + 1,78 0,20 9,50 
120 3,14 + 1,96 0,50 9,50 





Tabela 1 - Valores médios, mínimos, máximos e desvio padrão do diâmetro do colo e da altura de 
mudas de Cordia trichotoma em diferentes períodos de tempo. 


Analisando os quatro períodos de avaliação, é possível verificar que o diâmetro 
do colo foi menor aos 60 dias, quando comparado com o de 30 dias, o que ser um 
resultado da troca de ambiente neste período. É preciso salientar que inicialmente as 
mudas permaneceram em casa de vegetação com condições ambientais controladas de 
temperatura e luminosidade, sendo posteriormente transferidas para casa de sombra, 
que consiste em um ambiente sem o mesmo controle ambiental. Essa transferência pode 
resultar na morte de algumas mudas e consequentemente, na redução no número de 
observações neste intervalo de tempo. Contudo, é possível observar que aos 90 e 120 dias 
as mudas retomam o crescimento em diâmetro do colo, indicando condição de adaptação 
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e estabelecimento. Este comportamento não foi observado para a variável altura, que 
aparentemente não sofreu com a mudança de ambiente resultando em desenvolvimento 
contínuo. 


3.1 Anacova para Crescimento em Diâmetro do Colo 


O crescimento em diâmetro do colo das mudas, em um primeiro momento, foi 
verificado conforme uma análise de covariância, contemplando todos os fatores analisados 
(propágulo, adubação de cobertura e tempo) e suas interações (Tabela 2). 





Fonte de variação G.L. S.Q. Q.mM. F Prob.>F 
Modelo 29 43,5 1,5 6,1 < 0,0001 
Propágulo 2 27,9 13,9 56,3 < 0,0001 
Adubação de cobertura 3 2,5 0,8 3,3 0,0202 

Tempo 3 3,8 1,3 5,1 0,0018 

Propágulo vs. Adubação de cobertura 6 9,6 1,6 6,5 < 0,0001 
Adubação de cobertura vs. Tempo 9 2,5 0,3 1,1 0,3518ns 
Tempo vs. Propágulo 6 0,7 0,1 0,5 0,8313ns 
Erro 392 97,0 0,2 

Total corrigido 421 140,5 





Tabela 2 - Análise de covariância para o crescimento inicial em diâmetro do colo de mudas de Cordia 
trichotoma. 


Em que: G.L = graus de liberdade; S.Q. = soma de quadrados; Q.M. = quadrado médio; F = valor de F 
calculado e Prob. > F = probabilidade de significância. "s = não significativo. 


Analisando a Tabela 2 é possível verificar que o modelo de regressão ajustado foi 
significativo, com valor de F calculado igual 6,1 e probabilidade de significância menor que 
0,0001. Para essa variável, ao menos um dos tratamentos diferiu estatisticamente dos 
demais, denotado pela probabilidade de significância. Analisando as interações, verifica- 
se que houve interação entre a adubação e tipo de propágulo, com probabilidade menor 
que 0,0001. Isto sugere que o crescimento em diâmetro do colo das mudas de louro-pardo 
é afetado em relação ao tipo de propágulo em conjunto com a adubação de cobertura 
aplicada. Por outro lado, entre a adubação de cobertura e o tempo das avaliações, bem 
como para o tempo de avaliação e o tipo de propágulo não foi observada interação 
significativa. Isto indica que os tipos de propágulo e adubação de cobertura com o decorrer 
do tempo, tendem a não influenciar no crescimento em diâmetro do colo. 

Com o intuito de ajustar um modelo contendo apenas as variáveis e interações 
significativas para o diâmetro do colo, calculou-se novamente a análise covariância, porém, 
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com a exclusão das interações adubação de cobertura vs. Tempo, e Tempo vs. Propágulo 
(Tabela 3). O modelo ajustado apresentou valor de F calculado igual a 11,7 significativo 
com probabilidade menor que 0,0001. Este modelo de regressão apresentou coeficiente 
de determinação igual a 0,2874, erro padrão da estimativa de 0,4960 mm e coeficiente de 
variação igual a 21,38%. 





Fonte de variação G.L. S.Q. Q.mM. F Prob.>F 
Modelo 14 40,4 2,9 11,7 < 0,0001 
Propágulo 2 27,3 13,6 55,5 < 0,0001 
Adubação de cobertura 3 2,2 0,7 3,0 0,0295 
Tempo 3 4,2 1,4 5,6 0,0009 
Propágulo vs. Adubação de cobertura 6 9,5 1,6 6,5 < 0,0001 
Propágulo vs. Adubação de cobertura 407 100,1 0,2 

Total corrigido 421 140,5 





Tabela 3 - Análise de covariância para o diâmetro do colo de mudas de Cordia trichotoma com a 
exclusão das interações não significativas. 


Em que: G.L = graus de liberdade; S.Q. = soma de quadrados; Q.M. = quadrado médio; F = valor de F 
calculado e Prob. > F = probabilidade de significância. 


3.2 Anacova para Crescimento em Altura 


Aaltura das mudas foi verificada conforme uma análise de covariância, contemplando 
todas as fontes de variação analisadas e suas interações (Tabela 4). Foi possível verificar 
que o modelo de regressão ajustado foi significativo. Assim como para o diâmetro do 
colo, todas as fontes de variação apresentaram diferença significativa isoladamente. Há 
diferença entre os tipos de propágulo, entre as diferentes adubações de cobertura, bem 
como uma diferença no ritmo de crescimento em altura das mudas. Além disso, o tipo de 
propágulo apresentou interação significativa com a adubação. 








Fonte de variação G.L. S.Q. Q.mM. F Prob.>F 
Modelo 29 415,6 14,3 6,0 < 0,0001 
Adubação de cobertura 3 88,0 29,3 12,3 < 0,0001 
Tempo 3 59,9 20,0 8,4 < 0,0001 
Propágulo 2 106,2 58,1 22,3 < 0,0001 
Propágulo vs. Adubação de cobertura 6 46,5 Fora 3,3 0,0039 
Adubação de cobertura vs. Tempo 9 12,0 1,3 0,6 0,8279"s 
Tempo vs. Propágulo 6 7,8 1,2 0,5 0,8001"s 
Erro 392 930,5 2,4 

Total corrigido 421 1346,1 





Tabela 4 - Análise de covariância para o crescimento inicial em altura de mudas de Cordia trichotoma. 
Em que: G.L = graus de liberdade; S.Q. = soma de quadrados; Q.M. = quadrado médio; F = valor de F 
calculado e Prob. > F = probabilidade de significância. "s = não significativo. 
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Com o intuito de ajustar um modelo contendo apenas as variáveis e interações 
significativas para a altura das mudas de louro-pardo, calculou-se novamente a análise 
covariância, porém, com a exclusão das interações adubação de cobertura vs. Tempo, 
e Tempo vs. Propágulo (Tabela 5). O modelo ajustado apresentou valor de F calculado 
igual a 12,0 significativo com probabilidade menor que 0,0001. Este modelo de regressão 
apresentou coeficiente de determinação igual a 0,2917, erro padrão da estimativa de 
1,5323 cm e coeficiente de variação igual a 63,98%. 








Fonte de variação G.L. S.Q. Q.mM. F Prob.>F 
Modelo 14 392,8 28,1 12,0 < 0,0001 
Propágulo 2 101,4 50,7 21,6 < 0,0001 
Adubação de cobertura 3 81,3 271 11,6 < 0,0001 
Tempo 3 66,6 22,2 9,5 < 0,0001 
Propágulo vs. Adubação de cobertura 6 45,3 7,6 3,2 0,0042 

Propágulo vs. Adubação de cobertura 407 953,3 2,3 

Total corrigido 421 1346,1 





Tabela 5 - Análise de covariância para a altura de mudas de Cordia trichotoma com a exclusão das 
interações não significativas. 


Em que: G.L = graus de liberdade; S.Q. = soma de quadrados; Q.M. = quadrado médio; F = valor de F 
calculado e Prob. > F = probabilidade de significância. 


3.3 Teste de Média para o Tipo de Propágulo 


O maior crescimento em diâmetro do colo foi encontrado para as mudas originadas 
de propágulos de duas gemas (2G) e de uma gema (GU), com valor médio igual a 2,51 e 
2,40 mm, respectivamente (Tabela 6). Para diâmetro do colo o propágulo apical apresenta 
menor valor (1,75 mm). Para a altura das mudas não houve diferença entre o propágulo de 
duas gemas (2G) com 3,15 cm e apical (AP) com 2,62 cm, os quais diferiram da miniestaca 
constituída de uma gema que apresentou 1,82 cm. 








Tipo de propágulo Diâmetro do colo Altura 
(mm) (cm) 

2G 2,51 a* 3,15 a 

GU 240a 1,82b 

AP 1,75b 262a 





Tabela 6 - Teste de médias para o diâmetro do colo e altura de mudas de Cordia trichotoma para os 
diferentes tipos de propágulos. 


*Médias seguidas pela mesma letra não diferem estatisticamente pelo teste de Tukey (prob. = 0,05). 
Em que: 2G = 2 gemas; GU = uma gema e AP = gema apical. 
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O enraizamento adventício está relacionado com o gradiente hormonal do ramo, 
o status nutricional, o tamanho e o número de folhas mantidas na miniestaca podem 
afetar a competência ao enraizamento adventício, devido à produção de auxina endógena 
necessária para esse complexo processo (FACHINELLO et al., 2005; TAIZ; ZEIGER, 2017). 
Sendo assim o maior crescimento em diâmetro do colo e altura das mudas constituídas 
de duas gemas pode ser explicado em razão da maior competência ao enraizamento 
adventício e formação do sistema radicial. Devido as possíveis diferenças entre os tipos de 
miniestacas para o diâmetro e altura, a propagação vegetativa pode ser satisfatoriamente 
realizada a partir de miniestacas de uma ou duas gemas, e de apicais. 
3.4 Teste de Média para Adubação de Cobertura 


A maior média em diâmetro do colo foi obtida na fertilização A2 (1.000 g de 
sulfato de amônio e 300 g de cloreto de potássio aplicados a cada 30 dias), não diferindo 
estatisticamente das adubações A1 e A4 (Tabela 7). Esta adubação de cobertura foi 
superior a A3 (1.000 g de sulfato de amônio e 300 g de cloreto de potássio, a cada 15 dias), 
resultando em diferenças no crescimento de 10,2%. Para o crescimento em altura das 
mudas a adubação A3 (750 g de sulfato de amônio e 225 g de cloreto de potássio, a cada 
30 dias) e A1 (200 g de sulfato de amônio e 150 g de cloreto de potássio, a cada 30 dias) 
apresentaram maior valor, mas não diferiram do A2. Assim, o maior crescimento, tanto para 
o diâmetro do colo, quanto para altura foram observados nas adubações Ai e AZ. 








Adubação de Diâmetro do colo Altura 
cobertura! (mm) (cm) 

Ai 2,18 ab* 2,86 a* 

A2 2,35 a 2,40 ab 

A3 211b 3,02 a 

A4 2,24 ab 1,85b 





Tabela 7 - Teste de médias para o diâmetro do colo e altura de mudas de Cordia trichotoma nas 
diferentes adubações de cobertura. 


*Médias seguidas pela mesma letra não diferem estatisticamente pelo teste de Tukey (prob.=0,05). 
'Doses dos nutrientes sulfato de amônio [(NH )SO )] e cloreto de potássio (KCI), respectivamente, de: 
A1 =200e 1509, a cada 10 dias; A2 = 1.000 e 3009, a cada 30 dias; A3 = 750 e 2259, a cada 30 dias 
e A4 = 1.250 e 375 g de cloreto, a cada 30 dias. 


Considerando que no tratamento de adubação A2 (1.000 g de sulfato de amônio 
e 300 g de cloreto de potássio) o sulfato de amônio é fonte de nitrogênio (N) e enxofre 
(S) e o cloreto de potássio é o principal fornecedor de potássio (K) e cloro (Cl). Dentre 
esses nutrientes, o nitrogênio e o enxofre fazem parte dos compostos de carbono. O 
potássio é requerido como cofator e muitas enzimas, no estabelecimento do turgor das 
células, enquanto o cloro é requerido nas reações fotossintéticas (TAIZ; ZIEGLER, 2017). 
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Todos, exceto o cloro, são classificados como macronutrientes, ou seja, são requeridos 
em maior quantidade o que pode explicar os melhores resultados observados para os 
tratamentos com as maiores concentrações desses elementos. Já o cloro é classificado 
como micronutriente, mas também considerado um elemento essencial para o crescimento 
vegetal. 

Além disso, esse tratamento pode ser aplicado a cada 30 dias sem prejuízo ao 
crescimento, resultando em uma a redução da demanda de trabalho com essa atividade 
no viveiro, da quantidade de insumos, bem como os custos envolvidos nesta operação, 
quando comparado ao tratamento com aplicação a cada 10 dias. 

Inicialmente as formulações de adubações continham o nível mais alto possível, 
sem causar toxidade, contudo, deve-se ter preocupação com o uso racional dos recursos 
naturais em nossos estudos, pois a sugestão de fertilização a ser aplicada deve contemplar 
níveis que fornecem crescimento adequado, com mínimo uso de nutrientes, a fim de 
minimizar impactos ambientais. 


3.5 Teste de Média Para o Efeito do Tempo Sobre o Crescimento 


Observa-se, assim como era esperado, que a maior média de diâmetro de colo foi 
encontrada aos 120 dias, igual a 2,22 mm (Tabela 8). Este valor não difere estatisticamente 
dos mensurados aos 30 e 90 dias, contudo, foi superior a 60 dias. Para a altura a maior 
média também foi aos 120 dias, igual a 3,15 cm, diferindo dos demais tempos, exceto aos 
90 dias. 








Tempo (dias) das des colo atira 
30 2,22 ab* 2,08 €* 
em 2,06b 2.30 be 
sm 2,30 a 2,65 ab 
120 231a 3,15a 





Tabela 8 - Teste de médias para o diâmetro do colo e altura de mudas de Cordia trichotoma nos 
diferentes períodos de tempo avaliados. 


*Médias seguidas pela mesma letra não diferem estatisticamente pelo teste de Tukey (prob.=0,05). 


Para diâmetro do colo as médias para os 30, 60 e 90 dias não diferiram estatisticamente 
entre si, indicando que a troca de ambiente foi benéfica para o crescimento em diâmetro do 
colo das mudas, fornecendo as condições necessárias para o seu pleno desenvolvimento. 
O que pode ser ressaltado pela diferença entre o período de 60 e 120 dias. O período de 
avaliação das mudas deu-se nos meses de inverno, e de acordo com a Kôppen e Geiger a 
classificação do clima é “Cfa” (ALVARES et al., 2013), apresenta estações bem definidas, e 
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o mês de julho apresenta menor temperatura média, 13.1 ºC, e com menor radiação solar. A 
temperatura e irradiação solar são fatores importantes para o crescimento das plantas, pois 
exercem influência em todas as atividades fisiológicas, por controlar as taxas das reações 
metabólicas nas células. 


3.6 Efeito da Interação Entre a Adubação de Cobertura e o Tipo de Propágulo 


A interação entre a adubação de cobertura e o tipo de propágulo (Figura 1) demostra 
que a adubação composta por 1.000 de sulfato de amônio e 300 g de cloreto de potássio a 
cada 30 dias (A2) com a muda formada pela miniestaca de uma gema apresentou melhor 
resultado para o diâmetro do colo (2,72 mm), mas não diferiu das combinações A3 com uma 
gema (2,40 mm) e, da muda formada com miniestacas de duas gemas com as adubações 
de A1 (2,40 mm), A2 (2,61 mm), AS (2,57 mm) e A4 (2,38 mm). 








Diâmetro (mm) 








Al A2 A3 Aq 
Adubação de cobertura 








[Propágulo 26 EH AP E GU] 











Figura 1 - Boxplot do crescimento em diâmetro para a interação entre o tipo de propágulo (duas gemas, 
uma gema e apical) e a adubação de cobertura (A1, A2, A3 e A4) para o diâmetro do colo e a altura de 
mudas de Cordia trichotoma. 


Médias seguidas pela mesma letra não diferem estatisticamente pelo teste de Tukey (prob. = 0,05). Em 
que: 2G = 2 gemas; GU = uma gema e AP = gema apical. 


A interação entre os fatores para a altura (Figura 2), a maior média foi no tratamento 
com adubação de 200 g de sulfato de amônio e 150 g de cloreto de potássio, a cada 10 dias 
(A1) com duas gemas (4,03 cm). Mas esse não diferiu das adubações A2 (3,02 cm) e A3 
(3,80 cm) para duas gemas e para Aí (2,43 cm), A2 (2,54 cm) e A3 (3,12) para as mudas 
formadas com miniestacas apicais. 
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Figura 2 - Boxplot do crescimento em altura para a interação entre o tipo de propágulo (duas gemas, 
uma gema e apical) e a adubação de cobertura (A1, A2, A3 e A4) para o diâmetro do colo e a altura de 
mudas de Cordia trichotoma. 














Médias seguidas pela mesma letra não diferem estatisticamente pelo teste de Tukey (prob. = 0,05). Em 
que: 2G = 2 gemas; GU = uma gema e AP = gema apical. 


Mudas clonais de C. trichotoma formadas a partir de miniestacas de uma ou 
duas gemas ou apicais podem ser produzidas. O uso de 100 g de sulfato de amônio e 
300 g de cloreto para a adubação de cobertura a casa 30 dias não há perda do vigor 
de crescimento (Figuras 1 e 2). O crescimento inicial em diâmetro do colo e altura das 
mudas de C. trichotoma foi influenciado pela quantidade de sulfato de amônio e de cloreto 
de potássio utilizados (Figuras 1 e 2). As mudas são sensíveis às condições ambientais, 
porém, principalmente, para o diâmetro do colo, mas demonstram potencial de crescimento 


após 90 dias de estabelecimento em casa de sombra. 


41 CONCLUSÃO 
As mudas de C. trichotoma produzidas por miniestaquia de duas gemas, uma gema 
e apical, adubadas uma vez ao mês com 1000 g de sulfato de amônio e 300 g de cloreto de 


potássio apresentam bom crescimento em diâmetro do colo e altura. 
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RESUMO: A época do ano em que são coletadas 
as miniestacas para a produção de mudas 
exerce, em alguns casos, grande influência 
sobre o enraizamento adventício, e pode ser 
o principal fator de sucesso dessa atividade. 
Assim, objetivou-se verificar neste estudo o 
enraizamento adventício de miniestacas e 
selecionar clones de Cordia trichotoma (louro- 
pardo) coletadas em diferentes épocas do ano. 
O experimento é um bifatorial (época x clone) 
conduzido em câmara úmida, onde aos 30 dias 
de cultivo foram avaliadas a porcentagem de 
sobrevivência e enraizamento, o número de 
raízes, e o comprimento das três maiores raízes 
em miniestacas coletadas nas épocas de verão, 
primavera e outono. A época de coleta das 
brotações afetou o enraizamento adventício de 
miniestacas de louro-pardo, sendo a maior taxa 
de sobrevivência observada no outono, enquanto 
os melhores resultados para o enraizamento 
ocorreram durante as épocas de primavera e 
verão. Os clones 138M18, 13S8M10 e 13SM16 
foram menos afetados pelas diferentes épocas 
de coleta das miniestacas, porém com menor 
competência ao enraizamento adventício. Os 
clones 13SM17, 138M15, 13SM01, 13SM32, 
13SM09 e 135M26 foram selecionados por 
apresentar maior competência ao enraizamento 
adventicio na primavera, época que apresentou 
o maior índice de enraizamento de miniestacas 
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de louro-pardo. A época de coleta das miniestacas de louro-pardo afeta o processo de 
enraizamento adventício. Os clones avaliados apresentam diferença significativa entre si 
para o enraizamento adventício nas diferentes épocas avaliadas, o que demostra a variação 
existente entre os genótipos, possibilitando assim ganhos de seleção. 

PALAVRAS - CHAVE: Cordia trichotoma; melhoramento genético; propagação vegetativa. 


ROOTING OF MINI-CUTTINGS AND SELECTION OF CORDIA TRICHOTOMA 
CLONES IN DIFFERENT COLLECTION TIMES 

ABSTRACT: The collection time of mini-cuttings for plantlet production may affect the 
adventitious rooting competence, and it can be the main factor of plantlet production success. 
This work aiming to study the adventitious rooting of Cordia trichotoma mini-cuttings to 
maximize the genetic gain from selection for vegetative propagation in different collection 
times. The experiment was carried out in a humid chamber with mini-cuttings collected 
during the summer, spring and autumn seasons. The percentages of survival and rooting, 
number of roots, and length of three largest roots were evaluated at 30 days of cultivation. 
The collection time affected the adventitious rooting of C. trichotoma mini-cuttings. Autumn 
collections favored survival and spring and summer collections improved mini-cuttings 
rooting. The clones 13SM18, 135SM10 and 13SM16 were less affect by collection times, but 
they presented the lowest percentage of mini-cutting rooting. The clones 13SM17, 13SM15, 
13SM01, 135M32, 13SMO9 and 13SM26 were selected for rooting competence during spring 
collections, time that presented the highest rooting index of mini-cuttings of C. trichotoma. The 
collection time of the brown laurel minicuttings affects the adventitious rooting process. The 
clones evaluated showed a significant difference between themselves for the adventitious 
rooting in the different evaluated periods, which demonstrates the variation existing between 
the genotypes, thus enabling selection gains. 

KEYWORDS: Cordia trichotoma; genetic improvement; vegetative propagation. 


11 INTRODUÇÃO 


Com posição de destaque no desenvolvimento de uma economia de baixo carbono, 
o setor de base florestal diferencia-se por sua alta produtividade, tecnologia incorporada, 
e responsabilidade social. Em 2019 o setor cresceu 13,1%, contribuindo com 1,2% do 
Produto interno Bruto (PIB) nacional. As áreas destinadas às plantações florestais no 
Brasil totalizam 9 milhões de hectares. Dessa área, o eucalipto é o gênero mais plantado, 
representando cerca de 6,97 milhões hectares, enquanto, outras espécies tais como a 
seringueira e o paricá compõem apenas 0,39 milhões hectares (IBÁ, 2020). 

Afalta de informação sobre o plantio de espécies nativas para fins comerciais constitui 
um empecilho para a diversificação de matéria-prima no setor florestal, e indiretamente 
incentiva o cultivo de espécies exóticas que possuem cadeia produtiva definida, de alta 
confiabilidade econômica. No entanto, avanços tecnológicos aplicados ao melhoramento 
genético de espécies florestais nativas buscam suprir essa lacuna, proporcionando ganhos 
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genéticos significativos para diferentes características, tais como o enraizamento, qualidade 
de muda, crescimento dos indivíduos em campo, e densidade básica da madeira (PLATH 
etal., 2011). 

Através do melhoramento genético, é possível realizar a transformação do potencial 
genético das espécies nativas em recurso genético e, consequentemente, em produto com 
valor econômico no mercado. O estabelecimento adequado de um plantio florestal está 
diretamente relacionado à qualidade do material genético, tendo a produtividade maximizada 
quando bem adaptado às condições ecológicas e climáticas locais. Nesse sentido, o louro- 
pardo (Cordia trichotoma Vell.) é uma das espécies nativas mais promissoras para plantio 
no Sul, Centro-Oeste e Sudeste do Brasil, pois além de ocorrer naturalmente nessas 
regiões, apresenta rápido crescimento, tronco retilíneo, madeira de excelente qualidade e 
regeneração natural vigorosa (CARVALHO, 2003; ODA, et al., 2006). 

Para o estabelecimento de plantios florestais, a propagação vegetativa é um método 
especialmente atrativo, pois proporciona a operacionalização do processo de produção 
de mudas, e o estabelecimento de plantios florestais mais homogêneos. Contudo, a 
propagação vegetativa é afetada por diversos fatores que podem atuar isoladamente ou em 
conjunto na produção de mudas, tais como a espécie, o genótipo, as condições fisiológicas 
da planta matriz, a época e a posição de coleta, o grau de maturação das miniestacas, a 
presença de gemas e folhas, além de fatores ambientais, tais como temperatura e umidade 
relativa no ambiente de enraizamento (HARTMANN et al., 2011). 

A propagação vegetativa consiste no uso de partes de tecidos e células para a 
produção de indivíduos completos e idênticos à planta que lhe deu origem. A competência 
de regeneração das células ocorre através da totipotência celular, que se refere a 
capacidade na qual uma única célula carrega em seu núcleo o potencial de produzir todas 
as células diferenciadas do organismo de origem (HARTMANN et al., 2011). Dentre as 
técnicas de propagação vegetativa, a miniestaquia apresenta alta aplicabilidade e supera 
as dificuldades inerentes a produção de mudas seminais, por meio do aproveitamento do 
potencial juvenil dos propágulos para a indução do enraizamento adventício (ALFENAS et 
al., 2004). 

Alguns problemas que surgem em florestas clonais ocorrem pelo fato de o potencial 
de enraizamento das miniestacas não ser prognosticado antes da coleta dos brotos, devido 
a esse caractere ser altamente influenciado pelas mudanças sazonais e pelo manejo das 
matrizes, resultando em menor produtividade de miniestacas enraizadas por minicepa. A 
importância de se conhecer os fatores que afetam a formação de raízes e suas implicações 
está relacionada ao sucesso ou fracasso da produção de mudas (PENCHEL; LYRA, 1996). 

Em plantios clonais já estabelecidos, é recorrente o tombamento de indivíduos 
adultos após eventos climáticos severos. Isso ocorre, dentre outros fatores, pela má 
formação do sistema radicial das mudas em viveiro. Dessa forma, é de suma importância 


estudar as melhores condições de propagação de cada espécie, através do estabelecimento 
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de protocolos de propagação vegetativa. Além da obtenção de mudas de qualidade, 
os protocolos auxiliam no melhoramento genético para a miniestaquia, pois atenuam a 
influência ambiental e maximizam a expressão genotípica sob o fenótipo dos indivíduos 


(PENCHEL; LYRA, 1996). 
Nesse sentido, o desenvolvimento de um protocolo de propagação vegetativa 


para louro-pardo depende do entendimento e interação de fatores como as condições 
ambientais durante o processo de enraizamento. Como a estação do ano pode ser um fator 
decisivo para o sucesso da propagação vegetativa, o objetivo deste trabalho foi estudar 
o enraizamento adventício de miniestacas, bem como comportamento de diferentes 
genótipos de louro-pardo nas estações do ano, a fim de contribuir com o desenvolvimento 
de uma estratégia de seleção que maximize o ganho genético de seleção. 


21 METODOLOGIA 
O experimento foi conduzido em casa de vegetação do Núcleo de Melhoramento 


e Propagação Vegetativa de Plantas, do Departamento de Fitotecnia, da Universidade 
Federal de Santa Maria, campus Santa Maria. O minijardim clonal foi estabelecido em 
bandejas de polietileno (55 x 34 x 15 cm), onde foram cultivadas 12 minicepas em cada 
bandeja, no espaçamento 10 x 10 cm. O minijardim clonal foi conduzido no sistema fechado 
de cultivo sem solo, com areia grossa como substrato, adaptado de Bisognin et al. (2015). 

As minicepas de louro-pardo foram fertirrigadas com solução nutritiva de macro 
e micronutrientes duas vezes ao dia durante 15 minutos, com o auxílio de uma bomba 
submersa controlada por temporizador (KIELSE et. al., 2015). O pH da solução foi mantido 
entre 5,5 e 5,8 e a condutividade elétrica em 1 dS m”, ambos ajustados semanalmente. 

A coleta das miniestacas dos diferentes genótipos utilizadas para o estudo do efeito 
das estações do ano no enraizamento adventício de louro-pardo foi realizada nas estações 
verão, outono e primavera. Foram avaliados 11 clones (13SM17, 13S8M18, 13SMO07, 
13SM15, 13SM10, 13SM06, 13SM01, 13S8M32, 138M16, 13SMO9 e 135SM26). Apesar do 
louro-pardo ocorrer naturalmente em várias regiões climáticas do Brasil, a espécie cresce 
melhor em áreas sem geadas rigorosas, tornando sua recomendação restrita a locais 
com temperatura média anual maior que 18º C (CARVALHO, 1988). Em vista disso, a 
espécie diminui o crescimento no inverno não disponibilizando brotações para a coleta de 


miniestacas durante esse período do ano. 
As brotações coletadas de cada minicepa foram seccionadas em miniestacas de 


duas gemas contendo duas folhas com área original reduzida à 50%. As miniestacas 
tiveram sua porção basal imersa durante 10 segundos em solução hidroalcoólica de ácido- 
3-indolbutírico (AIB) de concentração 3.000 mg L'. Em seguida, as miniestacas foram 
cultivadas em tubetes de polipropileno de volume de 110 cmº com estrias longitudinais, 
contendo substrato comercial a base de casca de pinus e vermiculita na proporção de 2:1 
(VA). 
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O enraizamento das miniestacas foi realizado em câmara úmida com 85% 
de umidade relativa do ar, 27º C de temperatura média e nebulização automatizada 
realizada 12 vezes ao dia por 1 min. Foram avaliados a porcentagem de sobrevivência e 
enraizamento, o número de raízes e o comprimento das três maiores raízes aos 30 dias 
de cultivo. Foram consideradas enraizadas as miniestacas que apresentavam pelo menos 
uma raiz adventícia com comprimento igual ou superior a 0,1 cm. 

O experimento é um bifatorial (épocas x clones) conduzido em um delineamento 
inteiramente casualizado com 2 repetições contendo de 140 a 500 miniestacas cada, de 
acordo com a produtividade das minicepas. Para atender os pressupostos de normalidade, 
os dados de porcentagem foram transformados para arcosseno de (x/100)1/2 e de contagem 
e comprimento para (x+0,5)1/2. Os dados foram submetidos à análise de variância e as 
médias comparadas pelo teste de Tukey a 5% de probabilidade de erro. As análises foram 
realizadas com o auxílio do software R Core Team (2019). 


31 RESULTADOS E DISCUSSÕES 


As porcentagens de sobrevivência e enraizamento, o número de raízes e o 
comprimento das três maiores raízes apresentaram normalidade nos dados, além de 
significância para as épocas do ano estudadas. Tal comportamento demonstra que que ao 
menos uma época do ano difere estatisticamente para todas as variáveis analisadas. Além 
disso, foi verificado interação entre os fatores épocas do ano e clones de louro-pardo para 
os caracteres porcentagem de enraizamento e número de raízes. 

A época de coleta afetou todos os caracteres avaliados nas miniestacas (Tabela 1), 
sendo que o melhor resultado para a porcentagem de sobrevivência foi encontrado nas 
coletas de outono e primavera, diferindo significativamente daquelas de verão. Estes 
resultados podem estar associados as temperaturas mais amenas e menores taxas de 
irradiância no período de enraizamento, o que possibilitou reduzidas taxas de transpiração 
e, consequentemente, maiores taxas de sobrevivência das miniestacas (MARENCO et, 
al., 2001). Pimentel et. al. (2019) também encontraram maiores taxas de sobrevivência em 
miniestacas de llex paraguariensis A. St.-Hil coletadas no outono e inverno. 
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Estações do ano Sobrevivência (%) Enraizamento (%) NR C3MR 





Primavera 87,46 a* 46,76 a 1,68b 5,87 a 
Verão 80,40 b 29,19b 2,35 a 2,91b 
Outono 88,36 a 19,71 c 1,54b 3,16b 
Média 83,42 31,88 1,85 3,98 





Tabela 1 — Porcentagem de sobrevivência, de enraizamento, número de raízes (NR) e comprimento 
das três maiores raízes (C3MR) em miniestacas de Cordia trichotoma, após 30 dias de enraizamento 
em câmara úmida. 


*Médias seguidas por diferentes letra diferem estatisticamente pelo teste de Tukey (prob.=0,05). 


Já para a porcentagem de enraizamento, o melhor resultado foi encontrado com 
coletas na primavera, período em que a planta retoma o crescimento vegetativo, seguido 
do verão e do outono. Burin et. al. (2018) encontraram o mesmo comportamento para 
miniestacas de Cabralea canjerana (Vellozo) Martius, com porcentagem de enraizamento 
de 43,5%. Tal comportamento pode ser explicado pela capacidade de armazenamento de 
carboidratos durante a época de repouso vegetativo (inverno e outono), que se disponibiliza 
durante a primavera e o verão, época em que a planta inicia a indução do enraizamento 
adventício (BORTOLINI et al., 2008). 

O desenvolvimento das plantas está regulado por sinais ambientais, como, luz, 
temperatura e comprimento do dia (TAIZ; ZEIGER, 2013), que se alteram ao longo do ano, 
especialmente, em regiões com estações do ano bem definidas. Esses fatores ambientais 
alteram a fotossíntese, absorção de água e nutrientes, produção de fitormônios, entre 
outros, afetando o crescimento e desenvolvimento das plantas. 

Para o número de raízes, o melhor resultado foi observado para as coletas do 
verão, o que denota que esta estação favorece o processo rizogênico de miniestacas de 
louro-pardo. Além disso, a utilização de AIB, considerada a principal auxina sintética para 
o enraizamento de miniestacas, pode ter favorecido o processo de indução e formação 
radicial, promovendo o aumento da porcentagem, número e qualidade das raízes, bem 
como da uniformidade do enraizamento (FACHINELLO et al., 1995). 

Apesar das coletas da primavera terem apresentado baixo número de raízes em 
comparação ao verão, o comprimento de raízes foi significativamente superior ao observado 
quando as miniestacas foram coletadas no outono e verão (Tabela 1). Este resultado 
pode ser explicado pelo tempo da coleta das miniestacas ocorrer no início do período de 
crescimento vegetativo, momento em que se inicia a síntese de auxinas. As auxinas são 
fitormônios vegetais que estimulam o crescimento em elongação, divisão celular, formação 
de raízes adventícias entre outros efeitos fisiológicos (TAIZ; ZEIGER, 2013). 

A partir desses resultados observa-se que as miniestacas de louro-pardo coletadas 
no verão desenvolveram maior número e menor comprimento de raízes, o que pode ser 


explicado pela proximidade desta estação com a época de repouso vegetativo. Por outro 
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lado, quando a coleta é efetuada na primavera, momento em que se inicia o crescimento 
vegetativo, a indução do número de raízes é menor, mas o tempo de crescimento vegetativo 
perdura por mais tempo, possibilitando assim o alongamento das raízes. 

À medida que aumenta a disponibilidade de radiação solar no verão, e 
consequentemente, a temperatura do ambiente, elevam-se também as reações químicas 
de oxirredução no metabolismo da planta e o processo de divisão celular. Como a emissão 
de raízes ocorreu no período inicial do crescimento vegetativo, o processo de elongação 
correspondeu ao verão, período de maior crescimento vegetativo, o que favoreceu o 
alongamento das raízes (TAIZ; ZEIGER, 2013). Sá et al. (2018) também verificaram o efeito 
positivo no comprimento de raízes em miniestacas de Ilex paraguariensis coletadas no 
início do período reprodutivo (primavera). 

Com relação ao desempenho do enraizamento dos genótipos nas diferentes épocas 
de coleta, foram observadas diferenças significativas para a interação desses fatores (Tabela 
2). Os valores de enraizamento apresentaram variação de 47,6% na época de verão, de 
39,8% na época de outono e de 74,6% na primavera. Esse é um resultado relevante, 
que demostra a variação existente entre os genótipos para essa variável, possibilitando 
assim ganhos de seleção. De acordo com Burin et. al (2018) o número de miniestacas 
enraizadas seguido de porcentagem de enraizamento promove maior ganho de seleção e 
são características importantes para definir o potencial de um clone. 








Clones Verão Outono Primavera 
13SM17 28,00 bBCD* 41,01 DA 82,92 aA 
138M18 22,95 aCD 26,78 aAB 33,33 aBCD 
13SMO7 32,35 aABC 18,70 abBC 14,60 bCD 
138M15 32,67 bABC 19,37 cBC 55,38 aABC 
13SM10 16,67 aCD 17,42 aBCD 40,00 aBCD 
13SM06 14,28 LCD 14,17 bD 33,33 aBCD 
13SMO1 52,87 aA 17,46 bBCD 48,76 aABCD 
138M32 5,26 bD 15,53 bDBCD 60,00 aAB 
138M16 26,52 aCD 7,93 aCD 8,33 aD 
13SM09 51,42 aAB 21,66 bBC 64,28 aAB 
135M26 38,09 bABC 29,76 bAB 73,43 aAB 





Tabela 2 — Porcentagem de enraizamento em clones de Cordia trichotoma em três épocas de coleta, 
aos 30 dias de cultivo em câmara úmida. 


*Médias nas linhas seguidas por letras minúsculas e médias seguidas por letra maiuscula na coluna 
diferentes diferem a 5% de probabilidade pelo teste de Tukey. 
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A seleção para maior porcentagem de enraizamento indicou que o clone 13SM01, 
não diferindo estatisticamente dos clones 13SMOT7, 13SM15, 13SM0O9 e 135M26 apresentou 
o melhor resultado para as coletas do verão. Para o outono, o clone 13SM17 apresentou 
maior porcentagem de enraizamento, mas não diferiu de 135M18 e 135SM26. Nas coletas 
de primavera, o clone 135M17 apresentou a maior porcentagem de enraizamento, não 
diferindo dos clones 13SM15, 13SM01, 13S8M32, 13SM e 13SM26 (Tabela 2). Os clones 
13SM18, 13SM10 e 13SM16 não apresentaram diferença significativa para o enraizamento 
nas diferentes épocas de coleta de miniestacas. Apesar disso, esses clones apresentaram 
menores porcentagens de enraizamento comparado com os demais. 

Para Burin et al. (2018) a época de coleta de brotações afetou significativamente o 
enraizamento de miniestacas de Cabralea canjerana e o período de início de crescimento 
vegetativo resultou em maior porcentagem de enraizamento. Já para Pimentel et al. (2016) 
as maiores respostas ao enraizamento adventício de miniestacas de Ilex paraguariensis 
aos 60 dias de cultivo, foi obtido no outono e no inverno. Segundo Hartmann et al. 
(2011) algumas espécies apresentam indução de raízes quando coletadas em períodos 
específicos. Porém, para algumas, a indução de raízes adventícias ocorre em diferentes 
épocas do ano, mas são maximizadas em algum período específico. Além disso, Burin et 
al. (2018) ressalta que a seleção de clones superiores não deve ser afetada por mudanças 
de temperatura e fotoperíodo de crescimento vegetativo. 

Algumas espécies florestais, como o louro-pardo, passam por um período de 
repouso vegetativo, que consiste na redução da atividade morfofisiológica da planta 
durante as épocas mais frias (outono e inverno). Neste período, há a redução da área foliar 
e armazenamento de carboidratos que, reduzem a atividade metabólica até a chegada 
do inverno (FERRIANI et al., 2011). Isso pode acarretar a redução do enraizamento de 
miniestacas. Os melhores resultados para o enraizmento do louro-pardo em diferentes 
clones foram obtidos na primavera, seguido do verão. Isso denota que estas duas estações 
favorecem o processo rizogênico de miniestacas de louro-pardo. Hartmann et al. (2011) 
salientam que estacas coletadas na primavera e no verão tendem a ter maior facilidade 
de indução do enraizamento adventício em função de maior crescimento vegetativo nessa 
época. 

Quanto ao número de raízes, houve diferença significativa para as três estações do 
ano. O maior número de raízes foi observado no verão e na primavera. Em contrapartida, os 
menores resultados foram obtidos no outono, devido à proximidade de repouso vegetativo. 
Para o verão e outono não houve diferença significativa para o número de raízes dos clones 
estudos. Em contrapartida no outono o clone 135M18 apresentou menor número de raízes 
diferindo dos demais clones (Tabela 3). 
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Clones Verão Outono Primavera 





13SM17 3,28 DA* 2,44 bA 4,44 aA 
135M18 1.85 aA 1,88 aA 1,2aB 

13SM07 2,77 ah 1,90 aA 2,15 aAB 
135SM15 2,51 abA 2,10 bA 3,09 aAB 
135SM1O 1,422 aA 1,92 aA 1,75 aAB 
135M06 2,3 aA 4,00 aA 1,85 aAB 
13SM01 3,23 aA 2,09 bA 3,05 aAB 
135M32 2.50 aA 1,50 aA 2,00 aAB 
135M16 1,38 aA 3,00 aÃ 2,00 aAB 
13SM09 2,06 aA 2,15 aA 1,78 aAB 
135M26 2,60 aA 2.54 aA 2,55 aAB 





Tabela 3 - Número de raízes em clones de Cordia trichotoma em três épocas de coleta, aos 30 dias de 
cultivo em câmara úmida. 


*Médias nas linhas seguidas por letras minúsculas e médias seguidas por letra maiuscula na coluna 
diferentes diferem a 5% de probabilidade pelo teste de Tukey. 


Para Pimentel (2016) o menor número de raízes foi obtido na primavera, sem diferir 
estatisticamente do verão. O outono e o inverno foram as épocas que apresentaram maior 
formação de raízes nos propágulos vegetativos. Em contrapartida, para Somavilla (2018) o 
maior número de raízes aos 30 e 60 dias de cultivo, foi obtido no verão para o louro-pardo. 
Este comportamento está de acordo com os resultados obtidos para esta espécie neste 
estudo, que indica a primavera e o verão como estações favoráveis para a formação de 
promórdios radiculares, devido a maior atividade fisiológica e crescimento vegetativo nesta 
época quando comparado a outras épocas do ano (HARTMANN et a., 2011). 

De maneira geral, os resultados deste estudo mostram que as miniestacas de 
louro-pardo apresentaram maior sobrevivência quando as miniestacas foram coletadas no 
outono, enquanto, o enraizamento adventício ocorreu de forma mais eficiente durante a 
primavera e o verão. Assim, para realizar o enraizamento adequado de miniestacas de 
louro-pardo, visando ganhos com a seleção dos melhores genótipos, e diminuir o efeito 
do ambiente, a coleta de material propagativo deve ser realizada durante a primavera e o 


verão. 


41 CONCLUSÃO 


Aépoca de coleta das brotações das minicepas influencia o processo de enraizamento 
adventício das miniestacas de louro-pardo, sendo a maior taxa de sobrevivência observada 
no outono, enquanto os melhores resultados para o enraizamento ocorreram durante as 
épocas de primavera e verão. Os clones avaliados apresentam diferença significativa entre 
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si para o enraizamento adventício nas diferentes épocas avaliadas, o que demostra a 


variação existente entre os genótipos, possibilitando assim ganhos de seleção. 
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ABSTRACT: Drought impacis in Northeast Brazil, 
which tend to intensify due to climate change 
and overexploitation of the caatinga biome, 
have repeatedly brought famine, mass migration 
and social conflicts in this region. Drought 
prediction and monitoring, however, remain 
a central research theme. In water resources 
management in semiarid regions such as the 
Northeast of Brazil, itis fundamental to have tools 
to aid decision making. This work presents three 
components of the so-called SIGES (Drought 
Management System), the items related to 
drought prediction and monitoring, as well as 
many reservoir operation methodologies for water 
scarcity situations. Statistical models and artificial 
neural networks were used for drought prediction. 
Aiming at monitoring precipitation, several indices 
were adapted and incorporated into a drought 
basic characteristic monitoring system (duration, 
severity and intensity), so that different mitigating 
actions could be implemented in accordance with 
the values reached by these parameters. The 
following indices have been used for this purpose: 
Rainfall Anomaly Index (RAI); Bhalme & Mooley 
Drought Index (BMDI); Lamb Rainfall Departure 
Index (LRDI). Finally, new adaptive reservoir 
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NORTHEAST BRAZIL 


operation strategies for water scarcity situations 
are presented and discussed. The described 
components were applied to several States of 
Northeast of Brazil, specially, at Piancó-Piranhas- 
Açu River Basin (Paraiba and Rio Grande do 
Norte States) and Guaribas River Basin (Piauí 
State). 

KEYWORDS: water resources, drought index, 
drought, decision making, Northeast Brazil 


FERRAMENTAS PARA A TOMADA 
DE DECISÕES ESTRATÉGICAS NA 
GESTÃO DE RECURSOS HÍDRICOS EM 
SITUAÇÕES DE SECA E VARIABILIDADE 
CLIMÁTICA NO BIOMA CAATINGA, 
NORDESTE DO BRASIL 


RESUMO: Os impactos da seca no Nordeste 
do Brasil, que tendem a se intensificar devido 
às mudanças climáticas e à superexploração do 
bioma caatinga, têm trazido fome, migração em 
massa e conflitos sociais na região. A previsão e 
o monitoramento da seca, no entanto, continuam 
sendo um tema central de pesquisa. Na gestão 
de recursos hídricos em regiões semiáridas como 
o Nordeste do Brasil, é fundamental dispor-se de 
ferramentas que auxiliem na tomada de decisões. 
Este trabalho apresenta três componentes do 
denominado SIGES (Sistema de Gerenciamento 
de Secas), os itens relacionados à previsão e ao 
monitoramento de secas, bem como diversas 
metodologias de operação de reservatórios 
em situações de escassez de água. Modelos 
estatísticos e redes neurais artificiais foram 
usados para a previsão de secas. Com o objetivo 
de monitorar a precipitação, diversos índices 
foram adaptados e incorporados a um sistema de 
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monitoramento das características básicas da seca (duração, severidade e intensidade), para 
que diferentes ações mitigadoras pudessem ser implementadas de acordo com os valores 
alcançados por esses parâmetros. Os seguintes índices foram usados para este propósito: 
Rainfall Anomaly Index (RAI); Bhalme & Mooley Drought Index (BMDI); Lamb Rainfall 
Departure Index (LRDI). Por fim, novas estratégias adaptativas de operação de reservatórios 
para situações de escassez de água são apresentadas e discutidas. Os componentes 
descritos foram aplicados em diversos Estados do Nordeste do Brasil, em especial, na Bacia 
do Rio Piancó-Piranhas-Açu (Estados da Paraíba e Rio Grande do Norte) e na Bacia do Rio 
Guaribas (Estado do Piauí). 

PALAVRAS-CHAVE: recursos hídricos, índices de secas, secas, tomada de decisão, 
Nordeste do Brasil. 


11 INTRODUCTION 


Drought impacts in Northeast Brazil, which tend to intensify due to climate change 
and overexploitation of the caatinga biome, have repeatedly brought famine, mass 
migration and social conflicts. Drought prediction and monitoring, however, remain a central 
research theme for this region. Marengo (2006), Brito et al. (2017) conducted studies on the 
impacts of global climate change for many areas of Brazil, as Brazilian Amazon, Northeast, 
Pantanal and the Prata River Basin, showing precipitation and temperature anomalies, and 
water balance for the XXI century. The semiarid northeastern presenting short but crucially 
important rainy season in the current climate could, in a warmer climate in the future, 
become arid. 

Tothe northeast, Salati et al. (2007) assessing climate variability in the region, showed 
that the average temperature increased by 0.6ºC within the period 1991 to 2004, when 
referred to the period from 1961 to 1990. For maximum temperature values they indicated 
a 0.6ºC increase and for minimum temperature values a 0.5ºC raise. The precipitation 
decreased 153 mm, a 11.6% drop. Nobre et al. (2004) indicate that the future biomes 
distribution in South America may be affected by the combined impacis of climate and land 
use changes, which can take the system to the savannization of parts of the Amazon and 
the desertification of Northeast Brazil. 

Several studies have indicated the influence of numerous atmospheric phenomena 
on rainfall in Northeast Brazil (Moura and Shukla, 1981; Hastenrath, 1984; Freitas, 1996; 
Freitas and Billib, 1997; Uvo et al., 1998; Moscati and Gan, 2007). Also several climatological 
studies have indicated the existence of a strong relationship between sea surface 
temperature distribution (SST - sea surface temperature) along the tropical Atlantic basin 
temperature and the semiarid northeastern Brazil precipitation, as well as a decadal trend 
associated with changes in the meridional position of the ITCZ - Intertropical Convergence 
Zone (Moura and Schukla, 1981; Billib and Freitas, 1996). These phenomena are indicative 
to be related to climate variability and extreme droughts and floods in the region. 

Freitas (1996) presented a methodology for integrated regional drought analysis, 
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which briefly consists of the following topics: (1) definition of different types of droughts, (2) 
forecasting and monitoring: (3) water resources management and optimization, (4) effects 
evaluation and (5) planning of mitigating actions. Freitas (1996), and Freitas & Billib (1997) 
demonstrated the feasibility of using prediction models for the dry northeastern Brazil: 
statistical-probabilistic models and models based on neural networks. 


21 THE SIGES (DROUGHT MANAGEMENT SYSTEM) 


In water resources management in semiarid regions such as the Northeast of Brazil, 
itis fundamental to have tools to aid decision making. This work presents three components 
of the so-called SIGES (Drought Management System), the items related to drought 
prediction and monitoring, as well as many reservoir operation methodologies for water 
scarcity situations. Statistical models and artificial neural networks were used for drought 
prediction. 





Figure 1: Drought Management System - SIGES (Freitas, M. A. S. et al., 2019). 


Among the various SIGES applications, will be reported those related to shortage 
impact mitigation based on forecasting and monitoring information for these periods. In 
semiarid regions with high evaporation rates the water resources manager often faces the 
problem of determining minimum water losses conditions for surface reservoir operation, 
considering different uses, their priorities shortages or rationing situations. 


2.1 Strategic Management 


Aiming at monitoring precipitation, several indices were adapted and incorporated into 
a drought basic characteristic monitoring system (duration, severity and intensity), so that 
different mitigating actions could be implemented in accordance with the values reached by 
these parameters. The following indices have been used for this purpose: Rainfall Anomaly 
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Index (RAI); Bhalme & Mooley Drought Index (BMDI); Lamb Rainfall Departure Index 
(LRDI). Finally, new adaptive reservoir operation strategies for water scarcity situations are 
presented and discussed. The described components were applied to all federal States of 
Northeast of Brazil. 

Due to the difficulty of having a single drought phenomenon definition, drought index 
determination is also problematic. Drought severity is usually expressed as a function of the 
mean value of one or several climate parameters. The most usual drought indexes express 
the total precipitation in given region as a known average data percentage. Thus, one can 
compare the current value with the average value of weekly, monthly, semesterly or annually 
series. The first index, however, that was incorporated into the Drought Management 
System - SIGES has made it possible to compare rainfall deviations related to the normal 
conditions of various regions is the Rainfall Anomaly Index (RAI), described by Rooy (1965) 
and pioneered in Brazil by Freitas (2005). 
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Figure 2: Bhalme & Mooley Drought Index (BMDI) 


2.2 Tatic Management: Reservoir Operation Optimization using Climate 
Information, Negotiated Water Allocation and Resilience Strategies 

Among the various SIGES applications, we will report now those related to shortage 
impact mitigation based on forecasting and monitoring information for these periods. In 
semiarid regions with high evaporation rates the water resources manager often faces the 
problem of determining minimum water losses conditions for surface reservoir operation, 
considering different uses, their priorities shortages or rationing situations. One developed 
method for dealing with this problem is performed using the so-called Risk Aversion Curve 
(Freitas and Gondim Filho, 2007). 

In many hydrographic basins, due to the lack of fluviometric data, it is often necessary 
to obtain a series of affluent flows into the reservoir using rainfall-runoff or synthetic 
streamflow generation models. For the generation of synthetic flows in intermittent rivers 
typical of semiarid regions, Freitas (1995, 2010) presented the SAGE system - Stochastische 
AbflussGEnerierungsmodelle, consisting of the several models, such as: i) Two-tier model 
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(PAR (1) / GAR (1) of Fernandez & Salas (Freitas, M. A. S.; Freitas, G. B., 2019a); ii) 
Disaggregation model / AR (1) by Valencia & Schaake (Freitas, M. A. S.; Freitas, G. B., 
2019b), etc. Later, Freitas (2010) presented applications of the model ARRF - Alternating 
Renewal Reward — Fragment (Freitas, M. A. S.; Freitas, G. B., 2019c) to various basins in 
northeastern Brazil. 

Water deficits as result of non-sustainable water resources development, are the 
main reason for water related conflicts in water scarcity affected regions, especially due to 
strong competition between the different water users. According to Rusteberg and Freitas 
(2018), new approaches, methods and tools have been developed to support the integrated 
planning and management of water resources towards the sustainable development of the 
region. Non-conventional water resources, such as treated wastewater, brackish or imported 
water needs to be part of an integrated strategy for IWRM implementation. 

According to Lopes & Freitas (2007), water allocation in Brazil, historically, is 
characterized by a strong intervention of the public sector. Water allocation can be 
understood as an Integrated Water Resources Management (IWRM) measure that aims to 
provide water to current and future users of the water resources system, matching water 
supply and demand, even meeting environmental demands, being in line with strategic 
management objectives. In this sense, there are several mechanisms of water allocation, 
which operate from public authority guidelines, negotiation processes among water users or 
from technical concepts, such as the limits of the use of water bodies, or economic, such as 
charging for water use (Rusteberg and Freitas, 2018). 


31 PIANCÓ-PIRANHAS-AÇU RIVER BASIN 


For the so-called Negotiated Water Allocation a methodology is described in Freitas, 
M.A.S. (2003); ANA/GEF/PNUMA/OEA (2004); and improved by Freitas, M.A.S (2010), 
Martins et al., (2013). This methodology was largely internalized in the administrative 
procedures of the National Water Agency, through Technical Note No. 10/2015/COMAR/ 
SRE (2015) and has been applied throughout the semi-arid region of North-East Brazil. 

Figures 3 shows the Risk Aversion Curves - CAR proposed for the Curema-Mãe 
D'água system reservoirs (Freitas and Gondim Filho, 2007). The operation of reservoirs 
based on the Risk Aversion Curves is guided by a probability distribution of flows flowing 
into the reservoir system and is useful for determining the water allocation among users 
with different guarantees. The figures show the useful reservoir volume (%) on June 1st 
(recommended period for the Negotiated Water Allocation), for the different years (including, 
in red, the years of the 2012-2020 drought period), where it is possible to obtain the flow to 
be released. ANA, within its sphere of competence, has placed itself in articulation with other 
bodies of the Federal Government, state governments and water resources management 


institutions, supply companies, city halls and basin committees in search of solutions that 
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can minimize the effects of this severe drought. 


recommended month for 
negotiated water allocation 


Risk Aversion Curves for the 
Curema and Mãe D'água Reservoir 
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Figure 3: Risk Aversion Curves for the Curema and Mãe D'água Reservoir 


41 GUARIBAS RIVER BASIN 


The Bocaina dam, in the semi-arid region of the State of Piauí, with a capacity of 106 
hm? of water, was built with the objective of regularizing the Guaribas River, among other 
uses. Over the past three decades, use for irrigation and water supply in neighboring cities 
has increased significantly. The objective of this study is to analyze, by means of several 
multicriteria methods, the alternatives of operation of the afore mentioned reservoir to meet 
multiple uses. Five alternatives for reservoir operation, socioeconomic and environmental 
development were analyzed. For this, four different multicriteria techniques were used: i) 
weighted average method — WAM (Zadeh, 1963); ii) compromise programming — CP (Yu, 
1973); iii) Promethee method with weighted averages (Promethee WAM) and iv) Promethee 
method (Brans and Vincke, 1985). 

The decision-making process in water resources management currently involves 
multiple objectives and multiple decision-makers (river basin committees, consortia, etc.). In 
general, conflicts of interest are established between groups with different views about the 
goals to be adopted in the planning and management of water resources. For groups linked 
to pure and simple economic development, the net economic benefit should be maximized, 
since economic values express the interest of society. On the other hand, environmental 
groups preach the preservation of the environment in its natural form and are opposed to 
any intervention that may transform it. Between these two extreme possibilities. 
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The Bocaina Reservoir (Figure 4), with a capacity of 106 hm? of water, was built 
under the responsibility of DNOCS - National Department of Works Against Drought - and 
executed by the Engineering and Construction Battalion, with the objective, among other 
aspects, of regularize the Guaribas River. The estimated area by the viability studies of 
the Bocaina Reservoir reached 2,000 (two thousand) irrigable hectares, according to the 
Bocaina Dam - Hidroterra SA Executive Project - located in the municipalities of Bocaina, 
Sussuapara and Picos, with excellent physical and chemical conditions for intensive 
irrigation use, at a relatively low cost, considering that there was no cost for deforestation, 
systematization and settlement of producers. 

The study of alternatives (A1 to A5) was carried out in two stages: i) Use of rainfall- 
runoff hydrological models (Freitas, 2010), stochastic streamflow generation models 
(Freitas, 2010) and Monte Carlo analysis, as well as models for simulation / optimization 
of the reservoir operation; ii) Use of multicriteria techniques, such as weighted averages 
method, compromise programming (CP) method and Promethee method. 


4.1 Simulation of the Bocaina Reservoir Operation 


The current complexity of problems related to the management of water resources 
requires the use of instruments and techniques capable of assisting in decision making, 
especially in periods of scarcity. Among these techniques, two deserve mention: mathematical 
simulation and flow network modeling (Azevedo et al., 1997). Mathematical simulation 
is a very flexible technique and is therefore widely used, however, as a disadvantage of 
not offering users the chance to restrict decision-making space and, therefore, problem 
solving is achieved through a “trial and error” process. On the other hand, the network flow 
modeling, in addition to reducing the decision space to a set of viable solutions, can be 
optimized through linear programming. 

The water availability in a river is generally used to water supply to the population 
and industries; animal supply; irrigation of agricultural areas; ecosystem conservation; 
wastewater dilution; hydro-energy production and navigation. In the specific case of the 
Guaribas River, downstream ofthe Bocaina Dam, only the last two uses are not contemplated. 

The simulation method uses the continuity equation. In the reservoir operation 
simulation, the process is initiated considering the average accumulated volume of the dam. 
Then, several simulations are made to define regularized volumes, associated with various 
guarantees or risks. In relation to irrigation, simulations were made for several demand 
areas, based on the estimated area potential for the downstream valley, forecasting monthly 
demands based on percentages of annual demand, proportional to the monthly deficits 
determined based on the water balance in the area. 
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Figure 4: Bocaina Dam location, in the semi-arid region of the State of Piauí. 


The failure risk was calculated based on the number of failures in the period of 
interest divided by the number of time intervals for that period, that is, the simulated period. 
The simulation model of the reservoir operation used was the one shown in Figure 5, 
developed by Prof. Dr. Fontane, at the University of Colorado. The model uses Excel's 
Solver to maximize the annual demand met for the various uses in the basin (irrigation, 
human supply and animal feed) and to determine the useful (active) volume needed to meet 
a given demand (reservoir project). The minimum and maximum volume of the Bocaina 
Dam is 10 and 106 hmº, respectively. The maximum annual demand met for the simulated 
series was 22.27 hm. 


4.2 Multicriteria Analysis in the Alternatives Study 


For the second stage, the MCDA spreadsheet - Multi-Criteria Decision Analyzes 
(Figure 8), also developed by Prof. Dr. Fontane, from the University of Colorado. Five 
alternatives (Al to A5) were considered for dam operation / economic and environmental 
development in the region. The alternatives were evaluated according to five criteria: i) 
irrigation; ii) human and animal supply; iii) the environment; iv) safety and reliability of the 
dam's operation and v) economic development. Up to five sub-criteria were established 
for each of the above criteria. Figure 5 shows the values of the criteria and sub-criteria for 
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each alternative. Five alternatives of relative importance between the criteria were analyzed 
(G1 to G5). The relative weights between the criteria and sub-criteria were raised based on 
surveys with specialists in the area and interviews between users. 
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Figure 5: Decision making multicriteria models. 

















Four different techniques were used: i) weighted averages method (WAM); ii) 
compromise programming (CP); iii) Promethee method with weighted averages (Promethee . 
WAM) and iv) Promethee method. As a result, there are the “scores” and the ranking 
(hierarchy) among the various alternatives. In the Bocaina Basin example, the following 
priority was found: A1 5 AZ 5 A55A354A4. 


5 | CONCLUSIONS 


With the increase in demand for water resources, there is a tendency to intensify 
conflicts with other users, especially during periods of long-term and severe droughts, such 
as that of 2012-2020. In situations like this, the integrated management of water resources 
stands out as an alternative to the resolution of conflicts, contributing to minimize the 
impacts of long periods of drought. In this article, the application of some of the techniques 
and methods developed, especially for semi-arid regions, in the Piancó-Piranhas-Açu 
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hydrographic basin as well as Bocaina Dam / Guaribas hydrographic basin were presented. 
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RESUMO: A presente pesquisa foi desenvolvida 
com o objetivo de determinar o consumo do 
óleo diesel na produção de cereais e frutas; 
avaliar a consequente quantidade de dióxido 
de carbono (CO,) emitida para a atmosfera; e, 
ainda, relacionar projetos florestais que, por meio 
da fotossíntese, durante os ciclos de crescimento 
das árvores, contribuem para o sequestro de 
parte dessas emissões. Foram selecionadas 41 
culturas: 38 perenes e 03 temporárias, tendo 
sido resgatados dados de produção de 13 anos 
consecutivos, de 2000 a 2012. Assim, com o 
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modelo matemático resultante, foi possível obter 
as emissões de COZ de tais cultivos - individuais 
ou agrupados -; bem como, comparar os 
resultados das plantações perenes e temporárias. 
O estudo de caso selecionado apresenta, como 
exemplo, uma comparação entre o consumo 
de óleo diesel e as respectivas emissões de 
CO2 das lavouras de BANANA (Musa sp.) e 
SOJA (Glycine max), com base na produção 
brasileira. Em relação ao sequestro de carbono 
(C) foram elaboradas quatro alternativas de 
mitigação florestal: i) PROJETO DE PROTEÇÃO 
FLORESTAL; ii) REFLORESTAMENTO com 
espécies de rápido crescimento; iii) SISTEMAS 
DE AGROFLORESTAIS (1º ciclo); e (iv) 
PROJETO DE ARBORIZAÇÃO URBANA. Para 
cada opção, há a indicação do número de 
árvores a serem plantadas ou hectares a serem 
protegidos / enriquecidos para efeito de redução, 
ainda que parcial, da quantidade de CO2 emitida 
para a atmosfera. 

PALAVRAS - CHAVE: Agricultura & floresta; 
consumo de diesel; Emissão de COZ; sequestro 
de Carbono com projetos florestais. 


DIESEL CONSUMPTION ASSESSMENT 
MODEL IN AGRICULTURE, WITH 
CO2 EMISSION ESTIMATES FOR 

ATMOSPHERE AND FOREST PROJECTS 
CARBON HIJACKING. CASE STUDY: 
BANANA X SOY 


ABSTRACT: The present research was 
developed with the objective of determining the 
consumption of the diesel oil in the production 
of cereals and fruits; assess the consequent 
amount of CO, emitted into the atmosphere; and, 
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in addition, list forest projects that, through photosynthesis, during tree growth cycles, can 
contribute to the total or partial sequestration of these emissions. 41 crops were selected: 38 
perennials and 03 temporary, for which were rescued production data ffom 13 consecutive 
years, ffom 2000 to 2012. Thus, with the resulting mathematical model, it was possible to 
obtain CO2 emissions from all crops - individual or grouped -, as well as to compare the 
results of perennial and temporary plantations. The selectees study case show, for example, 
a comparation the diesel consumption and the respective CO2 emissions of the crops: 
BANANA (Musa sp.) and SOYBEAN (Glycine max), based on the Brazilian production. About 
regarding the sequestration of carbon (C) four forest mitigation alternatives were elaborated: 
i) FOREST PROTECTION PROJECT; ii) REFORESTATION with fast growing species; 
iii) AGROFLORESTAIS SYSTEMS (1st cycle); and (iv) URBAN ARBORIZATION PROJECT. 
For each option, there is the indication of the numbers of trees to be planted or hectares 
to be protected / enriched for a purpose of the reduction of the quantity of CO2 issued for 
atmosphere. 

KEYWORDS: Agriculture and forest; diesel consumption; CO2 emission; carbon sequestration 
with forests projects. 


11 INTRODUÇÃO 


A modernização do setor agrícola no Brasil e em outros países da América Latina 
começou na década de 1960, com investimentos governamentais em: infraestrutura, 
implantação de sistemas de comunicação; ampliação da assistência técnica, garantia 
de preços mínimos e liberação de linhas de crédito. Mas esta modernização não inclui a 
totalidade dos cultivos perenes e dos produtos vegetais do extrativismo. 

Além disso, na época, não havia preocupação com as emissões de gases de efeito 
estufa, mas, sim, uma visão desenvolvimentista predominante em todos os setores da 
economia. Apenas na década de 1980, as mudanças climáticas despertaram preocupação 
na comunidade científica internacional, devido às evidências que correlacionam o aumento 
das emissões de gases de efeito estufa (GEE) com o aquecimento global. Somente em 1990, 
com o primeiro relatório do Painel Intergovernamental sobre Mudanças Climáticas (IPCC) 
[1], foi alcançado um relativo consenso sobre este tema, mas não excluiu divergências 
metodológicas entre países desenvolvidos e em desenvolvimento, cujos entendimentos 
sobre as finalidades e capacidades de divulgação dos parâmetros agroflorestais são 
opostos. Há, obviamente, as diferenças entre as atividades praticadas dos pequenos e 
médios agricultores e as ações de marketing conduzidas pelos grandes investidores do 
agronegócio [2] [3]. Assim, a escolha da BANANA e da SOJA baseou-se no fato desses 
cultivos ocuparem os extremos no ranking do consumo de diesel e das respectivas emissões 
de dióxido de carbono (CO,). Além, dos tratos culturais serem de dois grupos distintos de 
cultivos: respectivamente, permenentes e temporários [4]. 

Provavelmente, por esse motivo e por outras razões estruturais, verifica-se, há 
muito tempo, com poucas exceções, a demora no aprimoramento da catalogação dos 
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dados da produção agrícola e florestal. Observa-se, por exemplo, que, no Brasil e em 
nuitos países do hemisfério sul, dificuldades para o levantamento dos dados primários, 
exigindo inferências diversas, que incluem, dentre outros, o cálculo das taxas de fixação de 
carbono de espécies arbóreas. Logo, poucos pesquisadores se aventuram a enfrentar as 
dificuldades inerentes deste tipo de estudo devido ao elevado grau de incerteza, causado, 
também, pelas consideráveis diferenças de solo, clima e vegetação das regiões terrestres. 

Dessa forma, apesar dessa constatação, a presente pesquisa agroflorestal teve 
como finalidade superar as seguintes dificuldades: 


* A escolha das culturas (perenes e temporárias) e das ferramentas de proces- 
samento digital; 


*-— Adecisão de agrupar em um mesmo componente do modelo de cálculo culturas 
agrícolas com processos produtivos relativamente similares; 


* A sistematização de diversos dados agroflorestais para subsidiar a pesquisa 
técnico-científica, com a elaboração de um modelo de cálculo de apoio suple- 
mentar baseado na Ciência Florestal; 


*-— Aquantificação estimada do volume de óleo diesel usado na mecanização e na 
produção e transporte das diferentes safras; e 


* | Conversão matemática de medidas (área e volumes) em parâmetros de ener- 
gia. 

Este estudo também mostra, na medida do possível, as grandes diferenças entre 
os modelos de produção das safras do agronegócio (SOJA, cana-de-açúcar e algodão) 
[5] [6] e os gastos menores feitos para a produção de BANANA, de outras frutas tropicais 
[7] [8] [9] e de produtos do extrativismo vegetal [10]. Consequentemente, em função das 
peculiaridades e características das lavouras e florestas, o modelo em apreço buscou 
gerar indicadores estatísticos de projetos de reflorestamento, de arborização e de ações 
de proteção agroecológica, com o intuito de estimular plantios de árvores, tanto por 
produtores rurais, como por instituições públicas, voltados para redução das emissões de 
dióxido de carbono (CO,). Devido à inconsistência dos dados técnicos e/ou de informações 
socioeconômicas, houve a necessidade de algumas inferências, enfatizando: 


* Nível de mecanização de um processo produtivo [11] [12] [13]; 
* | Aimportância atribuída a um tipo de tratamento cultural em relação a outros; e 


* Padrão de vida dos produtores rurais (pequenos, médios ou grandes produto- 
res [14]). 

Ao optar neste estudo pela comparação entre uma cultura tradicional e outra 
altamente mecanizada, a intenção é demostrar, a partir das experiências agrícolas 
brasileiras e latino-americanas, a grande diferença que existe entre as atividades dos 
pequenos e médios agricultores e dos grandes. investidores do agronegócio. 


Conceitos e Conhecimentos de Métodos e Técnicas de Pesquisa Científica em Engenharia 


Florestal 2 Capítulo 5 DE 


21 JUSTIFICATIVAS E OBJETIVOS 


Modelo matemático que visa sensibilizar a sociedade em geral e os empresários 
do setor agrícola, em particular, para a necessidade de investimentos em projetos de 
recuperação de paisagens florestais, que contribuem, adicionalmente, para a proteção dos 
solos, da biodiversidade; bem como, para a mitigação das emissões de CO,, derivadas das 
atividades agrícolas. 


2.1 Objetivos Gerais 


* | Propor modelo de sistematização de dados de produção agrícola para cálculo 
do consumo de óleo diesel e das respectivas emissões de CO,; 


*— Indicar opções florestais e agroflorestais para apoiar projetos de sequestro de 
carbono, apontando o número de árvores e hectares a serem plantados ou 
protegidos, usando como exemplos comparativos as safras de BANANA e de 
SOJA; e 


* Contribuir para a melhoria da reputação das empresas que, de fato, optarem 
pela neutralização das emissões de CO, , contribuindo, assim, para a gerar 
oportunidades de trabalho no meio rural e o aumento da renda dos pequenos 
produtores. 


31 METODOLOGIA 


3.1 Planilhas de cálculo - Dados de entrada 


Foram selecionados 41 produtos diferentes (agrícolas e extrativistas), sendo: 38 
perenes e 03 temporários (Tabelas 1), das quais resgatou-se dados de produção de 13 
anos consecutivos, de 2000 a 2012. 


BANANA 
Café e cacau 
Coco e palma 
Frutas Cítricas: Laranja, Limão e Tangerina 
Frutas Tropicais: Abacate, Anonáceas (atemóia, graviola, abacaxi, etc.), Caju, Goiaba, 
Cultivos Mamão, Manga 
Permanentes | Frutos subtropicais: caqui, figo, maçã, pêra + marmelo e pêssego + ameixa + 
nectarina 
Uva, Maracujá e Pimenta do Reino 
Extrativismo vegetal (palmeiras): Açaí, Babaçu, Buriti + Carnaúba, Palmito e Piaçava. 
Extrativismo vegetal (Madeiras nobres): Borracha, Castanha, Erva Mate, Guaraná, 
Mangaba, Pequi e Umbu. 


Algodão (herbáceo e arbóreo) 
Cana de açúcar 
SOJA 


Cultivos 
Temporários 





Tabela 1: Culturas pesquisadas e principais agrupamentos 
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A Tabela 2 mostra os dados da produção, produtividade e transporte interno das 


lavouras de soja e banana, fruta que é produzida em 130 países, que compõe a dieta 


alimentar de pessoas de todos os continentes [7]. 


BANANA 2000 
Dados de Entrada 


Área Plantada (1000 ha) [4] [15] [32] 


2001 


2011 


2012 





Produtividade (toneladas /ha) [15] 





Dias de transporte / ano -70%-função das zonas: produção e colheita [32] 





Distância: da área de plantio ao comércio / indústria (km) [32] 
Capacidade média do caminhão: “Toco (- toneladas) [16] 





Consumo estimado de óleo diesel por caminhão (kmh) 
Conversão: 1.0 tonelada de diesel 0,848 tEP 
Densidade do Óleo Diesel 0,840 t/m? 





Preço médio do Óleo Diesel na bomba 








13,988.4 





2011 


24,052.4 





2012 


27,121.5 





Área Plantada (Mil hectares) [4] [15] [32] 
Área de Colheita (Mil hectares ) [15] 
Produtividade (toneladas / ha) [15] 


13,854.1 
2,40 


24,050.2 
14,36 


27,716.6 
12,96 








Dias de transporte / ano -70%-função das zonas: produção e colheita [32] 
Distância: da área de plantio ao comércio / indústria (km) [32] 


136 
85 


136 


136 
85 





Capacidade média de carga de caminhões com 02 e 03 eixos (- t) [16] 
Consumo estimado de óleo diesel por caminhão (TRUCK - km/1) [16] 
Conversão: 1.0 tonelada de diesel = 0,848 tEP 


30,0 
3,5 
0,848 


30,0 
3,5 
0,848 





Densidade do Óleo Diesel = 0,840 t/mº 
Preço médio do Óleo Diesel na bomba 


TABELA 2: Dados Básicos de Produção e Transporte Interno das Culturas de Banana e Soja 





0,84 
0,821 











3.2 Critérios técnicos adotados para o cálculo do sequestro de carbono 





0,84 
2,319 


Para quantificar as emissões de CO, decorrentes do consumo de óleo diesel no 


plantio, colheita e transporte interno da produção agrícola foram adotados os critérios 


sugeridos pela Convenção Internacional do Clima: “Diretrizes para Inventário Nacional 


de Inventários de Gases de Efeito Estufa”. 


De acordo com o método top-down do IPCC [1], a primeira tarefa é converter para 


uma unidade comum (TERA JOULE) a quantidade de emissões de gases de efeito estufa 


e o consumo de energia. Somente depois disso, as emissões de CO, podem ser medidas, 


usando a seguinte equação: 











CC=CA x Fconv x tEP x Fcorr 























Onde: 

1 tEP (Tonelada Equivalente de Óleo) = 45,2 x 103 
TJ (Tera-Joule) = 1012 J 

CC = Consumo de energia (TJ); 

CA = Consumo de combustível (mº, litros ou kg /t) 
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Fconv = Fator de conversão da unidade física de medida da quantidade de 
combustível para tEP, com base no maior poder calorífico (PCS) do combustível! 
Fcorr = Fator de correção PCS para ICP (menor valor calorífico)?. 


Nesta pesquisa, de acordo com o IPCC (Metodologia Simplificada), para a conversão 
para a unidade comum foi utilizado o valor calorífico inferior (PCI), o qual considera a 
energia do combustível efetivamente aproveitável. Portanto, foi adotada a seguinte fórmula: 


QC = CC x Femiss x 10º 





Onde: 

QC = teor de carbono expresso em GgC =tC x 10% 
CC = consumo de energia em TJ; 

Femiss = fator de emissão de carbono (tC / TJ)º 


Para estimar a quantidade de dióxido de carbono (CO,) e de carbono (C) disperso na 
atmosfera, é indicado o uso da razão de peso molecular, expressa pela seguinte equação: 








Onde: 
ECO, = emissão de CO, e EC = emissão de C. 


3.3 Alternativas Florestais para Mitigar as Emissões de Dióxido de Carbono 
(CO,) Parâmetros Técnicos 


As 04 (quatro) propostas selecionadas para o sequestro de carbono consideram as 
seguintes premissas: 


* Estimativas anuais de crescimento das árvores; 
* Espécies de árvores adequadas; 


* | Ciclos completos de desenvolvimento florestal conforme determinado pela fixa- 
ção de carbono; 


* Ciclos de sistemas agroflorestais (SAF's); e 


* | Dados sobre regeneração de floresta degradada com enriquecimento de flo- 
resta. 


Os dados do Projeto FLORAM [19] indicam que cada hectare de floresta nativa 


1 De acordo com o Balanço Energético Nacional, divulgado pela Empresa de Pesquisa Energética - EPE / 2014). O 
valor Fconv do óleo diesel em 2014 foi de 0,848456 tEP / mº 

2 No presente estudo, o fator de correção utilizado, que corresponde à relação PCS / ICP, foi de 0,939534884 

3 Segundo o IPCC (1996), o valor do fator de emissão de carbono para óleo diesel é de 20,2 tonC / TJ); 
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no clímax acumula em média entre 7,0 e 10,0 toneladas de carbono / ano. Por outro 
lado, as plantações de eucaliptos (2,5m x 2,5m, com 1.600 árvores), acumulam cerca de 
144 toneladas de C por hectare, em um ciclo de 06 a 07 anos. Essas taxas explicam a 
preferência dessa espécie na maioria dos projetos florestais. As sugestões de parcerias 
operacionais e um resumo de cada alternativa são apresentados a seguir: 

1º OPÇÃO: REFLORESTAMENTO com espécies de rápido Crescimento 

Essa alternativa de neutralização deve ser feita com as lavouras de eucalipto por 
meio do estabelecimento de parcerias com empresas que realizam reflorestamentos com 
fins econômicos. De acordo com este modelo, os talhões devem respeitar o espaçamento 
de 2,5 x 2,5m (1.600 árvores por hectare) [20] [21] [22], durante 04 (quatro) ciclos de 06 
(seis) anos cada. 

2º OPÇÃO: PROTEÇÃO FLORESTAL e PROJETOS de ENRIQUECIMENTO 

Para essa alternativa são recomendadas parcerias com Secretarias Estaduais e 
Municipais de Meio Ambiente, empresas do setor elétrico e ONGs. As áreas escolhidas 
devem circundar e proteger os remanescentes de mata ciliar [23] [24], as nascentes de 
água doce e as áreas degradadas das margens dos reservatórios. A ideia principal é 
proteger os recursos hídricos, o solo e a biodiversidade regional. É importante identificar os 
líderes locais e realizar atividades educacionais para envolver outras partes interessadas. 

3º OPÇÃO: IMPLEMENTAÇÃO de SISTEMAS AGRFLORESTAL (SAFs - 1º Ciclo) 

Esta alternativa sugere a adoção de Sistemas Agroflorestais (SAFs). É fundamental 
buscar parcerias com Secretarias Estaduais e Municipais de Agricultura e Meio Ambiente 
que, em parceria com Sindicatos de Trabalhadores Rurais e agências de assistência 
técnica especializada e extensão rural, possam auxiliar na escolha de quais assentamentos 
de produtores rurais e / ou familiares devem ser auxiliados. 

No Brasil, tais modelos SAFs podem ser modificados em experimentos de base 
regional realizados pela EMBRAPA. Os agricultores familiares seriam responsáveis pela 
manutenção dos projetos nos anos subsequentes, com o auxílio de agências de fomento 
agroflorestal [25- 26]. Como as outras opções, os SAFs exigem a conclusão dos ciclos para 
contribuir, mais efetivamente, para o sequestro de carbono. 

4º OPÇÃO: PROJETOS DE ARBORIZAÇÃO URBANA (RUAS E PARQUES) 

Sugere-se, neste caso, que as Secretarias Municipais de Meio Ambiente 
estabeleçam parcerias com empresas, ONGs e órgãos públicos em áreas urbanas. O ciclo 
de crescimento das espécies urbanas foi estimado em 25 a 30 anos [27] [32]. Utilizamos 
01 (um) hectare de referência para facilitar o cálculo. Assim, a quantidade de árvores a 
ser plantada em uma área retangular convencional é indicada de acordo com o espaço 
disponível. 
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Plantio de mudas de árvores do Centro de Referência em 
Recuperação Florestal (CR-ad UNTIVASF - Petrolina, Nordeste do 
Brasil) para REFLORES TAR ce RECUPERAR áreas degradadas 


SISTEMAS  AGROFLORESTAIS árvores são 
intercaladas por culturas agrícolas (BANANA, feijão, milho, café, 
etc.) produzindo alimentos e renda 


onde as 


ros agricultores familiares 


municipio de Delfinópolis, Estado de Minas Gerais, Brasil; 
mecaniza Foto: EMATER / MG. 





PROTEÇÃO FLORESTAL com para estimular a 
regeneração natural em áreas de preservação permanente (APP's: 
roteção de nascentes, margens de nos. vales, etc.) 


cercas 


Exemplo de PROJETO DE ARBORIZAÇÃO URBANA, em 


Brasília / DF, Brasil, para paisagismo c amenização climática 
(Cassia sp.) 


no Centro-Oc 


óleo diesel Foto: Daniel Popov, Canal Rural. 


A Tabela 3 apresenta uma decomposição de dados para cada opção, indicando 


ciclos, premissas de plantio de árvores, dados sobre SAFs e crescimento de espécies 


nativas, bem como os requisitos para arborização urbana. 
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Table3 Forestry technical assumptions considered for cach CO; emission mitigation alternative [27]. 


1st option: reforestation with fast growing species 

Species used 

Number of cycles for commercial use/years 

Area (m”) 

Spacing 

Number of trees/ha 

Annual C fixation by tree (ton C/tree) 

Total amount of fixed C (ton C/ha) 

Conversion factor for CO,/ha 

Total C captured in the complete cycle (ton C/ha) 

Technical projechon for emissions neutralization (number of ha) 
Total number of trees to be planted 

2nd option: forestry protecton project with ennchment 

(1 ounding degraded areas with forests at an initial stage of 
Reference area (m?) 

Percentage of forests in | ha 

Amount of C accumulated after four years (ton) 

Estimate of total natural regeneration until the forest reaches climax stage (ton C) 
Spacing 

Ennchment with forest species and fruit trees (number of trees) 
Total yearly estimated carbon (ton C/tree/ha) 

Number of native trees or planted fruit specres 

(2) Enclosure to protect native forest remnants and'or secondary forests 
Percentage of forests in | ha 

Amount of C accumulated/year (ton C) 
Estimated cycle until climax (years) 
Estimated C asset per year (ton) 

Estimation of the amount of accumulated carbon at climax (ton Cy 
Estimation of annual € fixation (ton/ha) 


Total C captured during entire cycle (ton C) 
Projected area to neutralize CO, emission (ha) 
3rd option (famuly farmers financing) implant: 
(1) Forest par 

Reference area (m?) 

Ga for forest management (ha) 

Forest plantation area (m?) 

Proposed spacing 
Number of arboreal species 

Estimate of annual € fixation (tontree/ha) 
Timber extraction cycle 
Total C captured during the comp 
(2) Agronomic part 
Reference area (m?) 

% for agricultural production (ha) 
Number of day cycles (sprouts and replanung) 
Proposed spacing 


om of agrofores! 


Eucalypeus sp. 

4cycles » 6 years = 24 
10,000 
25m=25m=625m 
1,600 = 10,0006.25 
0.0450 

2878 

44/12=3,67 

1,055.3 = 287.8 = 3.67 
19.0 = 20,000 (ton C/hay 1,055 (ton €) 
30,400 = 19 1,600 


10,000 
oss 

2783 

114.8 = [(27.83/4) x 0.55 x 30 years)] 
40m=45Sm= I80mº 

555.6 = 10,00018 

4.70 = (0.00840 ton'tree/ha = 556) 
4214 


04s 
1148 
ET 
130 
640 
670=[(2783 x 0.55) +(LI48 x 0.45)] 
2,639 = [(4.67 x 30 years » 3.67) + (67 x 30 years) 
+ 1148] 
7.6 = 20,000 (ton C0:)2,639 (ton COv/ha) 
le) 


10,000 

os 

5,000 
3S5m=4m=I4m' 
7143 = 10,000/14 
634 = (0008871 = 714) 
3 cveles = 6 years = 18.0 
418.20 


10.000 

os 

3,000 

no 
20m=20m=4m 
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(Table 3 continued) 


Number of banana trees 

Estimation of the amount of dry matter/ha in the Ist cycle (> 500 to < 600 d) 
(ton) 

Estimated amount of Cha (ton) 

Estimation of C fixation per banana tree (ton/ha) 

Amount of C fixed in AFS planting (example) (ton) 

Total C capture during the complete cycle (ton) 

Planting of grains + sugarcane 

Area of various crops to produce grains (com, beans) and others (m?) 
Number of cycles (years) 

Estimation of dry matter/ha/cycle (ton) 

Estimation of the fixation of C/ha by temporary crops (ton) 

Amount of C fixed in the present planting of the agroforestry system (ton) 
Total C capture during the complete cycle (ton) 

Total C captured by the agroforestry system (ton C) 

Technical project for emission neutralization (ha) 


Number of trees for forest management/ennichment (considering suggested 
cvcles) 


750 = 3000/4 
48.79 


244= (488 = 0.4 default IPCC) 
0.026 = (19.5/750 tree) 

234 = 0,026 750 = 12 

858.8 = 234 x 3.67 


2,000 
180 

75 

182=75x0242 
3267=1.82x 18 years 
119,79=32.67x 3.67 

13968 =4182+8588+ 1198 
1432 = 20,000/1,396.8 


10,232 = 14,32 = 714 


4th option: Urban afforestation project (streets and public parks) with native species 


Reference area (m*) 

Estimated full development cycle 

Planting space in streets and avenues 

Planting space in public parks 

Amount of trees per referential ha (planting in the streets) 
Amount of trees per referential ha (planting in squares) 
Total trees'reference ha 

Estimation of annual C fixation per tree (ton C) 
Estimate of the annual fixation of C/reference ha (ton C) 
Estimation of total carbon fixation (ton C) 

Conversion factor for CO; 

Total C capture in the complete cycle (ton C) 

Technical project for emission neutralization 

(number of reference ha) 

Total number of trees to be planted 


10,000 

| cycle x 25 years = 25 
6m=6m=36m 
25m=25m=625m 
10,000/36 = 278 
10,000/6.25 = 1,600 
278 + 1,600 = 1,878 
0.00849 

15.944 = 0.00849 = 1,878 
39832 

44/12=367 

1,460.5 


137 
25,714 


Sources of options 1, 2, 3 and 4: templates developed by the author, using data from Floram project and information's studies done 


by EMBRAPA and Brazilian universities 





41 RESULTADOS E DISCUSSÕES 


As alternativas de sequestro de carbono apresentadas neste modelo matemático 
(Opções: 1, 2, 3 e 4) visam cobrir os diferentes aspectos da ciência florestal, desde o 
reflorestamento com espécies de rápido crescimento até projetos agroflorestais. Em todos 
eles há um desejo explícito de estimular o plantio de árvores para atender a demanda por 
madeira e outros produtos; bem como proteger: solo, biodiversidade e recursos hídricos, 
como forma de compensar as emissões causadas pelo consumo de óleo diesel [23] [24]. 

A maior dificuldade encontrada em termos de captura de carbono foi, sem dúvida, em 
relação à alternativa que trata da proteção e enriquecimento de remanescentes florestais 
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(Opção 2 — recuperação de áreas verdes), uma vez que cada área a ser recuperada, devido 
aos diferentes níveis de degradação, requer plantios de espécies nativas em quantidades 
muito variáveis [28] [29]. 

Apesar disso, a contribuição da regeneração natural não deve ser negligenciada, 
especialmente em regiões tropicais e subtropicais onde a alta incidência de luz solar 
combinada com outros fatores edafoclimáticos favorecem o rápido crescimento da 
cobertura vegetal. Portanto, considerou-se que determinadas áreas, com os plantios de 
enriquecimento, são uma importante ação para mitigar parte das emissões de COZ. 

Por outro lado, as opções 1 e 4, de reflorestamento e arborização urbana, que 
indicam, respectivamente, plantações de eucalipto e de espécies nativas ou exóticas que 
deveriam atingir a idade de clímax aos 25 anos, foram as de menor grau de complexidade. 

O Eucaliptus, por exemplo, possui diversos estudos de desenvolvimento florestal 
realizados no Brasil, o que facilitou a determinação de seus parâmetros. 

Por sua vez, pelo fato de a arborização urbana das cidades brasileiras utilizar um 
número reduzido de espécies florestais nativas, os dados médios de crescimento arbóreo 
foram inferidos com o apoio da bibliografia consultada. 

A opção 4, que propõe a implantação de sistemas agroflorestais, teve o proósito 
de estimular os pequenos produtores rurais a fazerem plantios integrados de alimentos 
e árvores para atender às suas necessidades. Nesse caso, o modelo apresenta um 
exemplo de associação de espécies de 1,0 hectares, sendo composto por: 0,5 ha com 
espécies florestais, 0,3 ha com banana; e 0,2 ha com cana-de-açúcar e grãos). Levando 
em conta um ciclo de 18 anos, estima-se que um montante de 860,0 toneladas de carbono, 
aproximadamente, podem ser sequestrados. Porém, esse resultado se aplica apenas a esta 
situação específica, pois a quantidade de carbono a ser capturada dependerá diretamente 
da composição do SAF. Mas, presume-se que a participação dos pequenos agricultores em 
larga escala tende a ser relevante, contribuindo para a segurança alimentar e o equilíbrio 
ambiental. 

Na tabela 3 verifica-se, detalhadamente, as informações básicas de cada uma das 
alternativas sugeridas. Os dados de crescimento das espécies florestais foram obtidos no 
Projeto FLORAM (USP). Os dados de emissões de COZ foram elaborados pela equipe 
IVIG/ COPPE / UFRJ, após a pesquisa bibliográfica, incluindo: dados de produção agrícola 
do IBGE, estudos da EMBRAPA e consultas na internet sobre os métodos de produtos 
agrícolas que fazem parte do presente modelo matemático (Tabela 1). 

As Tabelas 4A / 4B exemplificam os resultados das produções brasileiras de 
BANANA e SOJA, entre 2000 e 2012; a tabela 5 apresenta as sugestões de mitigação 
das respectivas emissões. Enquanto, a tabela 6 resume os resultados das alternativas de 
mitigação de todas as culturas agrícolas: permanentes e temporárias, realçando o fato da 
banana e da soja ocuparem os extremos do ranking, tanto no Brasil, como em todos os 


demais países da América Latina. 
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Por fim, a tabela 7 apresenta o ranking que mostra, em ordem decrescente, as 


culturas que mais consumiram óleo diesel e emitiram CO2 para a atmosfera, com a 


respectiva conversão para consumo em toneladas equivalentes de óleo (tEP) durante os 


13 anos. 


























TABELA 4A - BANANA: RESULTADOS DA PRODUÇÃO 2000 2001 2011 2012 Total 
PRODUÇÃO - Carga transportada por ano (Mil toneladas) [4] [15] [32] 7.451,5 6.838,7 7.262,9 6.786,7 90.580,5 
Média da Produção Anual [15] 6.967,7 
DADOS DE MECANIZAÇÃO e de TRATAMENTOS CULTURAIS 
Nº de hectares com Mecanização (Plantios, Manutenção e colheita - 1000 ha) 47,23 51,03 75,86 78,53 
Média Estimativa do Consumo de Óleo Diesel (Litros por hectare) 44,0 44,0 44,0 44,0 
Consumo Total de Óleo Diesel na Mecanização (Litros) 2.078.010,0 2.245.377,2 3.337.620,0 3.455.100,0 35.758.125,4 
TRANSPORTE DA ÁREA DE PLANTIO ATÉ O DEPÓSITO DE DISTRIBUIÇÃO 
Distância das Viagens de Transporte - ida e volta (km) 133.062 122.120 129.694 121.191 
Nº de Viagens por ano 532 488 519 485 
Carga Transportada por dia (1000 toneladas) 36,53 33,52 35,60 33,27 
Consumo Total de Óleo Diesel da Produção (Litros) 33.265,4 30.529,9 32.423,5 30.297,8 404.377,3 
EL ENEI Cc IM N TE 
Estimativa do Consumo Total de Óleo Diesel (Litros / ano) 2.111.275,4 2.275.907,1 3.370.043,5 3.485.397,8 36.162.502,7 
Estimativa do Consumo Total de Bieo Diesel (m* / ano) 2.111,28 2.275,91 3.370,04 3.485,40 36.162,5 
Estimativa do Consumo Total de Oleo Diesel toneladas / ano) 1.773,5 1.911,0 2.830,8 2.927,7 30.376,5 
Consumo de Oleo Diesel (Tonelada Equivalente de Petróleo - tEP) 1.504,7 1.621,2 2.400,5 2.482,7 25.760,/1 
ESTIMATIVA DAS EMI ES DE COZ (toneladas) 5.215,7 5 8.325,35 8.610,32 89.335,8 
MÉDIA GERAL DAS EMISSÕES DE COZ 13 anos - toneladas] 6.872,0 
CONSUMO TOTAL DE ENERGIA TJ) 1.206,2 


TABELA 4B - SOJA: RESULTADOS DA PRODUÇÃO 2000 


2001 


2011 2012 





PRODUÇÃO — Carga transportada por ano (Mil toneladas) [4] [15] [32] 38.823,61 
[15] 


Média da Produção Anual 


DADOS DE MECANIZAÇÃO e de TRATAMENTOS CULTURAIS 

Nº de hectares com Mecanização (Plantios, Manutenção e colheita - 1000 ha) 
Média Estimativa do Consumo de Óleo Diesel (Litros por hectare) 118 
Consumo Total de Óleo Diesel na Mecanização (Litros) 1.429.619.878,8 


TRANSPORTE DO PLANTIO ATÉ A INDÚSTRIA DE BENEFICIAMENTO 

Distância das Viagens de Transporte - ida e volta (km) 1.094.120 
Nº de Viagens por ano 186.000.451 
Carga Transportada por dia 241.350,07 
Consumo Total de Óleo Diesel da Produção 66.428.732,5 


12.324,31 


(1000 toneladas) 
(Litros) 


37.905,72 


12.589,52 
118 
1.460.383.844,4 


1.263.524 
214.799.059 
278.718,50 
76.713.949,6 


74.815,46 81.486,8 


21.647,17 24.949,35 
118 118 
.511.071.604,0 |2.894.124.600,0 


2.493.849 
461.758.756 
550.113,67 
151.412.237,6 


2.716.227 


599.167,68 
164.913.770,0 


28.020.466.083,6 


1.453.270.533,5 





DIES nd G ONS O 


Estimativa do Consumo Total de Óleo Diesel Litros / ano) 1.496.048.611,3 


1.537.097.794,0 


« |2.662.483.841,6 |3.059.038.370,0 


29.473.736.617,1 





Estimativa do Consumo Total de Óleo Diesel 'm?/ ano) 1.496.048,61 


1.537.097,8 


3.059.038,4 


29.473.736,6 





Estimativa do Consumo Total de Oleo Diesel (toneladas / ano) 1.256.680,8 
Consumo de Óleo Diesel onelada Equivalente de Petróleo - tEP) 1.066.238,4 
ESTIMATIVA DAS EMISSÕES DE COZ Ç 
MÉDIA GERAL DAS EMISSÕES DE COZ 





1.291.162,1 
1.095.494,3 








2.569.592,2 
2.180.185,9 





(13 anos - toneladas) 
J) 


CONSUMO TOTAL DE ENERGIA 
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24.757.938,8 


PLANTIOS 
CULTIVOS AGRÍCOLAS / ALTERNATIVAS DE SEQUESTRO DE CARBONO | 2000 a om | 2 | FLORESTAIS 
[e [| 
23 24 








* OPÇÃO: REFORESTAMENTO COM EUCALIPTUS SPP. EE 


* PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES (Nº de ha) Ea [E Nº de ha. 2,0 
BANANA Nº de 


& Nº DE ÁRVORES PLANTADAS ape Árvores 


% PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES (Nº de ha) 1.038,0 |..|.. Nº de ha. 
SOJA Nº de 


% Nº DE ÁRVORES PLANTADAS 1.660.873 | .:]/ie |, 2876859 Árvores 
P PROJETOS DE PROTEÇÃO FLORESTAL + plantios de enriquecimento 


* PROJEÇÃO: Número de hectares para ser isolado com plantios de 
BANANA recuperação florestal com árvores nativas 
46 5,0 coli 7,4 7,6 Nº de 








2.449.779 


% PROJEÇÃO: Número de hectares para ser isolado com plantios de hectares 
SOJA recuperação florestal com árvores nativas 


3.283,6 3.373,77 |. 5.843,7 6.714,1 


PÇÃO IMPLANTAÇÃO DE SISTEMAS AGROFLORESTAIS (1º Ciclo) 
Nº de 
% PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES (Nº de ha) o hectares 
Nº de 
* Nº DE ÁRVORES PLANTADAS a li Árvores 


Produtores 
* Nº PRODUTORES BENEFICIADOS (02 famílias / ha) e lie Beneficiados 


Nº de 
1.239,8 1.273,9 |. |. 2.206,5 2.535,2 
hectares 


Nº de 
& Nº DE ÁRVORES PLANTADAS 885.606 | 909.905 |..|.. | 1.576.092 | 1.810.838 ênúres: 


Produtores 
% Nº PRODUTORES BENEFICIADOS (02 famílias / ha io Egas E oa cad 5.070 | peneficiados 


[4º OPÇÃO | PROJETOS DE ARBORIZAÇÃO URBANA 
* PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES Nº ha de 
(Nº hectares de Referência) Eai [ss referência 
Nº de 
& Nº DE ÁRVORES PLANTADAS aailtas Árvores RAS 


* PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES Nº ha de 7.473 
(Nº hectares de Referência) 750 sã io 1.299,2 referência ” 
2.077.482 








* PROJETO DE NEUTRALIZAÇÃO DAS EMISSÕES (Nº de ha 
SOJA 








Nº de 











* Nº DE ÁRVORES PLANTADAS 1.408.467 ..|.. 2.439.677 Árvores 


TABELA 5: Opções sugeridas para mitigar as emissões de dióxido de carbono (CO,) / tabela 
comparativa entre as safras de banana e soja 
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TABELA 6: ALTERNATIVAS PARA MITIGAR AS EMISSÕES DE CO2, 
POR TIPO DE CULTIVO - SAFRAS PERMANENTES 


Média das Toneladas de 
Emissões de Óleos OPÇÃO 1 OPÇÃO2 | OPÇÃO3 OPÇÃO 4 
coz Equivalentes 
Nº 
(em 13 anos) Nº de Nº de Nº de Nº de Nº de Produtores 
Árvores Hectares Árvores Hectares Árvores —Beneficiados 
BANANA 6.872,0 1.981,5 3.004,0 2.548,0 1.646,0 


Café e Cacau 156.234,4 45.049,0 68.299,0 34 57.919,0 J 34417,0 





Cõco e Dendê 43.703,3 11.946,8 19.105,0 16.202,0 s 10.467,0 


Frutas Citricas 282.222,0 B1.376,5 123.374,00 | 104.625,0 67.590,0 
Frutas de Clima 
Tropicais 135.658,6 39.141,7 59.343,0 50.324,0 32.511,0 
Frutas de Clima 
Sub-Tropical 313.171,3 90.300,5 136.904,0 5, 116.099,0 75.003,0 
Uva, Maracujá 
e Pimenta do 34.245,6 9,874,5 14.9740 ol 12.696,0 s 8.202,00 
Reino 
Extrativismo 
vegetal 13.719,7 4.709,5 7140,0 ç 055, k. 3.912,0 
(palmeiras) 
Extrativismo 
egetal (Produtos) 30.180,7 10,359,9 15.707,0 ; 13.320,0 20, 8.605,0 
N- Madeireiros e 
iadeiras nobres 














SUBTOTAL 1.016.007,6 294,739,8 447.847,0 379.786,0 245.351,0 
Permanentes 

















CULTIVOS TEMPORÁRIOS (agronegócio) 
Média das Toneladas J E 
Emissões de de Óleos OPÇÃO 1 OPÇÃO 2 
coz Equivalentes 
Nº 
(em 13 anos) Nº de Nº de Nº de Nº de Nº de Produtores 
Árvores Hectares Árvores Hectares Árvores Beneficiados 


Algodão 

(Herbáceo e 310.668,7 89.578,9 135.910,0 8s, 115.171,0 178,0 74.403 208,0 
Arbóreo 

Açúcar [7] 3.570.142,7 | 1.029.976,0 | 1.561.541,0 976,0 | 1.324.231,0 2046,6 855.486 2.395,0 


(Temporários) 





a compensação pela produção de biodiesel e etanol anidro não foi deduzida dos percentuais de 
consumo de óleo diesel. Esses temas não fizeram parte dos objetivos desta pesquisa. 
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[*] No cálculo das opções para neutralizar as emissões de CO2 das lavouras de soja e cana-de-açúcar, 


Cana-de-Açúcar 3.570.142,7 1.029.976,0 
Frutas de Clima Sub-Tropical 313.171,3 90.300,5 
Algodão (Herbáceo e Arbóreo) 310.668,7 89.578,9 
Frutas Cítricas 282.222,0 81.376,5 
Café e Cacau 156.234,4 45.049,0 
Frutas de Clima Tropical 135.658,6 39.141,7 
Côco e Dendê 43.703,3 11.946,8 
Uva, Maracujá e Pimenta do Reino 34.245,6 9.874,5 


Extrativismo vegetal (Palmeiras): 30.180,7 10.359,9 


Extrativismo vegetal (Produtos N- Madeireiros e 13.719,7 4.709,5 
Madeiras nobres 














(*) Cultivos Temporários 


51 CONCLUSÃO 


O modelo matemático, em relação ao consumo de óleo diesel e emissão de dióxido 
de carbono (CO,), mostra - numericamente e comparativamente - a grande supremacia do 
cultivo da SOJA, em relação às demais plantas pesquisadas, visto que o uso de maquinário 
agrícola é muito elevado em todas as etapas do processo de produção. Em vez disso, 
nas plantações de BANANA, o processo de mecanização é muito baixo, tanto no Brasil 
quanto nos países do hemisfério sul, e pouco contribui para as emissões deste gás de 
efeito estufa, pois ainda há uma grande participação do trabalho braçal (energia humana). 
De acordo com a tabela 6, o cultivo da BANANA, no Brasil, durante 13 anos (2000 a 2012), 
lançou para a atmosfera, apenas 6,90 mil toneladas de CO,; enquanto a SOJA 5,60 milhões 
de toneladas. 

Este estudo, s.m.j., incentiva a política de agricultura de baixo carbono (ABC), 
elaborada pela EMBRAPA - [30][31]). Outrossim, dá subsídios aos debates voltados para 
a melhoria dos métodos de coleta de dados do setor agroflorestal. Devido a ausência de 
alguns dados básicos e de perceptíveis inconsistências, muitas inferências precisarm 
ser feitas [32] com o suporte das pesquisas bibliográficas mais amplas. Isto reforça a 
necessidade de maior recursos e apoio operacional do governo para a urgente adequação 
das pesquisas de campo, mesmo reconhecendo os esforço dos técnicos do IBGE que 
participam dos Censos Agropecuários. 
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Apesar disso, o modelo matemático apresenta uma síntese consistente da produção 
agrícola brasileira, associada aos levantamentos dos parâmetros determinados pela 
Ciência Florestal que são adotados nas diferentes regiões do país (Tabelas: 3 e 4). Ele 
fornece, também, informações relevantes sobre a fixação de carbono em árvores tropicais; 
bem como, as estimativas a respeito dos ciclos de crescimento das mesmas, as quais 
serviram de base para as opções de mitigação das emissões de CO, sugeridas e podem 
contribuir para as futuras decisões a serem tomadas no mercado de carbono [31]. 

Por fim, conclui que caso as propostas de plantios de árvores para sequestro de 
carbono sejam adotadas e/ou transformadas em políticas públicas, muitos benefícios 


adicionais seriam, direta ou indiretamente, alcançados: 
* Geração de trabalho e renda no meio rural; 
* Proteção do solo, água e biodiversidade; 
* Melhoria da arborização urbana; e 


* Estímulo à implantação de sistemas agroflorestais SAFs, contribuindo para a 
segurança alimentar de pequenos produtores rurais. 
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Fonte: Ministério do Meio Ambiente (MMA). 
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RESUMEN: El objetivo del presente estudio fue 
de Cuantificar el valor económico del potencial 
forestal en 4 concesiones mineras de la cuenca 
del río Madre de Dios. En cada concesión se 
estableció una parcela de 100 m x 10 m (0.1 
ha) y se evaluaron todos los individuos con un 
diâmetro a la altura del pecho mayor a 10 cm 
(DAP > 10 cm). Para la valorización del potencial 
forestal de las áreas evaluadas se utilizará el 
valor al estado natural (VEN) de la madera de 
rolliza (soles/m3), de todos los individuos y 
especies que se registrarán en los inventarios, de 
acuerdo a la “resolución de dirección de dirección 
ejecutiva Nº 241-2016-SERFOR-DE”, que indica 
el valor comercial de la madera al estado natural 
en nuevos soles por metro cubico”. La riqueza de 
especies (7-22 especies) encontradas fue inferior 
al reportado en otros en Madre de Dios. Se calculó 
que el valor de la madera al estado natural oscila 
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entre S/. 216/ha y S/. 1305/ha. Las diferencias en 
los valores al estado natural de la madera entre 
las concesiones estudiadas están influenciadas 
por la presencia de algunas especies con un 
elevado valor económico (categorias By C) y al 
elevado grado de extracción maderera en la que 
tienen valores más bajos. 

PALABRAS CLAVE: Riqueza de especies, valor 
al estado natural, madera. 


ECONOMIC VALUATION OF FORESTRY 
POTENTIAL IN MINING CONCESSIONS 
IN THE SOUTHERN PERUVIAN AMAZON 
ABSTRACT: The objective of this study was to 
quantify the economic value of forest potential 
in 4 mining concessions of the Madre de Dios 
basin. In each concession, a plot of 100m x 10m 
(0.1 ha) was established, and all individuals with 
a diameter at breast height greater than 10 cm 
(DBH>10 cm) were evaluated. For the valorization 
of the forest potential of the evaluated areas, the 
value to the natural state (VEN) of the round wood 
(soles / m3), of all the individuals and species that 
will be registered in the inventories, will be used, 
according to the “resolution executive direction 
address No. 241-2016-SERFOR-DE “, which 
indicates the commercial value of wood to natural 
state in new soles per cubic meter”. The species 
richness (7-22 species) found was lower than 
that reported in others in Madre de Dios. It was 
calculated that the value of wood to the natural 
state ranges from 216 soles/ha and 1305 soles/ 
ha. The differences in the values to the natural 
state of the wood between the concessions 
studied are influenced by the presence of some 
species with a high economic value (categories B 
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and C) and the high degree of logging in which they have lower values. 
KEYWORDS: Species richness, value to the natural state, wood. 


11 INTRODUCCIÓN 


Los bosques tropicales de la Amazonía son como uno de los más biodiversos del 
planeta, debido a que alberga una enorme riqueza de especies y recursos genéticos. 
Estos bosques brindan servicios ambientales de suma importancia a escala global, como 
la conservación de suelos y captura y proceso de CO, (CAYUELA e GRANZOW, 2012; 
WRIGHT, 2010). Por otro lado, se tiene conocimiento que más del 30% de sus especies 
son endémicas, con ecosistemas bastante frágiles, que vienen siendo sometidos por afios 
a actividades degradantes, como son: la extracción de madera, la minería y ampliación de 
la frontera agropecuaria (GOREMAD-IIAP, 2008). 

Un interesante enfoque basado en los servicios ecosistémicos es el uso de técnicas 
de valoración económica, el cual cuantifica los servicios que proveen los ecosistemas 
en valores monetarios (RESENDE et al., 2017), Aunque suene complicado asignarle un 
valor monetario a bienes y servicios que provee la naturaleza (LEÓN, 2007). El término 
“valoración económica” surgió a partir de la premisa de compensación de los dafios socio- 
ambientales, utilizando instrumentos de gestión ambiental adecuados, además de estar 
promocionado por el programa de “adaptación al Cambio Climático de la Organización de 
las Naciones Unidas” (LEÓN, 2007). 

Según Recavarren, Delgado, Angulo, León, y Castro (2011), está proyectado 
que para el afio 2030 un aproximado de 994 910 ha hayan sido deforestadas. Siendo 
el incremento poblacional, carretera interoceánica y la minería los principales agentes. 
Asimismo, el incremento poblacional generará una mayor presión por bienes y servicios, 
ocasionando mayores tasas de deforestación (disminución de sumideros de carbono). 

Por otro lado, se tiene conocimiento que la actividad minera aurífera en la región 
Madre de Dios está afectando severamente en todos los componentes de los ecosistemas 
forestales, como flora, fauna, suelo, servicios ecosistémicos, alterando la topografia e 
hidrología. Esto genera un cambio en la condición primaria de los bosques tropicales de 
Madre de Dios, debido a que estos también generan cambios importantes en la estructura 
edáfica por el uso de maquinaria pesada, que ocasiona la compactación de los suelos 
(BARRERA-CATANO e VALDÉS-LÓPEZ, 2007). 

Asimismo, se tiene conocimiento que, en las concesiones mineras del departamento 
de Madre de Dios, existe insuficiente información acerca de la composición florística, 
estructura horizontal y el valor económico de las especies maderables, lo cual dificulta su 
manejo y aprovechamiento sostenible. La importancia de realizar la presente investigación 
se fundamente en que: (1) servirá de ejemplo en la valoración económica de estos 


ecosistemas. (2) una vez cuantificado el valor económico de las especies forestales se 
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podría establecer políticas de estado que permita al concesionario minero realizar un pago 
por derecho de aprovechamiento de estas especies forestales de importancia económica. 

Por otro lado, el estado peruano ha identificado a los factores sociales, económicos 
e institucionales como fuente de degradación de los ecosistemas boscosos del país. 
Asimismo, el estado como estrategia prioriza el fomento del aprovechamiento sostenible 
y conservación de los recursos forestales y de fauna silvestre, debido a su potencial 
contribución a la inclusión a la economía peruana y a la adaptación y mitigación del cambio 
climático. En este sentido, el estado peruano mediante un sistema de áreas naturales 
protegidas, comunidades nativas y áreas destinadas a producción forestal contribuyen en 
la conservación y sostenibilidad de los bosques. 

Sin embargo, se tiene identificado una debilidad regional y estatal en lo referido al 
ordenamiento territorial, debilidad en el cumplimiento legislativo y falta de coordinación 
entre diferentes entidades estatales. En este contexto, se entorpece la aplicación de 
instrumentos de control sobre los factores socioeconómicos que degradan los bosques 
de Madre de Dios. El objetivo del presente estudio fue cuantificar el valor económico 
del potencial forestal que poseen los bosques comprendidos en 4 concesiones mineras 
de la cuenca del río Madre de Dios (Perú), otorgados por el estado peruano con fines 
de explotación de minerales no metálicos. Debido a que, estos bosques están siendo 
aprovechados sin contar con la información sobre su composición florística, diversidad 


arbórea y su importancia económica. 


21 MÉTODOS 


Área de estudio 


El estudio se desarrolló en 4 concesiones mineras de la cuenca de Madre de Dios, 
en la Amazonia sur del Perú (Figura 1). 
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Figura 1. Ubicación de las 4 concesiones mineras evaluadas en el presente estudio. 


Inventario de especies arbóreas en las concesiones mineras 


Para el inventario de las especies arbóreas se utilizará una parcela de 100 m x 
10 m (0.1 ha), en la que se registrará a todos los árboles con un diámetro a la altura del 
pecho mayor a 10 cm (DAP > 10 cm). Utilizando un formato de campo se registrará la 
especie, nombre común, DAP y altura; de cada uno de los individuos registrados durante 


el inventario. 


100 m 


A, AAA, 
PII AI: 
47 DP 9% A, 57 


a 
PR 7 


PEDI DDDA: 





Figura 2. Disefio de la parcela utilizada para la evaluación del potencial forestal de 4 concesiones 
mineras en Madre de Dios. 
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El área basal de cada árbol se obtendrá mediante la fórmula 1: 


(DJ? 





[1] 


Dónde: 

AB = Área Basal del árbol (m?) 
m=3.1416 (Constante) 

D = Diâmetro. 


El volumen de cada árbol, para estimar el potencial forestal maderable, se obtendrá 
mediante la fórmula 2: 


V=ABx+Hfx+Ff [2] 


Donde: 

V = Volumen del árbol (mº) 

AB = Área Basal (m?) 

Hf = Altura de fuste (m), distancia desde la base del árbol hasta el inicio de la primera 
ramificación. 

Ff = Factor de forma (0.65) para bosques tropicales. 


Con la finalidad de calcular los volúmenes maderables en el área implicada, se 


evaluará cada una de ellas, determinando sus áreas. x 


Valoración del potencial forestal del bosque 


Para la valorización del potencial forestal de las áreas evaluadas se utilizará el 
valor al estado natural (VEN) de la madera de rolliza (soles/m3), de todos los individuos y 
especies que se registrarán en los inventarios, de acuerdo a la “resolución de dirección de 
dirección ejecutiva Nº 241-2016-SERFOR-DE”, que indica el valor comercial de la madera 
al estado natural en nuevos soles por metro cubico” según especie (Tabla 1). 








Categoría s o Denominación 
A 55 Altamente valiosas 
B 12 Valiosas 
Cc 6 Intermedias 
D 4 Potenciales 
E 2 Otras especies (Valor económico futuro) 





Tabla 1. Categorización del valor de la madera al estado natural de acuerdo a la Resolución Nº 
241-2016-SERFOR-DE. VEN = Valores al Estado Natural de la Madera. 
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31 RESULTADOS Y DISCUSIÓN 


Riqueza de especies según concesión 


En las Tablas 2-5 se muestra el listado de especies y el volumen aportado por cada 
especie en el presente estudio. En las concesiones “Mary Elena Erica” y “Paolita II” se 
encontró la mayor riqueza de especies, con 22 y 20 especies, respectivamente. Mientras 
que en las concesiones “E&0-2” y “Paolita |” se registró una menor riqueza de especies, 
con 10 y 7 especies, respectivamente. Sin embargo, la riqueza de especies encontrada es 
inferior al reportado para Madre de Dios por otros estudios (DUENAS e GARATE, 2018; 
FAUSET et al., 2015; PITMAN et al., 2003; TER STEEGE et al., 2013). 

Considerando las 4 concesiones, Ficus insipida (Paolita |; 155,3 m3/ha) fue la 
especie con mayor aporte en volumen (Tabla 2), seguida por Sapium marmieri (Paolita |; 
56,6 m3/ha) y Brosimum lactescens (Paolita Il; 48.9 m3/ha). 

En “Mary Elena Erica”, las especies con mayor aporte en volumen fueron las 
pioneras Acacia loretensis y Ochroma pyramidale (Tabla 1). En las concesiones “E&O-2” 
y “Paolita |”, las especies con mayor aporte en volumen fueron Ficus insípida y Sapium 
marmieri, seguido por las pioneras Erytrina poeppigiana, Ochroma pyramidale y Cecropia 
membranacea (Tabla 1,2). Mientras que en Paolita Il, las especies con mayor aporte en 
volumen fueron Brosimum lactescens, Ceiba pentandra y Calycophyllum spruceanum. 








Especie Volumen (mº/ha) 
Acacia loretensis 20.8 
Ochroma pyramidale 14.8 
Ficus insipida 11.5 
Erythrina ulei 8.2 
Sapium marmierii 2.0 
Cecropia membranacea 1.2 
Hyeronima alchorneoides 0.9 
Cordia alliodora 0.9 
Terminalia oblonga 0.9 
Otoba parvifolia 0.8 
Machaerium sp 0.8 
lriartea deltoidea 0.7 
Pleurothyrium cuneifolium 0.6 
Mabea nitida 0.6 
Dendropanax arboreus 0.6 
Pseudolmedia laevis 0.6 
Bixa excelsa 0.5 
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Urera caracasana 0.5 


Virola calophylla 0.4 
Casearia pitumba 0.4 
Psychotria sp 0.4 
Tabernaemontana cymosa 0.4 





Tabla 2. Volumen (m3/ha) de las especies registradas en la evaluación de la concesión “Maria Elena 








Erica” 
Especie Volumen (mº/ha) 

Ficus insipida 42.5 
Erytrina poeppigiana 28.3 
Ochroma pyramidale 24.1 
Cecropia membranacea 22.2 
Cordia alliodora 18.6 
Schizolobium parahyba 10.6 
Sapium marmierii 3.1 
Muntingia calabura 1.6 
Trema micrantha 0.9 
Heliocarpus americanus 0.3 


Tabla 3. Volumen (m3/ha) de las especies registradas en la evaluación de la concesión “E&O-2” 





Especie Volumen (mº/ha) 

Ficus insipida 155.3 
Sapium marmieri 59.6 
Cecropia membranacea 22.1 
Inga oerstediana 19.6 
Guazuma crinita 15.4 
Acacia loretensis 15.1 
Jacaratia digitata 7.9 


Tabla 4. Volumen (m3/ha) de las especies registradas en la evaluación de la concesión “Paolita |” 








Especie Volumen (m3/ha) 
Brosimum lactescens 48.9 
Ceiba pentandra 32.5 
Calycophyllum spruceanum 27.1 
Spondias mombin 17.8 
Tapura juruana 13.9 
Dipteryx micrantha 13.4 
Terminalia oblonga 11.9 
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Sloanea guianensis 11.3 


Brosimum alicastrum 11.1 
Pouteria ephedrantha 9.9 
Pseudolmedia laevis 7.8 
Ficus yoponensis 5.6 
Pouteria krukovii 5.2 
Hyeronima alchorneoides 4.8 
Luehea cymulosa 3.8 
Zanthoxylum juniperinum 3.6 
Trichilia ado!fii 3.6 
Apeiba membranacea 1.3 
Sapium marmieri 0.8 
Guazuma crinita 0.2 





Tabla 5. Volumen (m3/ha) de las especies registradas en la evaluación de la concesión “Paolita II” 


Valoración económica de las especies comerciales en 4 concesiones mineras 
de la cuenca del río Madre de Dios 

En la Figura 1 se observa los valores al estado natural de la madera estimado de 
acuerdo a la Resolución Nº 241-2016-SERFOR-DE. Se encontró un menor valor al estado 
natural de la madera en “Mary Elena Erica” (S/. 216) y “E&O-2” (S/. 597), mientras que 
“Paolita II” (S/. 1305) obtuvo la mejor valoración de la madera al estado natural. 


nam 1305.6 
1200 
1000 Rei 
< 
E 
õ 800 
(ão) 597.1 
ui 600 
E 
400 
216.4 
. E 
0 
Mary Elena Erica E&O-2 Paolita 1 Paolita 2 


CONCESIONES 


Figura 3. Valor al estado natural de la madera de especies forestales en 4 concesiones mineras de 
Madre de Dios. 
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Valor al estado natural de la madera (soles/ha) por especie 


En la concesión Mary Elena Erica, las especies con mayor valor natural fueron 
Acacia loretensis (38.4 %) y Ficus insípida (21.3 %); seguido por algunas especies pioneras 
como Ochroma pyramidale y Erythrina ulei. Asimismo, las 5 especies con mayor valor de la 
madera al estado natural representaron el 85.8% (Tabla 6). 








Rank Especie Categoría VEN (Soles S/. /ha) % 
1 Acacia loretensis D 83.08 38.4% 
2 Ficus insipida D 46.12 21.3% 
3 Ochroma E 29.52 13.6% 
pyramidale 

4 Erythrina ulei E 16.46 7.6% 
5 Cordia alliodora B 10.56 4.9% 

Otras 30.66 14.2% 





Tabla 6. Las 5 especies con mayor valoración de la madera al estado natural (VEN) en la concesión 
minera Mary Elena Erica. 


En la concesión E&O-2, las especies con mayor valor natural fueron Cordia alliodora 
(37.5 %) y Ficus insípida (28.4 %); seguido por algunas especies pioneras como Ochroma 
pyramidaley Erythrina poeppigiana. Asimismo, las 5 especies con mayor valor de la madera 
al estado natural representaron el 90.9% (Tabla 7). 








Rank Especie Categoria VEN % 
1 Cordia alliodora B 223.7 37.5% 
2 Ficus insipida D 169.9 28.4% 
3 Erytrina poeppigiana E 56.7 9.5% 
4 Ochroma pyramidale E 48.2 8.1% 
5 Cecropia membranacea E 44.4 74% 
Otras 54.3 9.1% 





Tabla 7. Las 5 especies con mayor valoración de la madera al estado natural (VEN) en la concesión 
minera E&O-2. 
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En la concesión Paolita |, las especies con mayor valor natural fueron Ficus insipida 
(64.6 %) y Sapium marmieri (12.4 %); seguido por algunas especies pioneras como 
Cecropia membranacea. Asimismo, las 5 especies con mayor valor de la madera al estado 
natural representaron el 94.3 % (Tabla 8). 








Rank Etiquetas de fila Categoria VEN o 
1 Ficus insipida D 621.0 64.6% 
2 Sapium marmieri E 119.2 12.4% 
3 Guazuma crinita D 61.5 6.4% 
4 Acacia loretensis D 60.5 6.3% 
5 Cecropia membranacea E 44.2 4.6% 
Total 55.1 5.7% 





Tabla 8. Las 5 especies con mayor valoración de la madera al estado natural (VEN) en la concesión 
minera Paolita |. 


En la concesión Paolita Il, las especies con mayor valor natural fueron Calycophyllum 
spruceanum (25.4 %), Brosimum lactescens (22.5 %) y Dipteryx micrantha (12.3 %); 
resultados similares fueron encontrados por Perez (2013). Asimismo, las 5 especies con 
mayor valor de la madera al estado natural representaron el 75.4 % (Tabla 9). 








Rank Etiquetas de fila Categoría valor % 
1 Calycophyllum spruceanum B 332.1 25.4% 
2 Brosimum lactescens Cc 293.6 22.5% 
3 Dipteryx micrantha B 161.1 12.3% 
4 Ceiba pentandra D 129.9 9.9% 
5 Sloanea guianensis Cc 67.5 5.2% 
321.5 24.6% 





Tabla 9. Las 5 especies con mayor valoración de la madera al estado natural (VEN) en la concesión 
minera Paolita Il. 


Proporción de especies según categorización del Valor al estado natural de la 
madera (soles/ha) 


En el presente estudio, no encontramos especias consideradas altamente valiosas 
(Categoría A) (Tabla 10). En las 4 concesiones se encontró diferencias en la proporción de 
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especies, según la categorización del valor de la madera al estado natural. Así, la mayor 
proporción de especies categorizadas como Valiosas e Intermedias (categorías By C) se 
encontró en la “Paolita II” (Tabla 10). Mientras que en “Mary Elena Erica” se encontró una 
menor proporción de especias valiosas. Sin embargo, no se encontró especies categorizadas 
como Intermedias y Valiosas en “Paolita |”, asimismo, en “E&O-2” no se encontró especies 
categorizadas como Intermedias. Por otro lado, las especies categorizadas como “otras 
especies” con un valor económico futuros fueron encontrada en la mayor proporción en 
las concesiones “Mary Elena Erica” y “E&0-2” (Tabla 10). Diferentes resultados fueron 
encontrados por (AREVALO, 2015; PEREZ, 2013). 








Categoría Denominación mo E&0-2 Paolita 1 Paolita 2 
A Altamente valiosas - - - - 
B Valiosas 5% 10% - 15% 
Cc Intermedias 9% - - 38% 
D Potenciales 14% 20% 43% 38% 
E Otras especies (Valor 73% 70% 57% 62% 


económico futuro) 





Tabla 10. Proporción de especies según categoría en las 4 concesiones evaluadas. 


CONCLUSIONES 


Los valores de la madera al estado natural en las concesiones mineras estudiadas 
oscilaron entre S/216/ha y S/ 1305/ha. Las diferencias en los valores al estado natural de la 
madera entre las concesiones estudiadas están influenciadas por la presencia de algunas 
especies con un elevado valor económico (categorías By C). La riqueza de especies fue 
inferior al reportado por otros estudios en Madre de Dios. 
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RESUMO: A ocorrência e dinâmica populacional 
de insetos nos ecossistemas, tende a se manter 
em níveis estáveis, com pequenas flutuações 
constantes. Fatores endógenos (e.g., migração/ 
imigração, estresse, interações sociais, pool 
genético) e/ou exógenos (e.g., alimento, espaço, 
inimigos naturais, clima), no entanto, podem 
proporcionar o aumento pontual e explosivo 
na abundância de insetos e levar a surtos 
populacionais. Todavia, um surto populacional 
nem sempre representa um problema ecológico. 
A variação temporal e espacial em um recurso ou 
fator pode induzir um aumento exponencial de 
uma população, como resposta direta ao aumento 
na oferta. Devemos, portanto, compreender 
a concepção de diferentes ambientes em sua 
complexidade ecológica. Em áreas conservadas, 
a resiliência do ecossistema será capaz de 
sobrepor o potencial de dano de um surto 
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e retornar ao estado de equilíbrio, ou gerar 
um novo. Já áreas degradadas condicionam 
maior vulnerabilidade à ocorrência de surtos, 
com modificações significativas na riqueza e 
estrutura das comunidades de insetos locais. 
Os monocultivos são ainda mais favoráveis à 
existência de episódios de surtos, em função da 
concentração de uma única espécie vegetal e do 
manejo intensificado. A necessidade de controle 
dos surtos será dependente da composição, 
estrutura e do potencial de resposta ecológica 
dos ecossistemas naturais e plantados, e 
dos possíveis níveis de danos econômicos e 
ecológicos, considerando sempre a manutenção 
do equilíbrio dinâmico e não a erradicação do 
inseto-praga. 

PALAVRAS - CHAVE: espécie-praga, surto 
populacional, complexidade ecológica, análise 
sistêmica. 


AN ANALYSIS ON POPULATION 
DYNAMICS AND PEST INSECTS 
OUTBREAK 


ABSTRACT: The occurrence and population 
dynamics of insects in ecosystems tend to remain 
at stable levels, with small constant fluctuations. 
Endogenous (e.g., migration/ | immigration, 
stress, social interactions, gene pool) and/or 
exogenous (e.g., food, space, natural enemies, 
climate) factors, however, can provide a punctual 
and explosive increase in the abundance of 
insects and lead population outbreaks. However, 
population outbreaks not always represent an 
ecological problem. The temporal and spatial 
variation in a resource or factor can induce an 
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exponential increase in a population, as a direct response to the increase in supply. We must, 
therefore, understand the design of different environments in their ecological complexity. In 
conserved areas, the ecosystem's resilience will be able to overcome the damage potential 
of an outbreak and return to equilibrium or generate a new one. Degraded areas, on the 
other hand, conditioning greater vulnerability to the occurrence of outbreaks, with significant 
changes in the richness and structure of local insect communities. Monocultures are even 
more favorable to the existence of episodes of outbreaks, due to the concentration of a 
single plant species and intensified management. The need to control outbreaks will depend 
on the composition, structure, and potential for ecological response of natural and planted 
ecosystems, and the possible levels of economic and ecological damage, always considering 
the maintenance of dynamic balance and not the eradication of the insect pest. 

KEYWORDS: pest species, populational outbreaks, ecological complexity, systemic analysis. 


11 INTRODUÇÃO 


Quando se fala na relação inseto-planta, uma das primeiras concepções que 
emergem é a existência de “pragas”, muito em função dos cultivos produtivos humanos, 
que abordam a ocorrência de determinados insetos como um problema a ser resolvido. 
Coevolutivamente, o estudo e manejo de pragas construiu-se como um subconjunto 
da história da agricultura e, mais recentemente, da silvicultura (DENT, 2000). De forma 
que muitos dos atuais problemas no controle de flutuações explosivas nas dinâmicas 
populacionais de insetos, são oriundos de ações falhas do passado, tanto para proteção 
dos cultivos, quanto para aumento da produtividade (WAAGE, 1993). 

Aintensificação da agricultura e da silvicultura trouxeram novos e maiores problemas 
relacionados aos surtos de insetos-praga (DENT, 2000). Alterações na dinâmica, estrutura 
e composição das populações naturais de insetos, no entanto, sempre existiram. E 
talvez esteja na análise e compreensão destes mecanismos naturais, a chave para a 
solução de problemas futuros na produtividade e conservação dos sistemas naturais e 
agroecossistemas. 

Dada a importância do assunto, apresentamos uma breve revisão, elucidando 
fatores fundamentais sobre a dinâmica populacional e os surtos em insetos-praga: |) como 
é, em geral, a dinâmica populacional em insetos; Il) o que são surtos e quais as suas 
causas; Ill) qual o impacto dos surtos nos diferentes ambientes e IV) quando é necessário 
intervenção. 


21 METODOLOGIA 


O presente trabalho consiste numa revisão sobre dinâmica populacional e surtos de 
insetos-praga, com base no livro Entomologia Florestal, de autoria de Costa et al. (2014), 
e na literatura disponível na Plataforma Google Scholar, acessada em dezembro de 2018, 
tendo como termos de busca as palavras: “dinâmica populacional”, “surto” e “insetos-praga”. 
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31 DESENVOLVIMENTO 


3.1 Dinâmica populacional em insetos 


Todos os fatores que afetam a densidade populacional, a composição e a estrutura 
dos ecossistemas atuam de maneira a regularem comunidades e populações dentro de 
limites estreitos de densidade e equilibro de espécies (TRIPLEHORN; JONNSON, 2011). 
A tendência nos ecossistemas é que as populações se mantenham em níveis estáveis 
(RIBEIRO; PAGOTTO, 2011), em busca de um equilíbrio, mas com pequenas flutuações 
constantes. Esse equilíbrio, portanto, não é invariável, e nem permanece em um grau 
de entropia crescente (CAPRA; LUISI, 2014). Da mesma forma constitui-se a dinâmica 
populacional em insetos. O número de indivíduos é determinado, em primeiro lugar, pelo 
tamanho da comunidade, enquanto sua densidade é controlada por fatores endógenos e 
exógenos (TRIPLEHORN; JONNSON, 2011). 

Como fatores exógenos destacam-se o alimento, o espaço, a existência de inimigos 
naturais (fatores dependentes da densidade) e o clima (fator independente da densidade) 
(TRIPLEHORN; JONNSON, 2011). Já os fatores endógenos são, em sua totalidade, 
dependentes da densidade, e representam respostas intrínsecas das espécies e/ou 
indivíduos aos diferentes níveis de interação (migração/ imigração, aumento do estresse, 
interações sociais, pool genético, etc.) (TRIPLEHORN; JONNSON, 2011). Cabe ressaltar 
que, na prática, os fatores atuam de maneira conjunta, em um rearranjo complexo de 
caracteres, que impedem o isolamento de um fator como causa das alterações na dinâmica 
populacional, e dificultam a previsão e interpretação dos seus efeitos. É, contudo, em 
resposta a estes fatores que a seleção natural atuará como elemento determinante da 
viabilidade das populações, exigindo a manutenção dos nichos ecológicos para a expressão 
do potencial biótico das espécies e para sua manutenção em níveis relativamente estáveis 
(RIBEIRO; PAGOTTO, 2011). 


3.2 O que são surtos populacionais em insetos? O que causa um surto? 


Para Costa et al. (2014), os surtos populacionais correspondem a um aumento 
inusitado, pontual e explosivo na abundância das populações de insetos. Podem ser 
classificados em: eruptivos — quando apresentam alta densidade populacional, com 
manutenção espacial e temporal; gradientes — quando dependem de eventos extremos, 
tanto endógenos quanto exógenos; e comportamento temporal — quando ocorrem 





periodicamente, sendo dependentes da dinâmica populacional de inimigos e hospedeiros 
(COSTA et al., 2014). Essa definição de surto, vinculada diretamente ao espaço-tempo, 
é especialmente importante para evitar caracterizações ou diagnósticos errôneos. O 
aumento populacional gradual, ao longo de um prolongado espaço de tempo, apesar de 
comum em grupos e espécies de insetos, não significa o surgimento de um surto, mas sim 


um rearranjo, readaptação comportamental, uma resposta fisio-genética às alterações de 
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nicho, ou ainda, o processo de equilíbrio dinâmico em ação (TRIPLEHORN; JONNSON, 
2011; CAPRA, 2014; BARKER et al., 2017). 

Em alguns modelos, os surtos só são possíveis na existência de um antagonista 
especializado (modelo hospedeiro/ patógeno) (DWYER etal., 2004). Ou só ocorrem quando 
os predadores generalistas falham (modelo generalista/ predador) (DWYER et al., 2004). 
Nenhum destes, no entanto, é suficiente para, por si só, caracterizar o comportamento 
episódico de surto (DWYER et al., 2004). Para compreender como é possível a existência 
destes surtos, é fundamental elucidar aspectos eco/fisiológicos da vida de um inseto. Em 
geral, os insetos possuem uma admirável capacidade reprodutiva, competitiva e eficiente. 
O número de indivíduos viáveis produzidos a cada geração é dependente, essencialmente, 
de três fatores: do número de ovos férteis depositados pela fêmea, da duração de uma 
geração (tempo de vida) e da proporção de fêmeas (para produção da próxima geração) 
(TRIPLEHORN; JONNSON, 2011). 

A alta capacidade reprodutiva, somada a fatores comportamentais, fisiológicos 
e ambientais favoráveis, permitem que, dentro da zona limiar, populações aumentem 
exponencialmente de tamanho e possam vir a causar danos econômicos e ecológicos 
(DENT, 2000). O conhecimento dos fatores que influenciam a dinâmica populacional em 
insetos é fundamental para prever surtos, com destaque para fatores críticos como o 
ambiente e os critérios de seleção do hospedeiro (MOUSSALEM; SANTOS-SILVA, 2007). 
As inúmeras modificações que um inseto pode sofrer ao longo de um ciclo de vida (em 





aspecto: ovo, pupa, larva, adulto; em hábitos alimentares: consumo de folhas, caules, 





raízes, flores, frutas; em comportamento: com resposta a estímulos ambientais, estado 
interno, tipo de experiência) (TRIPLEHORN; JONNSON, 2011), também contribuem para a 
possibilidade de surtos. E a condensação da fase predadora no tempo e espaço facilita o 
surgimento de surtos associados ao momento de vida. 

Os surtos podem ser causados, portanto, em decorrência de transformações 
no ambiente físico e biótico, como alterações no uso do solo, eventos meteorológicos 
extremos, degradação dos ecossistemas, alterações e mutações genéticas (naturais ou 
não), interações tróficas mutualistas ou antagonistas, mudanças quali e quantitativas na 
estrutura, composição e dinâmica do ecossistema, redução drástica da população de 
algum agente da cadeia trófica, entre outras (COSTA et al., 2014). 


3.3 Qual o impacto ecológico dos surtos? Até que ponto um surto de insetos 
é um problema ecológico? 

A destruição dos nichos ecológicos decorre na quebra de cadeias tróficas, 
principalmente em decorrência da perda de espécies e consequente simplificação 
dos ecossistemas. Apesar disso, surtos de espécies-praga não são pressupostos da 
simplificação dos ecossistemas, apesar de mais prováveis neles (NAIR, 2007; WHYLE; 
SPEIGHT, 2012). A monocultura, ou dominância ecológica, também é frequente em 
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ecossistemas naturais, onde os mecanismos operam dentro do equilíbrio esperado. E 
surtos ecológicos também são descritos em áreas teoricamente complexas, biodiversas e 
em bom estado de conservação (DENT, 2000). 

Borges et al. (2003), analisando flutuações populacionais de insetos em diferentes 
sistemas de cultivo de lIlex paraguariensis A.St.-Hil., concluíram que as flutuações 
associadas aos insetos com potencial predador estão presentes mesmo no sistema de 
cultivo nativo, em áreas em bom estado de conservação. Barbosa et al. (2012), Bjorkam 
e Niemela (2015), também demonstraram que, em sistemas florestais, surtos cíclicos de 
pragas ocorreriam associados às mudanças naturais na condição fisiológica do hospedeiro, 
tendo o clima como fator preponderante. 

Dessa forma, um surto populacional nem sempre representa um problema ecológico. 
A própria variação temporal e espacial em um determinado recurso ou de um fator pode 
ser indutora de um aumento exponencial de uma população, como reposta interativa direta 
ao aumento na oferta. Entender como essas interações são capazes de gerar mudanças 
significativas na abundância, composição ou estrutura das espécies tem sido, justamente, 
um dos principais objetivos da bioecologia nos últimos séculos (BRONSTEIN et al., 2003). 

Calore et al. (2013), ao analisarem pomares de Psidium guajava L., observaram que 
existe uma correlação climática com os aumentos populacionais, mas que os principais 
picos coincidem com a frutificação. E que o pico populacional do principal predador 
não coincide com o da presa, indicando um provável mecanismo de fuga e dinâmicas 
ecológicas de coexistência, mesmo em áreas plantadas. Borges et al. (2003), apesar de 
terem constatado aumento populacional no período de disponibilidade de recurso e/ou 
período reprodutivo dos insetos considerados praga da cultura de llex paraguariensis, tanto 
em áreas nativas quanto em plantios adensados, também não observaram grandes danos 
nas áreas de cultivo nativo, onde havia maior diversidade de espécies associadas. 

O problema ecológico surge, portanto, quando as alterações na estrutura, composição 
e dinâmica, ultrapassam os limites ecológicos de resiliência do ecossistema, fazendo com 
que as populações se desenvolvam fora de controle, para além da capacidade de equilíbrio 
de outras espécies, promovendo através dos surtos populacionais, perdas ecológicas e 
econômicas. 


3.4 O impacto de um surto é igual em diferentes ambientes? 


Primeiramente devemos compreender a concepção de diferentes ambientes em 
sua complexidade ecológica, não em limites políticos e barreiras fisiográficas. O ambiente 
consiste de um intrínseco e complexo emaranhado de relações, interações e iterações, 
não estocásticas, em sua maioria não previsíveis e irreprodutíveis (THOMPSON, 2005; 
MARTINS, 2012). Logo, diferentes estruturas, processos e composições, geram respostas 
ecológicas distintas (THOMPSON, 2005). 

Em áreas conservadas, estruturadas e equilibradas, mesmo que alguma alteração 
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pontual seja capaz de provocar um surto demográfico em uma espécie, a resiliência do 
ecossistema, sustentada na existência de diversidade, qualidade ambiental, múltiplas 
interações e dinâmicas complexas (THOMPSON, 2005; MARTINS, 2012), será capaz 
de sobrepor o potencial de dano e retornar ao estado de equilíbrio, ou gerar um novo 
estado de equilíbrio (CAPRA; LUISI, 2014; THOMPSON, 2005). Já em áreas degradadas, a 
alteração danosa do ponto de vista da composição, estrutura e riqueza de espécies basais, 
condicionam maior vulnerabilidade a invasão de espécies exóticas e a surto de insetos, 
com modificações significativas na riqueza e estrutura das comunidades de insetos locais 
(KNOPS et al., 2002). 

Lopes et al. (2005), constataram diferença na dinâmica e estrutura de populações 
de besouros predadores (Histeridae) ao compararem quatro ambientes com graus 
decrescentes de complexidade estrutural da vegetação, com tendência de decréscimo na 
diversidade e equidade, e substituição de espécies com a diminuição da complexidade do 
ambiente. Já em 1987, Powell e Powell observaram uma redução no número de indivíduos 
visitantes e na composição de espécies de abelhas Euglossinae em áreas fragmentadas. 
Becker et al. (1991), estudando os mesmos fragmentos após 5 anos, observaram um 
aumento da abundância de abelhas. Em 1998, esse fenômeno controverso observado 
por Becker et al. (1991) foi denominado comunidade saturada. Revelando uma importante 
dinâmica ecossistêmica, o aumento pontual da diversidade em resposta à perturbação e 
posterior queda brusca na diversidade e em outros descritores ecológicos. Thomazini e 
Thomazi (2000), descrevem alterações negativas visíveis e quantificáveis nas populações 
de insetos em áreas degradadas, se comparadas as áreas conservadas. 

Os monocultivos, neste cenário, também são responsáveis por adicionar fatores 
favoráveis a existência de episódios de surtos: a concentração de uma única espécie/ 
variedade, abundância alimentar/ nidificação nos cultivares de alta produtividade, redução 
de inimigos naturais, menores intervalos entre disponibilização de recursos, aumento da 
exposição e, conseguinte, da resistência a agrotóxicos, e introdução de espécies exóticas 
(DENT, 2000). Nestes ambientes agrícolas ou silviculturais, apesar da existência de uma 
estrutura formada por espécies de base, ocorre drástica simplificação do ecossistema, 
diminuindo sua capacidade de resistir ou até mesmo de se recuperar de prováveis danos 
(LINDENMAYER et al., 2016; SALES et al., 2018). 


3.5 A partir de que ponto é necessário a intervenção? 


A resposta ecológica, tanto relacionada a resistência quanto a resiliência à surtos 
populacionais, será diretamente dependente da composição e estrutura dos ecossistemas. 
Logo, a necessidade de monitoramento, intervenção e controle, também. Qualquer ação 
sobre o ecossistema demanda um conhecimento profundo dos mecanismos de controle 
natural de populações e da ecologia sistêmica, buscando prever surtos com confiabilidade, 
e controlando-os de maneira eficiente (TRIPLEHORN; JONNSON, 2011). Sem, no entanto, 
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desencadear outros desequilíbrios ecológicos (RIBEIRO; PAGOTTO, 2011). 

Apesar de diversificada literatura sobre o Manejo Integrado de Pragas, que deveria 
orientar o processo de análise, prevenção e controle de surtos populacionais. Na prática, 
as ações realizadas são, em sua maioria, paliativas (BARROS et al., 2006). A utilização 
indiscriminada de agrotóxicos de maneira profilática vem agravando os problemas 
ambientais e promovendo novos desequilíbrios (SILVA; COSTA, 2002; BARROS et al., 
2006). Isto porque, apesar de existirem alternativas, o impacto e a influência das grandes 
indústrias ainda obscurecem o poder da pesquisa de fontes alternativas. 

É fundamental elucidar que não existem insetos-praga se separados do contexto, 
pois, por definição, pragas são populações de organismos capazes de causar injúrias ou 
modificações estruturais apenas quando a ponto de causar danos econômicos (PEDIGO; 
RICE, 2006) ou ecológicos. Além de que, todas as espécies existentes são fundamentais 
para a manutenção dos processos e serviços ecossistêmicos. O objetivo das estratégias 
de intervenção jamais deve perpassar a possibilidade de erradicação da espécie no local 
(BARBOSA et al., 2012), mas apenas a manutenção dos níveis populacionais dentro 
dos limites do Nível de Dano Econômico (NDE) e do Nível de Controle (NC) (SILVA; 
COSTA, 2002). A intervenção será necessária, portanto, apenas quando estes limiares 
forem atingidos, principalmente tratando-se de áreas produtivas, onde o objetivo não é a 
complexificação do ecossistema. Por isso a importância do monitoramento como principal 
ferramenta de controle (SHONE, 2002). 

Uma importante ferramenta de prevenção e monitoramento é a modelagem 
matemática, que permite análises mais seguras dos fenômenos e enriquece o processo 
de controle (RICKLEFS; RELYEA, 2016; SHONE, 2002), principalmente por permitir 
diagnósticos personalizados, contextualizados e situacionais, independente das expectativas 
ou concepções pré-estabelecidas. Considerando-se áreas produtivas, diversas ferramentas 
e modelos matemáticos para cálculo do NDE e NC já foram desenvolvidos, além de 
programas e estratégias de intervenção e monitoramento. Os mesmos pressupostos nem 
sempre poderão ser aplicáveis em áreas naturais (sejam elas degradadas, conservadas 
ou em processo de restauração), pois nestes locais, o objetivo primordial não é a máxima 
produtividade com o menor dano possível, e sim, o retorno dos processos ecossistêmicos 
em um estado de equilíbrio dinâmico. Para estas áreas, uma ferramenta que desponta é 
a análise sistêmica, através da modelagem matemática de redes de interação (CAPRA; 
LUISI, 2014; CRESTANI et al., 2019; BOVO et al., 2018; GENRICH et al., 2017), que pode 
possibilitar diagnósticos precisos, baseados em fatores de estruturação do ambiente, 
apontando caminhos para o processo de restauração e conservação dos recursos naturais. 


41 CONSIDERAÇÕES FINAIS 


Nas populações naturais, a estrutura das comunidades de insetos é variável de 
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acordo com as espécies e ambientes interantes, mas se mantêm estável, com pequenas 
variações constantes, num estado denominado de equilíbrio dinâmico. Quando algum fator, 
endógeno ou exógeno, como o aumento na oferta de recursos ou a eliminação de um 
predador natural, interfere nesse equilíbrio, pode ocorrer um aumento pontual e explosivo 
na abundância das populações de insetos, decorrendo em surtos. Como a manutenção 
do equilíbrio nas populações depende do estado de equilíbrio do ambiente, as áreas 
com monocultivos, em decorrência da concentração de uma única espécie vegetal e da 
intensidade de manejo, são mais favoráveis à existência de surtos episódicos do que os 
ambientes naturais diversos e conservados. 

A necessidade e a intensidade de intervenção sobre episódios de surto de insetos 
são dependentes da capacidade de resistência/resiliência do ecossistema e, em áreas 
de produção comercial, do nível de dano econômico. Contudo, qualquer ação sobre os 
ecossistemas, naturais ou plantados, requer amplo conhecimento dos mecanismos de 
controle natural de populações e da biologia e ecologia dos insetos, além de uma análise 
criteriosa de custos e benefícios, para o setor produtivo e para os ecossistemas. Neste 
sentido, o que seguramente compreendemos sobre dinâmica populacional de insetos 
e surtos de insetos-praga é ainda a ponta do iceberg, em relação à complexidade dos 
ecossistemas. 
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RESUMO: O setor florestal contribui com a 
geração de emprego e renda no Brasil, sendo 
um setor bastante desenvolvido, contudo, 
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principalmente em áreas declivosas, uma serie 
de melhorias ainda devem ser implementadas 
e analisadas com relação às condições de 
saúde e segurança dos trabalhadores. Dentre 
as principais operações florestais em áreas 
declivosas, está a colheita florestal, uma das 
mais perigosas e com alto risco de acidentes. 
Diante deste contexto, esta revisão bibliográfica 
teve por objetivo avaliar os fatores de riscos que 
interferem nas atividades de colheita florestal em 
áreas declivosas, a fim de encontrar correções 
ou modificações necessárias que auxiliem na 
melhoria da qualidade de vida e trabalho para 
os trabalhadores. Para o desenvolvimento do 
trabalho foi feita uma análise bibliográfica, sendo 
constatado que grande parte dos acidentes 
ocorre no momento do corte e derrubada de 
árvores. Os resultados demonstraram que a 
colheita florestal oferece diferentes riscos, com 
destaque para os riscos físicos, ergonômicos e 
de acidentes. Para a revisão, foi realizada uma 
análise desses principais riscos, com base nas 
diretrizes das Normas Regulamentadoras Nº 
6, 17 e 31. Denota-se a falta de equipamentos 
adequados e a utilização de equipamentos e 
máquinas inseguros destinados a esse setor 
da atividade, entre eles a motosserra. Também 
foram verificados que muitos acidentes de 
trabalho acontecem pela não qualificação, falta 
de treinamentos e falta de profissionalização dos 
trabalhadores. Por fim, a condição de declividade 
do relevo montanhoso pode aumentar a 
quantidade e a gravidade dos riscos envolvidos 
nesta atividade, o que demonstra a necessidade 
de outras análises de riscos da colheita florestal 
para as mais diversas condições. 
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PALAVRAS - CHAVE: Norma Regulamentadora; Setor Florestal; Declividade. 


RISK ANALYSIS ASSOCIATED WITH FOREST HARVESTING IN DECLIVING 
AREAS IN BRAZIL 

ABSTRACT: The forestry sector contributes to the generation of jobs and income in Brazil, 
being a well-developed sector, mainly in the declared areas, a series of improvements still 
have to be implemented and analyzed in relation to the health and safety conditions of 
workers. Among the main forestry operations in sloping areas, is the forest harvest, one of the 
most dangerous and with a high risk of accidents. Given this context, this bibliographic review 
aimed to assess the risk factors that interfere with forest harvesting activities in sloping areas, 
to find necessary corrections or modifications that help to improve the quality of life and work 
for workers. For the development of the work, a bibliographic analysis was made, and it was 
found that most accidents occur when trees are cut and felled. The results showed that the 
forest harvest offers different risks, with emphasis on the physical, ergonomic and accident 
risks. For the review, an analysis of these main risks was carried out, based on the guidelines 
of Regulatory Norms No. 6, 17 and 31. There is a lack of adequate equipment and the use 
of unsafe equipment and machinery for this sector of activity, among others the chainsaw. It 
was also verified that many work accidents happen due to non-qualification, lack of training 
and professionalization of workers. Finally, the sloping condition of the mountainous relief can 
increase the quantity and severity of the risks involved in this activity, which demonstrates 
the need for further analysis of the risks of forest harvesting for the most diverse conditions. 
KEYWORDS: Regulatory Standard; Forestry Sector; Declivity. 


INTRODUÇÃO 


O setor florestal tem crescido na última década, e sua importância econômica se 
apresenta em números expressivos, com representatividade de 6,2% do PIB Industrial no 
país (IBA, 2020). No entanto, as florestas plantadas no Brasil têm avançado em áreas 
que não são consideradas agricultáveis, como as áreas declivosas, o que representa um 
grande desafio para a mecanização (LEITE et al., 2011) e todas as atividades ligadas ao 
setor. 

O setor florestal é considerado um dos mais perigosos mundialmente, despertando 
a atenção de pesquisadores e autoridades para a criação de medidas de segurança e 
proteção aos trabalhadores (ASSUNÇÃO; CAMARA, 2011). Segundo Souza et al. (2010), 
anualmente, ocorrem inúmeros acidentes no setor florestal, em que, grande parte deles 
ocorre na fase de colheita, nas operações de derrubada, desgalhamento, destopamento e 
toragem. A colheita florestal é vista como uma atividade de alto risco de acidentes, devido 
à grande exigência física requerida ao trabalhador, principalmente quando realizada de 
forma manual ou semimecanizada (ASSUNÇÃO; CAMARA, 2011). 

Na colheita florestal, os riscos de acidentes mais comuns são aqueles que incluem 


as máquinas e equipamentos, principalmente na fase de tombamento da árvore, e também 
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os riscos de picadas de animais peçonhentos (CANTO et al., 2007). Os riscos ergonômicos 
estão associados ao transporte manual de peso e aos ritmos excessivos, com jornadas de 
trabalho prolongadas em turno noturno, monotonia e repetitividade. Além desses riscos, há 
o alto nível de ruído e vibração das máquinas (BARBOSA et al., 2014). 

Em consonância com estes riscos, evidencia-se que o trabalho, em grande maioria, 
é realizado ao ar livre, ou seja, o trabalhador fica exposto às intempéries do clima e suas 
consequências, sofrendo com o calor ou frio, com a umidade, os ventos, entre outros. 
Entre outras causas de acidentes nesse setor, pode-se citar a falta de experiência e 
treinamentos de trabalhadores, tanto para a realização das atividades, quanto para a 
utilização das máquinas, ferramentas e até mesmo equipamentos de segurança, bem como 
o desconhecimento dos perigos que essa atividade oferece (VIEIRA, 2013). De acordo 
com Assunção e Câmara (2011), a adoção de medidas preventivas visando a resguardar a 
saúde e a segurança dos trabalhadores, sem dúvida nenhuma, revertem-se em benefícios 
sociais e econômicos para esses trabalhadores e para a sociedade em geral. 

Ao relatar os riscos eminentes aos trabalhadores, pode-se notar que a colheita 
florestal requer medidas preventivas extremamente rigorosas, visando a segurança e 
saúde ocupacional. No entanto, estudos científicos acerca da atividade de colheita florestal 
em condições de áreas declivosas são ainda incipientes no Brasil, principalmente com 
foco na segurança do trabalho desta atividade, em busca da prevenção de acidentes e 
condições de trabalho ideais. 

Diante deste cenário, objetivou-se realizar um levantamento e análise dos principais 
riscos de acidentes e as possíveis melhorias na etapa de colheita florestal, em áreas 
declivosas e com emprego de trabalho manual e ou semimecanizado. 


MATERIAL E MÉTODOS 


O presente trabalho se caracteriza como uma revisão de literatura com o intuito de 
reunir o que há de mais atual em termos de conhecimento científico relacionado ao tema 
de riscos de acidentes associados à atividade de colheita florestal em áreas declivosas. 

A pesquisa foi realizada por meio de observações e considerações das Normas 
Regulamentadoras NR-6, NR-17 e NR-31 e de pesquisas publicadas em revistas cientificas 
de relevante impacto e anexadas nas bases de dados Google Scholar, Scopus, Periódicos 
Capes e Isi Web of Knowledge limitadas no período de 2000 a 2021. 


DESENVOLVIMENTO 


Colheita florestal 


No setor florestal, a colheita é considerada a fase mais importante do ponto de 
vista econômico, uma vez que esta representa em torno de 50% ou mais do custo de 
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produção da madeira, com uma alta participação no custo final do produto (MACHADO, 
2014; CONEGLIAN et al., 2010). Além disso, caracteriza-se como uma fase do processo 
produtivo de fundamental importância para o potencial de crescimento do setor florestal, 
visto que pode influenciar significativamente, não só nas questões econômicas, mas como 
também nas questões ambientais e sociais (CARVALHO; SOARES; VALVERDE, 2005). 

A colheita florestal pode ser definida como um conjunto de operações efetuadas no 
povoamento florestal, que compreende, desde o preparo, até a destinação da madeira ao 
local de transporte ou utilização, fazendo-se o uso de técnicas e padrões estabelecidos, 
com a finalidade de transformá-la em produto final (MACHADO, 2014). O processo 
produtivo de colheita da madeira, em tese, é composto pelas etapas de corte (derrubada, 
desgalhamento e processamento ou traçamento), descascamento, quando executado em 
campo, extração e carregamento (SILVA, 2003). 

Atualmente, há uma grande diversificação dos sistemas de colheita de madeira, cuja 
“implementação” varia de acordo com a topografia do terreno, com o uso final da madeira, 
com as máquinas, equipamentos e recursos tecnológicos disponíveis (MINETTI etal., 2008). 
A modernização das operações florestais no Brasil, principalmente aquelas que fazem parte 
do processo de colheita e transporte de madeira, intensificou-se significativamente no final 
da década de 1990, com a abertura do mercado pelo governo brasileiro à importação de 
máquinas e equipamentos de países com tradição na mecanização florestal (MACHADO, 
2014). 

Somado a isso, o aumento da produtividade das florestas e dos custos de mão de 
obra e a necessidade de executar o trabalho com maior segurança operacional e de redução 
nos custos de produção, culminaram numa intensificação dessa modernização (MINETTE 
et al., 2008). Essas circunstâncias levaram empresas brasileiras a passarem da colheita 
manual ou semimecanizada, para sistemas totalmente mecanizados, com máquinas de alta 
tecnologia, produtividade e elevados custos (LEONELLO; GONÇALVES; FENNER, 2012). 

Esta necessidade de mecanização da colheita florestal acarretou na evolução 
significativa dos sistemas e máquinas empregados nesta atividade, visando diminuir os 
custos da matéria prima e aumentar a produtividade (MINETTE et al., 2008). Entretanto, 
a demanda por sistemas de colheita de madeira mais adequados quanto à segurança do 
trabalho ainda é sentida pelo setor. 


Determinação e classificação dos riscos ocupacionais 


Define-se “risco” como o grau de probabilidade de ocorrência de um determinado 
evento e o “fator de risco” está associado ao aumento de probabilidade de ocorrência 
do agravo à saúde, sem que o referido fator tenha que interferir, necessariamente, em 
sua causalidade (BRASIL, 2021b). Os riscos estão presentes nos mais diversos locais de 
trabalho e em todas as demais atividades humanas, comprometendo a segurança, a saúde 
das pessoas e a produtividade de empresas (BRASIL, 2021b). 


Conceitos e Conhecimentos de Métodos e Técnicas de Pesquisa Científica em Engenharia à 
Florestal 2 Capítulo 8 


Os fatores de risco podem ser sinais ou causa de eventos indesejados, mas em 


qualquer circunstância e indistintamente, devem ser identificados ou observados antes 


da ocorrência do evento prenunciado (BRASIL, 2021b). A Norma Regulamentadora Nº 09 


considera riscos ambientais, como os agentes físicos, químicos, biológicos, ergonômicos 


e, de acidentes, com potencial de afetar diretamente a saúde do trabalhador, expondo-o a 
doenças e acidentes de trabalho (BRASIL, 2021b). 
Ao observar o Quadro 1, pode-se compreender por meio da exemplificação da 


classificação e agrupamento dos principais riscos ocupacionais, considerando-se a 


natureza e a padronização das cores correspondentes. 
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Quadro 1 — Classificação dos riscos, de acordo com a sua natureza e a padronização das cores 


correspondentes - NR 9. 


Table 1 - Classification of risks, according to their nature and the standardization of the corresponding 
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colors - NR 9. 
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Riscos associados à colheita florestal no Brasil 


O setor florestal brasileiro vem ao longo dos anos sendo considerado um dos mais 
perigosos nos aspectos de segurança do trabalho do mundo (SOUZA et al., 2010). Os 
maiores riscos, segundo alguns autores, estão presentes nas operações de corte das 
árvores, manejo dos equipamentos e máquinas para o transporte dos materiais (SILVA, 
2013). 

Diante deste contexto e ao observar o Quadro 2, pode-se notar que a maioria dos 
artigos científicos encontrados apontam que os principais agentes de riscos presentes no 
ambiente de trabalho florestal estão relacionados à colheita florestal, principalmente em 
áreas declivosas. Estes ocorrem fundamentalmente em razão da crescente exposição do 
trabalhador a situações de tensão e estresse no ambiente de trabalho e também ao perigo 


pelo manuseio de máquinas, equipamentos e ferramentas, muitas vezes sem o treinamento 
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Quadro 2 — Principais estudos sobre riscos de acidentes de trabalho no Setor Florestal Brasileiro. 


Table 2 — Main studies on occupational accident risks in the Brazilian Forest Sector. 


No Brasil, a colheita florestal é, em sua grande maioria, realizada de forma manual 
ou semimecanizada, a qual é variável em função das condições de topografia dos terrenos 
e de fatores econômicos, ambientais e sociais, o que de fato a caracteriza como uma 
atividade pesada e perigosa (COELHO; LIMA; FONSECA, 2017). Essas atividades exigem 
dos operadores o manuseio de cargas elevadas e máquinas e ferramentas perigosas, 
com esforços repetitivos e posturas inadequadas durante praticamente toda a jornada de 
trabalho (VOSNIAK et al., 2010). 

Além disso, as condições e o ambiente de trabalho têm aspectos particulares, pois 
os mesmos são temporários e os trabalhadores atuam expostos a condições climáticas 
adversas, o que aumenta o risco de acidentes (MINETTI et al., 2008). Trabalhando ao ar 
livre, o empregado fica exposto às intempéries do clima e suas consequências, sofrendo 
com o calor ou frio, com a umidade, os ventos, entre outros (MEDEIROS; JURADO, 2013). 

Muitas vezes, o local de trabalho fica distante de sua residência, obrigando o 
trabalhador a dispender tempo e energia no trajeto, correndo o risco de sofrer acidentes 
(MEDEIROS; JURADO, 2013). Devido ao isolamento do local de trabalho, geralmente 
faltam facilidades para o atendimento médico e de primeiros socorros com a urgência 
necessária. 

Canto et al. (2007) evidenciaram a necessidade de conscientização das empresas e 
dos trabalhadores sobre o risco de acidentes inerentes à colheita florestal e a importância 
da utilização de equipamentos de proteção individual, bem como de instrução sobre como 
o trabalhador deve realizar a sua atividade de forma segura. Segundo Jacovine et al. 
(2005), a melhoria da qualidade do processo pode ser alcançada com investimento em 
treinamentos dos trabalhadores florestais. 

Apesar do avanço nos últimos anos, a melhoria das condições de trabalho ainda 
tem sido muito modesta no setor florestal, uma vez que as pesquisas estão mais voltadas 
para os aspectos de otimização do trabalho, redução de custos e aumento da produtividade 
(SOUZA et al., 2010). A falta de experiência nessas operações e a falta de treinamento 
pessoal são as principais causas de acidentes no trabalho, em especial, destacam-se os 
acidentes ocorridos com os operadores de motosserra em áreas declivosas (CANTO et 
al., 2007). Ademais, para evitar os acidentes no processo da colheita florestal, algumas 
técnicas adequadas devem ser adotadas, como também medidas preventivas, além de 
serem aplicadas as normas e leis vigentes (PESCADOR et al., 2013). Esses fatores 
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integrados irão determinar o melhor desempenho das atividades, bem como a melhor 
utilização dos recursos disponíveis de forma segura (CANTO et al., 2007). 


Riscos físicos associados à colheita florestal 


Na colheita florestal, o ruído é apontado como um dos principais problemas para 
o operador de máquinas, principalmente em áreas declivosas com o uso de motosserra 
(CUNHA, 2001). Os problemas auditivos causados pelo ruído são determinados pelo nível 
de pressão sonora, frequência e tempo de exposição (BRASIL, 2021a). 

A exposição intensa e prolongada ao ruído atua desfavoravelmente sobre o estado 
emocional do operador, com consequências imprevisíveis sobre o equilíbrio psicossomático, 
além de diminuir o desempenho do trabalhador nas suas funções e aumentar a probabilidade 
de ocorrência de acidentes no trabalho (KROEMER; GRANDJEAN, 2005). 

Diversos estudos realizados no setor florestal e madeireiro apontaram níveis 
elevados de exposição ao ruído com potencial de causar danos à saúde dos trabalhadores 
(SOUZA et al., 2010; FIEDLER et al., 2010). Minette et al. (2007) ao avaliarem os níveis 
de ruído em diversas máquinas de colheita florestal, constataram que todas estavam em 
desconformidade com a legislação brasileira. 

Assim, nota-se que a Associação Brasileira de Normas Técnicas (ABNT) estabelece 
os níveis de ruído para conforto acústico, sendo o nível de ruído de 85 dB permitido por 8 
horas de exposição (BRASIL, 2021a). Para cada aumento de 5 dB no nível de ruído acima 
deste limite, o tempo de exposição deve ser reduzido pela metade. 

Além do ruído, a vibração é considerada um dos maiores riscos físicos de 
acidente na Colheita florestal. A vibração ocorre em função do movimento oscilatório do 
corpo, proveniente de forças de componentes rotativos ou alternados de máquinas ou 
equipamentos (CUNHA, 2001). 

Estudos têm apontado que há uma maior probabilidade de perda auditiva em 
trabalhadores expostos ao ruído e a vibração de maneira conjunta (SORANSO et al., 2019). 
Nos estudos, realizados por Turcot et al. (2015), foram verificadas maiores alterações 
auditivas ao longo do tempo entre os trabalhadores florestais que também sofreram da 
doença dos dedos brancos, induzida pela vibração de mãos e braços. 

Riscos ergonômicos associados à colheita florestal 


Os riscos ergonômicos são aqueles relacionados com os fatores fisiológicos e 
psicológicos inerentes à execução das atividades laborais (BARBOSA et al., 2014). Estes 
fatores podem produzir alterações no organismo e estado emocional dos trabalhadores, 
comprometendo a sua saúde, segurança e produtividade. Entre os principais fatores 
ergonômicos, incluem-se a monotonia, a posição e o ritmo de trabalho, a fadiga, a 
preocupação, os trabalhos repetitivos, entre outros (VOSNIAK et al., 2010). 

Como na atividade de colheita florestal, as atividades são realizadas com a utilização 
contínua de máquinas, equipamentos e ou ferramentas, estas exigem um elevado esforço 
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físico e a adoção de posturas, muitas vezes inadequadas (VOSNIAK et al., 2010). Por 
esse motivo, as atividades realizadas nesse setor, na sua grande maioria, podem ser 
classificadas como moderadamente pesadas (BARBOSA et al., 2014). 

Além dos pesos e ritmos excessivos, existem agravantes, como os controles rígidos 
da produtividade, as jornadas de trabalho prolongadas, a monotonia e a repetitividade, o 
estresse físico e ou psíquico (BRASIL, 2021d). As lombalgias por excesso de peso e má 
postura são comuns nessa atividade. 

A declividade do relevo é um fator determinante no dimensionamento dos riscos e 
sua gravidade, com ênfase no aumento dos riscos ergonômicos (BARBOSA et al., 2014). 
Segundo Fiedler et al. (2010), boa parte dos acidentes de trabalho na colheita florestal é 
decorrente dos traumas por esforços excessivos. Eles são responsáveis pela maior parte 
de afastamentos dos trabalhadores, em consequência das doenças e lesões no sistema 
músculo-esquelético (IDA; BUARQUE, 2016). 

Nesse contexto, a realização de estudos ergonômicos é imprescindível, pois este 
auxilia na harmonização do sistema de trabalho, adaptando-o ao ser humano, por meio 
da análise da tarefa, da postura e dos movimentos do trabalhador, assim como de suas 
exigências físicas e psicológicas (BARBOSA et al., 2014). O conhecimento desses fatores 
é de fundamental importância para que a área de trabalho, o seu arranjo, as máquinas, 
equipamentos e ferramentas sejam bem adaptados às capacidades psicofisiológicas, 
antropométricas e biomecânicas do ser humano (VOSNIAK et al., 2010). 


Riscos de acidentes associados à colheita florestal 


As atividades florestais, por características próprias, já despertam atenção constante 
no que tange à proteção de seu trabalhador (FIEDLER et al., 2010). O risco de acidentes 
é iminente em todas as atividades, tendo em vista a necessidade frequente do uso de 
diversos tipos e condições de máquinas, ferramentas e ou equipamentos (SOUZA et al., 
2010). 

Na fase de derrubada das árvores, os processos de desgalhamento e traçamento 
são também de alto risco, visto que a motosserra funciona em sua rotação máxima e sem 
apoio, estando sujeita a resvalos, podendo atingir o trabalhador. Quanto ao empilhamento 
manual, é uma atividade extremamente pesada, sobrecarrega a coluna lombar e exige do 
trabalhador a admissão de posturas estereotipadas (SOUZA et al., 2010). 

Em outra vertente, de acordo com Coelho, Lima e Fonseca (2017) muitos acidentes 
com motosserras acontecem pela falta de qualificação e profissionalização adequada dos 
trabalhadores, visto que o conhecimento para o exercício da função é repassado por um 
colega de trabalho ou por instrutores não devidamente treinados e ou qualificados. A isso, 
somam-se como necessárias as atualizações e reciclagem periódicas dos trabalhadores da 
colheita florestal, quanto às normas regulamentadoras, prevenção de acidentes e técnicas 
mais seguras de trabalho (BARBOSA et al., 2014). 
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Na colheita florestal, a topografia do terreno pode aumentar a quantidade e a 
gravidade dos riscos de acidentes envolvidos no trabalho (BARBOSA et al., 2014). Apesar 
da evolução tecnológica na colheita florestal em área declivosas, onde se utilizam máquinas 
sofisticadas, muitas importadas, a motosserra ainda continua sendo largamente utilizada, 
isto devido ao corte de árvores com motosserra permitir uma boa produtividade individual e 
poder ser feito em locais de difícil acesso (PESCADOR et al., 2013). 

Entretanto, já se encontram no mercado, máquinas profissionais com dispositivos 
de segurança para esse tipo de terreno, tais como, freio da corrente, pino pega corrente, 
sistema antivibratório, protetores da mão esquerda e direita, trava de segurança do 
acelerador, direcionador de serragem e escapamento com silencioso e direcionador de 
gases, sendo estes dispositivos essenciais para que as empresas florestais contribuam 
para a segurança do operador (CANTO et al., 2007). 

Em sua quase totalidade, o local de trabalho no setor florestal encontra-se ainda 
distante das residências dos trabalhadores, obrigando-os a despender tempo e energia 
no trajeto, aumentando assim o risco de sofrer acidentes, dado aos meios de locomoção 
acessíveis aos trabalhadores (COELHO; LIMA; FONSECA, 2017). Diante desse cenário, 
Machado (2014), menciona que as empresas devem considerar como itens essenciais 
para a realização das atividades florestais: a formação de operadores abrangendo técnicas 
de operação, manutenção e segurança do trabalho; a observância de distância entre os 
operadores para evitar a queda de árvores sobre alguém; a sinalização nos limites e nas 
proximidades dos talhões; a disponibilidade de material de primeiros socorros e veículo 
para locomoção de feridos na área de corte; e meios de comunicação eficientes. 


Normas regulamentadoras 
Atualmente, existem 37 normas regulamentadoras vigentes, disponíveis no 

Ministério do Trabalho e Previdência Social. Essas normas regulamentadoras, além de 
estabelecerem as condições mínimas de segurança para a realização das atividades, 
servem como referência para a fiscalização do trabalho. As normas importantes para o 
desenvolvimento de atividades florestais, em geral, são: 

a) NR - 6 Equipamentos de Proteção Individual (EPI); 

b) NR - 17 Ergonomia; 


c) NR - 31 Segurança e Saúde no Trabalho na Agricultura, Pecuária, 
Silvicultura, Exploração Florestal e Aquicultura. 


NR - 6 Equipamentos de proteção individual — EPI 


Equipamento de Proteção Individual (EPI) é todo dispositivo ou produto de uso 
individual que deve ser utilizado pelo trabalhador, com a finalidade de protegê-lo de 
situações que possam ameaçar sua saúde e segurança no trabalho (BRASIL, 2021). 
Portanto, o EPI não evita a ocorrência de acidentes, mas ameniza suas consequências, 
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evitando lesões ou reduzindo sua gravidade (SILVA, 2013). 

O uso de EPI está previsto na legislação trabalhista, ou seja, na Consolidação das 
Leis do Trabalho (CLT), mais precisamente, a Lei 6514 de 1977, que define: 

“Cabe ao empregador quanto ao EPI, adquirir o adequado ao risco de cada 
atividade, exigir seu uso, fornecer ao trabalhador somente o aprovado pelo órgão nacional 
competente, orientar e treinar o trabalhador sobre o uso adequado, guarda e conservação.” 

Segundo Meireles e Oliveira (2016), os funcionários quando não são bem instruídos 
e treinados no uso do EPI, em geral, afirmam que os riscos a que se expõem são pequenos, 
que já estão acostumados e sabem como evitar o perigo, sendo que o uso de EPIs é 
incômodo e limita os movimentos. Segundo a CLT: 

“Cabe ao empregado quanto ao EPI, usar, utilizando-o apenas para finalidade a que 
se destina, responsabilizar-se pela guarda e conservação, comunicar qualquer alteração 
que o torne impróprio para o uso, cumprir as determinações do empregador sobre o uso 
adequado.” 

A não utilização dos Equipamentos de Proteção Individual para cada atividade 
específica ocorre, muitas vezes, pelo desconhecimento, ou por tornar-se desconfortável 
para o trabalhador, além da ocorrência de negligência por parte das empresas em treinar e 
fornecer os equipamentos de proteção para seus trabalhadores (SILVA, 2013). Entretanto, 
é imprescindível a sua utilização, uma vez que ao realizar as atividades de forma segura, O 
trabalhador estará protegido contra possíveis riscos de acidentes (ASSUNÇÃO: CAMARA, 
2011). 

Na colheita florestal, o capacete é essencial para a proteção do impacto proveniente 
da queda de galhos (CANTO et al., 2007). Segundo os mesmos autores, o protetor facial 
é destinado à proteção dos olhos e da face contra lesões ocasionadas por partículas de 
madeiras resultantes das atividades de derrubada, desgalhamento e traçamento. Ademais, 
as luvas são necessárias principalmente para proteção das mãos contra agentes abrasivos, 
escoriantes, cortantes, perfurantes, como farpas de madeiras nas atividades de extração e 
carregamento manual e semimecanizado. 

Botas de cano longo ou botinas com caneleira com malha de aço, além de serem 
fundamentais para a proteção de operadores de motosserras, são importantes para 
trabalhadores do campo, que estão sujeitos a acidentes com animais peçonhentos (CANTO 
et al., 2007). O uso do protetor auricular também se faz útil para proteger o sistema auditivo 
contra níveis de pressão sonora nas atividades de corte, pois os operadores de motosserra 
estão sujeitos à influência de um nível de ruído superior ao permitido pela legislação 
brasileira, o qual, para uma jornada de trabalho de oito horas, é de 85 dB, conforme a NR- 
15 (BRASIL, 2021c). 
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NR — 17 Ergonomia 


A Norma Regulamentadora Nº 17 visa adaptar as condições de trabalho às 
características psicofisiológicas dos trabalhadores, proporcionando um maior conforto, 
segurança e desempenho (BRASIL, 2021d). Assim, ela procura reduzir a fadiga, estresse, 
erros e acidentes, proporcionando segurança, satisfação e saúde aos trabalhadores, 
durante o seu relacionamento com esse sistema produtivo. 

Segundo esta norma, o empregador é responsável por fazer uma análise ergonômica 
do trabalho, a fim de adaptá-lo aos trabalhadores, devendo a mesma abordar, no mínimo, 
as condições de trabalho, conforme especifica a norma (BRASIL, 2021d). 

Dutra, Leite e Massad (2012) afirmam que a adaptação do homem ao trabalho é 
bastante complexa, sendo que a ergonomia poderá contribuir a partir do conhecimento 
do ser humano, projetando o ambiente de trabalho e ajustando-o às suas capacitações 
e limitações. Por isso, na realização de uma análise ergonômica é importante iniciar pelo 
levantamento dos fatores humanos e das condições de trabalho, o que permite conhecer 
o perfil dos trabalhadores e identificar as condições de realização do trabalho (BRITTO et 
al., 2015). 

Um local de trabalho deve ser sadio e agradável, proporcionando o máximo de 
proteção, sendo o resultado de fatores materiais ou subjetivos, que previne acidentes, 
doenças ocupacionais, além de proporcionar melhor relacionamento entre a empresa e 
o empregado (BRITTO et al., 2015). Por isso, para que as organizações possam alcançar 
o sucesso, devem buscar constantemente propostas para a melhoria do aperfeiçoamento 
dos métodos e das técnicas de trabalho, assegurando condições mais confortáveis, 
seguras e saudáveis ao ser humano e, consequentemente, aumentando a produtividade e 
a qualidade do trabalho (SANT'ANNA; MALINOVSKI, 2002). 

Dentre os principais fatores ergonômicos relacionados às atividades de colheita 
florestal, os biomecânicos (envolvendo as posturas, as forças aplicadas, a carga de trabalho 
físico e os movimentos repetitivos) têm influência direta sobre a saúde do trabalhador e, 
consequentemente, sobre a eficiência da operação, os quais podem ser minimizados por 
meio de modificações no método de trabalho e treinamentos específicos, com a finalidade 
de adoção de posturas mais seguras, saudáveis e confortáveis (IDA; BUARQUE, 2016). 

David, Fieldler e Baum (2014) concluíram que na colheita semimecanizada no 
Brasil, os principais riscos à segurança e saúde ocupacional correspondem às patologias 
ligadas à coluna vertebral, como lombalgia e hérnia de disco. A mecanização proporciona 
melhores condições no trabalho, porém, os riscos à saúde e segurança permanecem 
devido à repetitividade e monotonia da atividade (BARBOSA et al., 2014). Sua má condução 
propicia estresse e doenças psicossociais ao operador de máquinas, ainda estando sujeito 
a patologias musculoesqueléticas, como LER/DORT (Lesões por Esforços Repetitivos/ 
Distúrbios Osteomusculares Relacionados ao Trabalho) (SILVA et al., 2013). 
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NR-31 Segurança na agricultura, pecuária e exploração florestal e arquitetura 


Conforme consta na NR-31 (BRASIL, 2021e), cabe ao empregador ou equiparado 
“realizar avaliações dos riscos para a segurança e saúde dos trabalhadores e, com base nos 
resultados, adotar medidas de prevenção e proteção para garantir que todas as atividades, 
lugares de trabalho, máquinas, equipamentos, ferramentas e processos produtivos 
sejam seguros e em conformidade com as normas de segurança e saúde”. Essa norma 
é responsável por estabelecer os preceitos que devem ser abordados na organização e 
no ambiente do trabalho, voltada para as atividades da agricultura, pecuária, silvicultura, 
exploração florestal e aquicultura, podendo ser aplicada ainda em atividades de exploração 
industrial feita em estabelecimentos agrários. 

Segundo a NR-31, cabe aos empregadores cumprir e fazer cumprir os seguintes 
dispositivos legais (BRASIL, 2021e) 

* | Realizar avaliações dos riscos para a segurança e a saúde dos trabalhadores e, 


com base nos resultados obtidos, adotar as medidas de prevenção e proteção 
adequadas; 


* Analisar com a Comissão Interna de Prevenção de Acidentes no Trabalho Rural 
— CIPATR, as causas do acidente e as doenças decorrentes do trabalho, bus- 
cando prevenir e eliminar as possibilidades de novas ocorrências; 


* Assegurar que sejam fornecidas aos trabalhadores instruções compreensíveis 
em matéria de segurança e saúde, bem como toda orientação e supervisão 
necessária ao trabalho seguro; 


* | Adotar medidas de avaliação e gestão dos riscos, seguindo na ordem de eli- 
minação dos riscos, controle de risco na fonte, redução do risco ao mínimo, 
adoção de medidas de proteção pessoal no caso de persistirem os riscos. 

A prioridade, portanto, é a eliminação de riscos à saúde, não o fornecimento de 
EPI's. Este deve ser a última medida a ser tomada para a neutralização dos riscos, caso 
não possam efetivamente, ser eliminados, controlados na fonte ou reduzidos ao mínimo. 

As normas de segurança do trabalho, às vezes, são ignoradas por alguns 
trabalhadores, mas é de obrigação cumprir as determinações, adotar as medidas de 
proteção, e submeter-se aos exames médicos, colaborando assim, para que a NR-31 
seja cumprida. Também é um direito do trabalhador ser informado e consultado sobre as 
medidas de prevenção que serão adotadas pela empresa, podendo escolher, dentre as 
alternativas, aquela que lhe for mais conveniente (BRASIL, 2021e). 


CONSIDERAÇÕES FINAIS 


Dentre as atividades analisadas, as que apresentaram maiores riscos de acidentes 
na atividade de colheita florestal foram o corte e a derrubada de árvores, em especial com 


o uso de motosserra em áreas declivosas, agravando os riscos ergonômicos. Outros riscos 
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presentes na maioria das atividades ligadas à colheita florestal são os riscos físicos, de 
ruído e vibração sendo a declividade do relevo um fator determinante no dimensionamento 
dos riscos e sua gravidade. 

Amitigação desse cenário baseia-se na elaboração de planos de ação, determinando 
os riscos inerentes às atividades, de modo a: auxiliar em ações de treinamento, 
capacitação, fiscalização e atualizações periódicas, com foco no trabalho seguro dos 
envolvidos na colheita florestal; orientar para a utilização efetiva de todos os equipamentos 
de proteção necessários, seja individual, quanto coletiva, por todos os trabalhadores e; 
aplicar investimento sistemático em medidas de segurança e saúde dos trabalhadores e na 
disseminação de uma cultura prevencionista. 

Por fim, ainda são necessários mais estudos sobre avaliações quantitativas 
imediatas para os riscos físicos de vibração e ruído, de acidentes e análises qualitativas 
dos agentes ergonômicos, a fim de refinar as melhores ações mitigadoras. 
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RESUMO: A utilização de relações biométricas e 
volumétricas permitem determinar com exatidão 
a altura das árvores e realizar estimativas da 
prognoses da produção madeireira em uma 
floresta. O objetivo desse trabalho foi determinar a 
altura e o volume de seis espécies de Eucalyptus 
através de modelos de regressão. Os dados foram 
coletados na Estação Experimental de Ciências 
Florestais em Itatinga-SP, correspondendo uma 
área de 28,96 ha de 10 talhões de Eucalyptus 
spp. Quatro modelos hipsométricos foram pré- 
selecionados e ajustados. A seleção do melhor 
modelo seguiu os seguintes critérios: coeficiente 
de determinação ajustado (R2aj.), erro padrão 
da estimativa (Syx%), F calculado (Fc), 
coeficiente de variação (CV%) e análise gráfica 
da distribuição de resíduos. O melhor modelo foi 
o de Prodan (E= BO+Bl+*d+p2 + d?) para todos 
os talhões. Para os modelos ajustados de altura, 
o R2aj., Syx%, F... € CV % variou de (0,92 a 0,99); 
(2,18 a 22,19); (7,81 a 4073,88) e (4,12 a 31,47), 
respectivamente. A distribuição de resíduos 
mostrou uma diferençal para o modelo 5 em 
comparação aos outros modelos. O maior volume 
se encontra no talhão 79, que foi determinado 
pelo modelo ajustado: v = 0.0003 + d8298hh11712 
apresentando um volume de 486,13 m3.ha”. 
PALAVRAS - CHAVE: Modelos hipsoméftricos; 
Modelos volumétricos locais; incremento 


ABSTRACT: The use of biometric and volumetric 
relations allows to accurately determine the 
height of the trees and to make estimates of the 
prognoses of the wood production in a forest. 
The objective of this work was to determine the 
height and volume of six species of Eucalyptus 
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using regression models. The data were collected in the Forest Science Experimental 
Station in Itatinga-SP, corresponding to an area of 28.96 ha of 10 stands of Eucalyptus spp. 
Four hypsometric models were pre-selected and adjusted. The selection of the best model 
followed the following criteria: adjusted coefficient of determination (R2aj), standard error of 
estimation (Syx), calculated F (Fcalc.), Coefficient of variation (CV%) and graphical analysis 
of the distribution of residues. The best model was the Prodan model (E = B0O+B1l+d+p2x* ad? 
) for all stands. For the height-adjusted models, the R2aj., Syx, Fcalc. and CV (%) ranged 
from (0.92 to 0.99); (2.18 to 22.19); (7.81 to 4073.88) and (4.12 to 31.47), respectively. The 
distribution of residues showed a subtle difference for model 5 in comparison to the other 
models. The largest volume is in stand 79, which was determined by the adjusted model: 
v = 0.0003 + dt82º8hh11712 »resenting an volume of 486.13 m3.ha”. 

KEYWORDS: Hypsometric models; Local volumetric models; increment 


INTRODUÇÃO 


Uma das grandes vantagens do Eucalyptus é sua facilidade de cruzamento 
entre diferentes espécies, no qual o processo é denominado hibridação, resultando da 
combinação entre o Eucalyptus grandis e o Eucalyptus urophylla, o chamado Eucalyptus 
urograndis, um dos híbridos mais conhecidos e usados no Brasil. 

A principal finalidade da maioria dos inventários é obter a informação quanto 
ao volume, principalmente para fins comerciais, como os seus produtos madeireiros 
e a produção de madeira dos povoamentos florestais (Batista et al. 2014, Machado e 
Figueiredo Filho, 2014). A altura constitui-se em importante característica da árvore e pode 
ser medida ou estimada. Sua medição ou estimação é muito importante para o cálculo 
do volume, de incrementos em altura e, em determinadas situações, pode servir como 
indicadora da qualidade produtiva de um local (Silva et al., 2012). 

O método indireto de obtenção da variável altura (através de relações hipsométricas) 
das árvores é o mais utilizando pelos inventários florestais (Ferreira, 2018). Tais relações 
hipsométricas são equações que prediz a altura através da variável DAP, dos demais 
indivíduos dentro das parcelas, que não foram medidas as suas alturas em campo (Thiersch 
et al., 2013; Kershaw et al., 2016). Os modelos hipsométricos são equações ajustadas que 
expressam a relação altura-diâmetro da árvore (Sanquetta et al., 2014). Por meio desses 
pode-se obter estimativas acuradas da altura do indivíduo arbóreo. Com isso, o inventário 
diminui seu tempo na medição das árvores na parcela, o custo e as dificuldades na hora 
da operação, que torna os inventários mais precisos e econômicos (Atanazio et al., 2017). 

A hipótese aqui testada é de que independente da espécie e idade de Eucalyptus 
spp. aqui avaliadas o modelo matemático para determinação da altura total e do volume 
seja o mesmo para a área de estudo. 
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METODOLOGIA 


Para o estudo foi selecionado 10 talhões de plantio da espécie de Eucalyptus spp. 
(Tabela 1) correspondentes a uma área total de 28, 96 ha. O espaçamento dos talhões de 
Eucalyptus spp., variaram de acordo com o plantio dos mesmos. O povoamento pertence 
a Estação Experimental de Ciências Florestais do Departamento da Escola Superior de 
Agricultura Luiz Queiroz (ESALQ/USP), situado no município de Itatinga-SP (23º10' S e 
48º40' W e altitude de 857 m). 

O clima é classificado como Cwa, mesotérmico úmido, possuindo temperaturas que 
variam de 16,2 “C nos meses mais frios e 28,6 “C, nos meses mais quentes, conforme a 
classificação de Kôppen (Sette Júnior et al., 2010). A precipitação média anual nos períodos 
de outubro a março para EECFI corresponde a 1300 mm (Ferraz & Poggiani, 2014). O 
solo é do tipo Latossolo Vermelho Amarelo Distrófico típico A moderado com textura média 
(LVAd). O relevo se caracteriza por ser predominantemente plano, com uma declividade de 
2% a 10%, principalmente nas extremidades do talhão (Gonçalves et al., 2012). 

O procedimento de amostragem da área adotado foi estratificado por talhões, e 
as parcelas dentro dos talhões foram alocadas de forma sistemática. Foram realizadas 
a distribuição de parcelas nos respectivos talhões selecionados, oscilando de 387 a 534 
m2, pois as parcelas têm uma variação quanto ao tamanho, além das espécies dentro do 
mesmo, devido ao seu espaçamento. Foram alocadas de 1 a 3 parcelas compreendendo a 
área da parcela, de acordo com a tabela 1, no entorno da microbacia do córrego Tinga, de 
acordo com os talhões. 











Talhão Área (ha) Espécies Ra EM (m) RA o Ea LA. 
70 0,81 E. grandis 21 2,44x2,44 3 11,48 3 
62º 2,68 E. grandis — Coff”s Harbour 20 2,18x2,18 3 8,89 3 
64 142 E. grandis 20 2,19x2,19 3 815 0,37 1 
76º 54 E. urograndis I9 2,08x2,08 3 0,74 111 1 
79 9,56 E. urograndis x E. grandis I9 213x2,13 3 0,37 1,48 1 
139 174 E. grandis 8 220x2,20 3 Lilo - l 
86b 0,45 E. grandis 5 2,27x2,27 l Lil 1441 1 
151 0,86 E. grandis 5 2,14x 2,14 2 2,78 l 
166 232 E. grandis 4 213x2,13 3 10,00 1,85 1 
168 3,72  E.grandise E. clonal(1-144) 4 213x2,13 3 12,559 1,85 1 

Total 28,96 - - E 24 = E 





Tabela 1 - Espécies selecionadas da área da influência da MBH do córrego Tinga da EECF de Itatinga- 
SP e suas respectivas áreas (ha), idades (anos), espaçamentos médios (m), taxa de desbaste (%), 
taxa de mortalidade (%) e taxa de falha (%). 


Em que: Nº = número; D= desbaste; F= falha; M= mortalidade, T= talhão; I= idade (anos); 
E.M=espaçamento médio; 1.A. = intensidade amostral (%). Fonte: Autor (2018). 


As parcelas foram de formato retangular, sendo compostas por 6 linhas de 15 


indivíduos em cada linha, totalizando 90 indivíduos na parcela, que possuíam duas linhas 
de 15 indivíduos como área de borda, e entre as parcelas havia um intervalo de uma 
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parcela fantasma e a área de cada parcela é variável devido ao espaçamento de cada 
plantio e da regularidade do alinhamento do mesmo. 

Os dados diretamente coletados foram a altura e o DAP (diâmetro à altura do peito) 
a 1,3 m do solo, utilizando os equipamentos Vertex IV e fita métrica para a coleta dos 
mesmos respectivamente. Os modelos utilizados foram pré-selecionados de acordo com o 
uso mais atual na literatura florestal (Tabela 2). 








Modelos Modelos matemáticos Autores 
, h=Bo+Byj+d+B;+d? Trorey 
2 h = Bo +B;+*Lnd + B; + Lnd? Backman Modificada 
2 
3 C=fo+fird+fosd? Prodan 
d? 
4 >= = Bo +Bixd+Bp*d? Finger (1992) 


vh— 1,30 


Tabela 2 - Modelos hipsométricos para estimar a altura total (m) de Eucalyptus spp. da microbacia do 
córrego Tinga na EECF em Itatinga — SP. 





Em que: h = altura estimada em m; Ln = logaritmo neperiano; d = diâmetro a 1,30 m do solo em cm; B,, 
B, B, e B,= coeficientes de regressão 


Para a determinação do volume foi utilizado modelos volumétricos ajustados em 
relação ao local de origem, precisamente no Estado de São Paulo. Pela tabela 3, observa- 
se os modelos ajustados por cada talhão e região para determinar o volume com seus 
respectivos locais e autores (Filho et al., 2014). 





Talhão Modelos Coeficientes R? Syx(Y%) Local Fonte 
70 v=Bid2+e P10,000848 0,986 13,43 Itatinga, SP Silva (1996) 
62º v=fBid2+e [1000053 0,991 10,14 Itatinga, SP Silva (1996) 
64 v=Bid2+e B10,000633 0,987 12,32 Itatinga, SP Silva (1996) 
76º, 79, BO 0,00003 isca Mctague, 
139 v=PBodfihb2+e Pl 18298 0,65 a Batista e 
e 168 p2 1,1712 z a Steiner (1989) 
86º v=PBid2+e; P10,000563 0,988 Itatinga, SP Silva (1996) 
151 v=PBid+e B10,000591 0,992 Itatinga, SP Silva (1996) 
BO -13,153 Itapetininga, Angatuba, 
166  v=Bd2+e, B1-0,52295 0,988 11,84 Itatinga, Bofete, ore 
B2 89242 Botucatu e Pardinho, sp (1987) 





Tabela 3 - Modelos volumétricos ajustados para estimativa do volume total (mº) para as espécies de 
Eucalyptus spp. e suas respectivas estatísticas. 


Em que: v= volume estimado em m?; h = altura total (m), Ln = Logaritmo neperiano, d = diâmetro 
a altura do peito em cm; = B,, B,, B, coeficientes de regressão; e= exponencial; £, = termo de erro 
aleatório; R2= coeficiente de determinação; Syx = erro padrão da estimativa (%) e (m3). 


Na tabela 4 para seleção do modelo matemático que melhor descreve a estimativa 
do volume e da relação hipsométrica, foram avaliados o erro padrão residual (Syx%), os 
coeficientes de determinação (Rºaj.) , a dispersão gráfica dos resíduos e o valor F. 
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O erro padrão da estimativa em porcentagem informa a qualidade do ajuste e 
o quanto, relativamente, o modelo erra em média ao estimar a variável dependente. O 
coeficiente de determinação ajustado em porcentagem da %), O qual o coeficiente de 
determinação expressa à quantidade de variação total explicada pela regressão. Quanto 
mais próximo de 1 forem os valores dos coeficientes de determinação melhor será o ajuste 
da linha de regressão. Sendo assim, por este critério, seleciona-se o modelo matemático 
que apresentar o maior valor de R2, (Schneider e Schneider, 2008). O teste de F é uma 
estatística obtida por meio da análise de variância da regressão e indica se há significância 
no ajuste a determinado nível de probabilidade, através da comparação do valor de F 
calculado com o valor de F tabelado (Draper et al., 1989), se o valor de F calculado for 
maior que o F tabelado, a regressão analisada é considerada significativa ao nível de 
probabilidade desejada. 

A distribuição dos resíduos foi analisada plotando-se os desvios da regressão em 
gráfico, com o eixo das ordenadas centradas em zero. Observa-se a distribuição dos 
pontos e verifica-se se apresenta ou não tendências. Quando a distribuição dos resíduos 
da regressão apresenta-se homogênea, é considerado que a equação pertinente possui 








bom ajuste. 
Critério Equação Em que 
no ea SQ Regressão 
Coeficiente de determinação 2 = “SOToal 
Q Tota HEsf SQ = soma de quadrados; 
Coeficiente de determinação ajustado R2aj.= 1 — Ja —R2) (| n= número de observações; 
np p = número de coeficientes da 
p2 regressão; 
Erro padrão da estimativa Syx= Eh) h= de (m); 
a P o h = altura estimada (m); 
Erro padrão da estimativa (%) Syx% = idas h = altura média observada (m); 
E h sh = desvio padrão da altura observada. 
Coeficiente de variação CV% = E * 100 





Tabela 4 - Critérios de seleção para a escolha da melhor equação para estimativa da altura em 
Eucalyptus spp. dá microbacia do córrego Tinga na EECF em Itatinga — SP. 


O incremento médio anual (IMA) foi determinado pela fórmula (Finger ,1992): 
IMA = 2 
1 


Em que: IMA= Incremento Médio Anual (mº.ha''.ano); v= Volume (mº.ha""); I= Idade 
(ano). 
A intensidade amostral foi calculada pela fórmula (Soares et al. 2011): 


NS N-n 
Et 
NE2+t2Sy N 





Em que: n — número ótimo de parcelas; Y — número total de parcelas que cabem 
na área; Sr? —variância populacional; 4 — limite de erro admitido (0,10); t — valor de “t” 
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tabelado, da distribuição de Student, em função dos graus de liberdade. 

Foi realizado o Teste Tukey a fim de verificar a diferença estatística do volume e 
incremento médio anual. As análises foram feitas usando os softwares Excel, versão 16 
para a realização do cálculo estatístico, o SigmaPlot, versão 12 e Minitab 18 para gerar os 
gráficos. 


RESULTADOS E DISCUSSÃO 


Por meio dos valores mínimos e máximos de Hte DAP (tabela 5) é possível perceber 
que os maiores valores de IMA estão nos povoamentos mais jovens por estarem na fase 
inicial de crescimento. Os menores IMA em DAP estão nas maiores idades indicando a 
necessidade de intervenções silviculturais. O teste de tukey indicou diferenças significativas 


para ambas as variáveis analisadas. 


Variável Média Variância CVY% IMA Mínimo Mediana Máximo Assimetria Curtose 
Talhão 70 - Eucalyptus grandis 
DAP (cm) 19,05º 43,295 34,53 0,908 10,064 21,290 33,439 «0,02 -1,45 
Ht(m) 27,27 46,576 25,03 1,299 17,385 30,110 40,074 0,15 -1,58 
Talhão 62º - Eucalyptus grandis — Coff”s Harbour 
DAP (cm) 26,43º 32455 21,55 1,322 12,134 26,083 40,159 0,16 0,04 
Ht(m) 30,79 8,676 9,57 1,540 26,263 30370 38,518 0,66 0,29 
Talhão 64 - Eucalyptus grandis 
DAP (cm) 2547" 35,115 23,26 1274 12,006 24,682 38217 0,06 -0,71 
Him) 3485? 24,962 14,33 1,743 17,582 35419 41,780 «0,92 0,56 
Talhão 76º - Eucalyptus urograndis x Eucalyptus grandis 
DAP (cm) 20,15º 33,892 28,88 1,061 10,510 19,809 39,172 0,95 119 
Ht(m) 25,57: 24,718 19,44 1346 9,335 25,889 36,519 -0,3 0,36 
Talhão 79 - Eucalyptus urograndis x Eucalyptus grandis 
DAP (cm) 211º 50,245 33,54 1112 10,032 19,745 54,618 1,2 3,03 
Ht(m) 27,82º 21,505 16.66 1,465 19252 27,276 44,042 0,46 0,15 
Talhão 139 - Eucalyptus 
DAP (cm) 13,62º 4,258 15,14 1,703 10,027 13,687 19,099 0,28 «0,42 
Ht(m) 28,90" 8,492 10,08 3,613 23,551 29,092 36,098 0,13 «0,58 
Talhão 86º - Eucalyptus grandis 
DAP (cm) 24,17" 46,91 28,33 4834 1136 25,020 38,77 0,1 «0,89 
Ht(m) 2540" 9528 38.43 5,080 9,86 25,730 38,1 «0,22 -1,43 
Talhão 151 - Eucalyptus grandis (estudo adubação TIMAC) 
DAP (cm) 15,94º 4,199 12,85 3,189 10,250 15,884 21,306 0,09 0,15 
Ht(m)  26,48º 0,804 3,39 5,297 26,044 26,208 33,330 5,23 35,44 
Talhão 166 - Eucalyptus grandis 
DAP (cm) 14,85º 7,642 18,61 3,714 6,146 14,650 2484] 0,26 0,7 
Ht(m) 21,38º 4,549 9,97 5.347 9,009 21,725 24,850 -1,94 7,78 
Talhão 168 - Eucalyptus grandis e Eucalyptus clonal (1 - 144) 
DAP (cm) 13,13º 8,539 22,24 3,285 5350 13,503 23,599 “0,09 1,6 
Ht(m) 2232º 13,34] 16,36 5.581 6,797 23,730 24314 -2,72 6,97 


Tabela 5 — Estatística descritiva para as variáveis DAP e Ht (m) 
Em que: DAP: diâmetro a altura do peito; Ht: altura total; IMA: Incremento médio anual. 


Médias seguidas pela mesma letra para a variável DAP e Ht não diferem entre si a um nível de 95% de 
probabilidade de confiança. 
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Na distribuição altimétrica, com os talhões apresentando diferentes caracteristicas 
como idade, espaçamento e especie. Os talhoes 70, 64, 139, 86b e 151 apresentaram 
curtose leptocurtica (coeficiente de curtose < 0,263), ou seja a distribuição é menos 


achatada (mais afilada) do que a curva normal de mesma área. Os demais talhões 62a, 


76a, 79, 166 e 168 apresentaram curtose platicurtica (coeficiente de curtose >0,263), ou 


seja, a distribuição é mais achatada do que a curva normal de mesma área (Figura 3). 


Aaltura apresentou valor médio de 21,76 m, variando entre 6,8 m (talhão 168) e 44,04 
m (talhão 79); a variação da média das alturas foi de 15,39 metros, variando entre 21,39 m 
(talhão 166) e 34,85 m (talhão 64). Para a variável altura, Silveira et al. (2014) encontraram 


resultados de 16m para Eucalyptus dunnii aos cinco anos de idade, significativamente 


diferente dos valores de altura deste experimento (Figura 2). 


Média 
33,72 
2949 
25,61 
2172 
2212 


Frequência 








Centro de Classe de Altura Total (m) 


20 


Variável 
E Talhao 86b 
EE Talhao 139 
EEE Talhão 151 
EEE Talhão 166 
EE Talhão 168 


DesvPad N 
B474 16 
5.750 44 
1,466 30 
3,558 44 
3,526 44 


Variável 
E Talhão 70 


2815 
3041 
35.54 
2743 
30,44 


Frequência 





30 
Centro de Classe de Altura Total (m) 


EEB Talhao 623 
FEZES Talhao 64 
EE Talhas 764 
EEE Talhão 79 


Média DesuvPad N 


8,773 42 
4693 43 
6,129 44 
6436 41 
5,893 40 


Figura 1 — Frequência do centro de classe de altura total para os talhões 86b, 139, 151, 166 e 168s 70, 


62a, 64, 76a e 79. 
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Para a variável DAP, o talhão 168 apresentou o menor valor, contando com a 
presença de indivíduos com DAP de 5,35 centímetros. Dentro da mesma variável, o talhão 
168 também apresentou o menor valor para a média do diâmetro a altura do peito, com 
valor médio de 13,13 centímetros. Ainda se tratando do DAP, o talhão 79 apresentou o 
maior valor para DAP (54,62 cm). Reiner (2014) encontrou para a variável DAP aos quatro 
anos de plantio, valores médios de 14,08 centímetros. No mesmo contexto, Vidaurre et al., 
(2015) também encontraram resultados diferentes de DAP (17,98cm) e altura total (24,64m) 
com o experimento em Guarapuava, aos cinco anos. Miguel et al., 2014 em um plantio 
de Eucalyptus grandis, no município de Niquelândia/GO obteve valores para diâmetros 
do povoamento, apresentando média de 13,60, mediana de 13,79. Valores inferiores ao 
encontrado no presente trabalho, com exceção do talhão 168 que apresentou os menores 
valores . 

Quanto à curtose para distribuição diamétrica, os talhões 70, 64, 76a, 79 e 139 
apresentaram curtose leptocurtica enquanto os demais talhões (62º, 86b, 151, 166 e 168) 
a curtose foi apresentada como platicurtica (Figura 2). 
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Variável 
E Talho 70 
EEB Talhao 625 
FEZES Talhao 64 

Talhao 764 


EE Talhão 79 


Média DesvPad N 
1906 6,580 164 
2644 5697 62 
2548 5,926 96 
2016 5622 99 
2113 7.088 125 


Frequência 








225 30,0 37,5 450 525 
Centro de Classe Media DAP (cm) 


Variável 
EE Talhoo b6b 
EEB Talhao 139 
EEE Talhão 151 
ESEE Talhão 166 
EEE Talhão 168 


Média DesvPad N 
2417 6849 43 
19.63 2063 216 
15,95 2049 99 
1486 2765 192 
1314 2922 232 


Frequência 





25 30 
Centro de Classe Media DAP (cm) 


Figura 2: Frequência do centro de classe média para o diâmetro a altura do peito para os talhões 70, 
623, 64, 76:e 79 86º, 139, 151, 166 e 168. 


Os valores dos critérios estatísticos (Rºaj., Syx%, F... e CV (%)) e os modelos 
ajustados com seus respectivos coeficientes para validação dos modelos hipsométricos 
encontram-se na tabela 6. De acordo com a tabela 6, o melhor modelo para estimar a altura 


foio 3 (modelo de Prodan) para todos os talhões em termos de coeficiente de determinação 


calc. 


ajustado e erro padrão da estimativa. Já o coeficiente de variação obteve menores valores 
para os modelos de 1 e 2 para todos os talhões e o F calculado para todos os talhões 
desempenhou melhores resultados para o modelo 4, apesar deste modelo ter os maiores 
valores de R2aj., este não é selecionado para estimar altura, devido seu Syx% e CV%, 
possuírem valores altos. 

Para Sousa et al. (2013), os valores de R2aj. encontrados foram semelhantes (entre 
0,40 a 0,43) para os modelos 1 e 2 do talhão 62a, em relação aos modelos hipsométricos 
ajustados para a espécie de Eucalyptus urophylla. O modelo de Prodan foi o melhor para 
o presente trabalho, de forma igual o mesmo foi selecionado em um inventário florestal 
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de Eucalyptus urophylla S. T. Blake, de acordo com Leal et al., 2011; Miguel. 2009. Pode 
se verificar que no trabalho de Atanazio et al. (2017), o R2aj. variaram de 0,55 a 0,74, 
corroborando com os resultados dos modelos de a 1 e 2 dos respectivos talhões 79 (0,73%), 
139 (0,715%) e 166 (0,78%). 


Coeficientes 
































Talhão Nº Bo B: po Rºaj. Syx% Fcalc. CV% 
l 2,6662 1,6332 -0,0155 0,91 9,28 211,95 29,83 
7% 2 10,7777 -7,6084 4,5750 0,91 931 210,39 29,82 
3 -0,2863 0,5058 0,0101 0,97 4,34 603,99 30,01 
4 -3,6852 1,7570 0,1116 0,99 14,01 3104,45 31,47 
l 29,7509 -0,3787 0,0157 0,45 10,80 17,54 9,99 
6a 2 103,8906 -57,4073 10,7374 0,44 10,85 17,15 9,93 
ii 3 -7,8148 1,0127 0,0056 0,95 7,71 380,93 10,88 
4 -23,4338 3,0196 0,1037 0,99 22,19 1808,53 11,32 
l 0,3768 1,9865 -0,0232 0,69 9,44 48,26 14,45 
64 2 -24,4896 17,9269 0,1905 0,68 9,61 45,81 14,34 
3 6,4795 -0,1325 0,0230 0,93 5,61 293,05 14,75 
4 19,2620 -0,3360 0,1528 0,99 18,22 1418,56 15,58 
l -1,4718 1,8843 -0,0238 0,89 7,15 164,28 22,21 
Toa 2 -22,0020 16,4912 -0,0831 0,89 7,85 158,84 22,18 
3 2,5177 0,2234 0,0200 0,97 5:51 688,99 22,22 
4 3,4298 0,9659 0,1415 0,99 16,27 3382,4 23,38 
1 11,0168 0,9732 -0,0069 0,73 10,03 54,23 16,72 
79 2 16,3265 -6,0023 3,2835 0,73 10,05 53,87 16,71 
3 -1,4745 0,5333 0,0134 0,98 6,05 776,56 17,30 
4 -8,7808 1,9430 0,1203 0,99 18,11 4073,88 17,75 
86b l -1,8010 2,0527 -0,0254 0,92 7,06 79,83 23,85 
2 -32,9243 24,1734 -1,1076 0,92 7,36 73, 23,78 
3 3,0719 0,1097 0,0203 0,96 5,47 186,19 23,85 
4 5,8267 0,6071 0,1418 0,99 17,56 841,78 25,10 
l 4,0621 23147 -0,0334 0,71 10,40 54,04 16,95 
139 2 3,7812 0,8822 3,4197 0,72 10,39 54,1 16,96 
3 -0,4022 0,3727 0,0093 0,92 2,18 234,79 17,34 
+ -2,5180 1,2756 0,1097 0,98 6,42 1247,93 18,04 
l 9,3989 1,5279 -0,0333 0,38 4,94 10,07 4,12 
151 2 -141,7874 112,4356 -18,7878 0,38 4,95 9,97 7,89 
3 10,8391 -1,1511 0,0686 0,96 2,34 310,01 5,58 
4 29,2804 -3,1040 0,2823 0,99 6,56 1178,32 8,50 
l -4,5950 2,1725 -0,0662 0,77 7,88 71,04 14,45 
166 2 -52,0528 45,5326 -6,7090 0,79 7,54 79,48 14,73 
3 3,9311 -0,2199 0,0427 0,96 4,09 521,8 14,75 
4 7,7314 -0,2795 0,2042 0,99 10,67 | 2273,28 15,65 
l -8,5728 3,8641 -0,1102 0,38 6,67 10,07 12,85 
168 2 -86,5347 77,2956 -13,4520 0,38 7,55 9,97 12,42 
3 6,9682 -0,8859 0,0693 0,96 2581 310,01 13,22 
4 17,0766 -2,1881 0,2786 0,99 7,32 1178,32 14,34 





Tabela 6 - Estatísticas e coeficientes de regressão para os modelos hipsométricos para espécies de 
Eucalyptys spp. 


Em que: Nº= número; h= altura total (m); d = diâmetro a 1,30 m do solo em cm; R2aj.= coeficiente de 
determinação ajustado (m); Syx= erro padrão da estimativa (%); Fcal. = F calculado e CV = coeficiente 
de variação (%). Fonte: Autor (2018). 


Na análise gráfica de resíduos observa-se que o modelo de Prodan foi o que 
apresentou melhor ajuste ao conjunto de dados dos povoamentos, evidenciando suas 
seleções como a melhor para cada um dos clones testados (Figura 3). 
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Figura 3: Distribuição de resíduo em função da altura total estimada (m) dos modelos hipsométricos 
ajustados. 


Na Figura 4, pode-se observar a relação entre as alturas estimadas e observadas 
das para as espécies estudadas de acordo com o conjunto de dados utilizado nos ajustes. 
Arelação gerou coeficiente de correlação de Pearson de 0,764; 0,901;0,821;0,781 e 0,843 
para E. grandis (Talhão 70), E. grandis (Talhão 62º), E. grandi (Talhão 64), E. urograndis 
(Talhão 76º) e E. urograndis x E. grandis(Talhão 79. Também os valores de 0,671; 0,892; 
0,911;0,581 e 0,977 para E. grandis (Talhão 86"),E. grandis (Talhão 139), E. grandis(Talhão 
151), E. grandis e E. grandis (Talhão 166), e E. clonal (Talhão 168), respectivamente. Esse 
valor é considerado moderado a alto, pois, conforme Cohen (1988), valores entre 0,5 e 1 
podem ser interpretados como grandes. 

As curvas estimadas apresentaram comportamento esperado para tal relação, 
verificando-se ausência de tendenciosidade. A inclinação da curva fornece indícios de que 
este povoamento está em declínio nos talhões 70, 62º, 64, 76º e 79, pois visto que segundo 
Bartoszeck et al. (2004), quando a inclinação da curva é achatada, trata-se no geral de 
povoamentos mais velhos, onde o incremento corrente anual, principalmente em altura, 
está em declínio. 
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Figura 4: Comportamento da altura observada e estimada para o Eucalyptus spp. — para o modelo 3 
para os talhões 70, 622, 64, 76: e 79, 86”, 139, 151, 166e 168. 


Na Tabela 7, o maior volume médio se encontra no talhão 79 (486,13 m3.ha”!), e 
menor volume para o talhão 151 (184,26 m3.ha”). O IMA para o talhão 622, obteve um 
resultado inferior aos demais talhões, sendo de 9,38 m3.ha”'.ano”*, fato este, que é explicado 
pelo desbaste da área e o maior se deve aos talhões mais jovens. A média nacional que 
é de 35,7 m3.ha” segundo IBÁ (2017). O volume médio variou de 350,17 a 459,50 m?.ha” 
no trabalho de Behrenz et al. (2016) para a espécie de Eucalyptus grandis pertencente 
à região do Médio Alto Uruguai do Estado do Rio Grande do Sul. Santos et al. (2017) 
apresentaram o valor médio de 293,5 mº.ha”, valor próximo ao obtido neste trabalho que foi 
de 297,93 m3.ha”. O incremento médio anual por hectare aumentou com o adensamento do 
espaçamento para 80% dos estratos observados, os talhões 62º e 76º não seguiram essa 
tendência (tabela 7). Isto ocorreu devido um aumento no número de árvores por área em 
menores espaçamentos, sendo o mesmo comportamento observado por Múller e Couto 


(2006) que estudaram diferentes espécies e idades de Eucalyptus. 





Talhão Area(ha) Idade (anos) EM (m) Volume médio(mº.ha”!) IMA (mº.ha'!) 
151 0,86 5 2,14x2,14 184,26º 36,85 
62º 2,68 20 2,18x 2,18 187,78" 9,38 
166 3,72 4 2,13x2,13 212,19” 53,05 
168 2,32 4 213x2,13 258,56º 64,64 
139 1,74 5 2.27 X2;2] 319,35º 39,92 
64 1,42 20 2,19x 2,19 319,954 15,24 
76º 54 19 2,08 x 2,08 327,42º 17,23 
86b 0,45 8 2,20 x 2,20 330,66º 22,02 
70 0,81 21 2,44 x 2,44 353,01 16,81 
79 9,56 19 2,13x 2,13 486,13 25,59 


Tabela 7 - Dados de Volume médio, Incremento médio anual, idade e as áreas que abrange cada 
talhão para as espécies de Eucalyptys spp. 


Em que: E.M = espaçamento; IMA=incremento médio anual (m3.ha”!) Fonte: Autor (2018). 
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No trabalho de Satiler (2014), o IMA foi de 68,63 m3.ha”, para o Eucalyptus urophylla 
x Eucalyptus grandis na região Norte de Mato Grosso, com idade de 4 anos, sendo 
semelhante ao do trabalho, precisamente para o talhão 168. Para Santos et a/., (2012) que 
estudou clones de Eucalyptus com idade de 48 e 72 meses em Minas Gerais encontraram 
o IMA variando de 40 e 44,51 m3.ha-?, onde o talhão 139 resultou dentro desta faixa de IMA. 
Para Sattler (2014), os resultados do trabalho não corroboraram, com o estudo de um teste 
clonal realizado em Juiz de Fora e Ubá (MG) 60 e 72 meses, com IMA de até 90,9 e 82,3 
m3.ha-!.ano-!, e seus volumes médios de 63 e 64,6 mº.ha-?.ano-!, respectivamente, porém 
o volume do talhão 168 se aproximou do trabalho citado. Cortez et al., (2008), afirma que 
a produtividade de 60 a 80 m3. ha-!. ano-! são de espécies melhoradas geneticamente, de 
produção de híbridos e clonagem, sendo que apenas o talhão 168 correspondeu a faixa de 
produtividade. 


CONCLUSÃO 


Diante dos resultados obtidos para cada espécie, a equação que apresentou 
resultado superior de ajuste e precisão foi a equação de Prodan sendo esta recomenda 
para estimativa de altura para todas as espécies, idades e talhões estudados na EECFI 
validando a hipótese aqui testada. 

Os maiores incrementos em DAP (cm) e Ht (m) se deram nos talhões 151,166 e 168 
composto pela espécie Eucaliptus grandis e Ecaliptus clonal. Sendo estes indicados para 
o plantio na região, sendo que a altura das árvores apresentou pouca variação entre os 
espaçamentos estudados (sendo maior no espaço médio de 4,7 m?). Também a maior IMA 
em volume (mº) se deu nos talhões 168 seguido do 166. 
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RESUMO: Com a crescente busca por novas 
fontes de óleos vegetais para a síntese do 
biodiesel, estudos de valorização e aplicação 
dos subprodutos resultante do processo de 
extração do óleo, como a torta e/ou farelo, tem 
despertado o interesse das indústrias de nutrição 
animal devido ao seu elevado. Surge então como 
alternativa, a aplicação desse subproduto como 
ração na alimentação de ruminantes. Essa é 
uma prática sustentável para o aproveitamento 
de um subproduto que é comumente descartado 
incorretamente, além de ser uma alternativa de 
baixo custo e eficaz de enriquecimento nutritivo 
na ração. O objetivo dessa revisão bibliográfica 
foi avaliar o potencial da torta de Pachira aquatica 
Aubl. como fonte de nutrientes em ração de 
ruminantes comparando com subprodutos de 
oleaginosas já empregadas na nutrição animal. 
PALAVRAS - CHAVE: Biodiesel, Pachira 
aquatica aubl., energias renováveis, biomassa 
florestal, munguba, ruminantes. 
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RUMINANTES 


ABSTRACT: With the growing search for new 
sources of vegetable oils for the synthesis of 
biodiesel, valuation studies and application of 
by-products resulting from the oil extraction 
process, such as cake and / or bran, has aroused 
the interest of the animal nutrition industries 
due to the its high. Then, as an alternative, the 
application of this by-product as a feed in the 
feeding of ruminants arises. This is a sustainable 
practice for the use of a by-product that is 
commonly discarded incorrectly, in addition to 
being a low-cost and effective alternative to 
nutritional enrichment in the feed. The purpose of 
this literature review was to assess the potential 
of Pachira aquatica Aubl pie. as a source of 
nutrients in ruminant feed compared to oilseed 
by-products already used in animal nutrition. 
KEYWORDS: Biodiesel, Pachira aquatica aubl., 
renewable energies, forest biomass, munguba, 
ruminants. 


11 INTRODUÇÃO 


O uso de combustíveis fósseis como fonte 
energética tem aumentado consideravelmente 
nas últimas décadas. No entanto, estudos de 
prospecção realizados no final do século XX 
revelaram que as principais reservas de petróleo 
do planeta deverão se esgotar em cerca de 100 
anos (RAMOS et al., 2011). Segundo Howard e 
Elisabeth Odum (2013), os combustíveis fósseis 
não se formam tão rapidamente quanto são 
consumidos. Vale salientar que os combustíveis 
fósseis além de ser uma fonte esgotável, é um 


dos principais responsáveis pela emissão de 
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gases que aceleram o aquecimento global. Então, em 1997 foi assinado o Protocolo de 
Quioto, que propõe que mais de 90 países se mobilizem para promover uma ação conjunta, 
onde o intuito principal é estabilizar na atmosfera a concentração dos gases causadores 
do efeito estufa (Greenpeace International, 2003). Contudo, há o constante interesse na 
utilização de fontes renováveis de energias, como a energia solar, eólica e biomassa. 
Sendo a última a principal fonte de biocombustíveis, dentre eles o biodiesel. 

O biodiesel é um biocombustível produzido a partir de óleos vegetais e/ou gordura 
animal (triglicerídeos), e pela sua natureza orgânica, é considerado um combustível limpo e 
renovável. Além de representar uma alternativa essencial para atender à crescente demanda 
energética da sociedade moderna de forma sustentável (RAMOS etal., 2011). Staiss e Pereira 
(2001) afirmam que levando em consideração apenas a biomassa proveniente de atividades 
agroindustriais, calcula-se que o potencial combustível desse material seja equivalente 
a, aproximadamente, 6.587 milhões de litros de petróleo ao ano. Com a visibilidade e 
vantagens agregadas ao biodiesel, foi deferido que o óleo diesel comercializado em todo o 
Brasil deveria conter 10% de biodiesel (Resolução nº4/2021, CNPE). 

Uma das biomassas que vem sendo explorada para essa finalidade é a Pachira 
aquatica aubl.(munguba). A munguba é uma espécie arbórea da familia das Malvaceae 
(Bombaceae), com crescimento rápido e característico de vegetações pantanosas 
(CAMACHO et al., 2017), é uma espécie nativa da Floresta Amazônica e com ampla 
distribuição nas Américas (LORENZI, 2000; SANTOS et al., 2007; LIMA et al., 2012; 
CAMACHO et al, 2017). O fruto dessa espécie contém amêndoas organolépticas, ou 
seja, com alto teor lipídico (LORENZI, 2008; SILVA et al., 2010; JORGE e LUIZA, 2012; 
SILVA et al., 2012; CORREIA, 2019) (Fig.1). Com isto, essa espécie tem se destacado 
tanto na indústria de biodiesel quanto em pesquisas científicas para a obtenção de mais 
conhecimentos sobre suas propriedades e características. 





Figura 1. Fruto aberto da Pachira aquatica aubl. com as amêndoas expostas. 


Fonte: Anjos et al., 2020 
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Com o aumento do uso de óleos vegetais para fins de biodiesel, ocorre maior 
demanda para aplicação do subproduto resultante do processo de extração do óleo, como 
as cascas dos frutos e tortas de sementes. Geralmente, a torta ou farelo gerado na extração 
do óleo não passam por processo de agregação de valor porque são desconhecidas as suas 
potencialidades nutricionais e econômicas, salvo algumas exceções como soja, algodão e 
girassol (ABDALLA et al., 2008). Esses subprodutos são caracterizados pelo elevado teor 
de proteína e podem ser aplicados na alimentação animal. 

Há alguns estudos desenvolvidos a respeito de tortas resultantes da extração do 
óleo de sementes oleaginosas, como a da mamona por exemplo. Segundo Fernando Pivetti 
(2013), a torta de mamona, é um subproduto do processo da extração de óleo da semente 
da mamona, possui alto teor protéico e aparece como alternativa para a alimentação animal 
e para a produção de materiais biodegradáveis, destacando que a mamona contêm uma 
toxina (Fernando Pivetti, 2013). Porém, no caso da torta da munguba, mesmo apresentando 
potencial para aplicação na ração de ruminantes, ainda não há resultados específicos sobre 
a influência desta no desenvolvimento e na nutrição desses animais. 

A extração dos óleos vegetais é realizada por dois métodos muito conhecidos, 
a extração mecânica ou extração química. A forma na qual ocorre a extração do óleo 
influencia diretamente no subproduto. Por extração mecânica tem-se como subproduto a 
torta, e por extração química, o farelo. A extração mecânica consiste em utilizar uma prensa 
(seja ela continua ou hidráulica) para a extração do óleo. As prensas contínuas são dotadas 
de uma rosca ou parafuso sem fim que esmaga o material, liberando o óleo (SILVA, 2018). 
As prensas hidráulicas (prensagem descontínua) apresentam um cilindro perfurado onde 
se desloca um êmbolo que faz pressão na matéria-prima que fica dentro de um saco de 
pano ou lona (SILVA, 2018). Quando a extração é realizada por este método, é constatado 
que na torta ainda fica um teor lipídico residual que pode ajudar na melhor disgestibilidade 
do animal (CARDOSO, 2013). Na extração química, também conhecida como extração por 
solvente orgânico, os grãos são triturados para facilitar a penetração do solvente (hexano 
— derivado de petróleo, éter etílico, etanol, metanol, entre outros) (SILVA, 2018). Os óleos 
migram das sementes para o solvente por terem maior afinidade com este. Esse método de 
extração se mostra mais eficiente em retirar o máximo de óleo possível, no entanto requer 
solventes químicos que elevam o custo do processo e alteram a composição da torta. 

Contudo, a aplicação desses subprodutos na alimentação animal se torna uma 
alternativa de baixo custo e eficaz para valorização destes, ao invés do seu descarte que 
muitas vezes ocorre de forma incorreta, tornando esta uma prática sustentável. 


21 MATERIAL E METÓDOS 


Esta é uma pesquisa teórica do tipo bibliográfica, baseada na vasta revisão de 
literatura com consulta de dados obtidos de outros autores. 
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A pesquisa foi realizada envolvendo artigos, teses, dissertações e dentre outros 
materiais acadêmicos, em plataformas como o Scielo, Google Acadêmico e outras bases 
a partir do Periódico CAPES. Com a finalidade de produzir um referencial teórico sobre a 
utilização de subprodutos de tortas oleaginosas na nutrição animal, tendo como enfoque 
principal, o potencial da torta de Pachira aquatica aubl. na alimentação de ruminantes. 


31 REVISÃO DE LITERATURA 


3.1 Características das espécies oleaginosas 


É essencial o conhecimento sobre a biomassa que será utilizada, tanto para o 
processo de biodiesel, quanto para aplicação do subproduto na alimentação de ruminantes, 
cada espécie tem suas particularidades. Em relação a composição das espécies oleaginosas, 
há a variação da quantidade de lipídeos, proteínas, fibras, carboidratos e dentre outros 
componentes. Além dos aspectos relacionados à composição e digestibilidade, também é 
importante observar a oferta frente a localização geográfica e a economia (SOUSA et al., 
2020). 


Subprodutos de espécies oleaginosas utilizadas na nutrição animal. 


Babaçu 


O babaçu (Orbygnia speciosa) é uma espécie de palmeira da família das Palmae, 
nativa da Região Norte e é comummente encontrada em regiões de cerrado (SOUSA et 
al., 2021). Consegue atingir até 20 m de altura, o que a classifica como de grande porte. 
É uma espécie muito explorada por conseguir utilizar-lá por completo, desde suas folhas 
até suas amêndoas, sendo as amêndoas de maior interesse para produção de biodiesel 
(SOLER et al., 2007). 

Com a extração do óleo de babaçu, vários subprodutos são produzidos, sendo um 
deles o farelo de mesocarpo de babaçu (FMB), também denominado farelo de amido de 
babaçu (SOUSA et al., 2021). Antes, esse subproduto seria descartado sem agregação 
de valor. Porém, com o tempo, devido a disponibilidade dessas espécie regionalmente, 
tornou-se uma fonte de enriquecimento nutricional na alimentação de ruminantes, mesmo 
sem respaldo científico (PORRO, 2019; SOUSA et al., 2021). 

Baseado na aplicação desse material na alimentação de ruminantes com apenas o 
conhecimento empirico, alguns pesquisadores realizaram estudos para melhor entendimento 
e exatidão da influencia desses subprodutos na alimentação dos ruminantes. FREITAS 
et al., (2014) realizou um experimento, onde foram aplicados diferentes quantidades do 
farelo de babaçu utilizada na alimentação de ruminantes, sendo 0%, 7,34%, 14,54%, e 
21,70%. O principal objetivo do estudo era avaliar o uso de farelo do coco babaçu como um 
suplemento na alimentação de ruminantes. Na Tabela 1 a seguir, é possivel verificar como 
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se deu a inclusão do FMB e de que forma se relacionou com a digestibilidade dos animais. 





Inclusão Farelo de 


babaçu IMS (g/kg ") DapMS (%) DapMO(%) 
(% ração) 
0 0,084 b 61,77 a 65,29 a 
7,34 0,089 ab 57,39 a 60,30 ab 
14,54 0,089 ab 58,27 a 64,18 ab 
21,70 0,092 a 58,09 a 59,70 b 
Blocos 
1 Rodada 0,089 a 59,45 a 62,30 a 
2 Rodada 0,087 a 58,30 a 62,44 a 
CV (%) 4,24 4,57 5,38 


Tabela 1. Inclusão de farelo de babaçu (% ração), ingestão de matéria seca (IMS), digestibilidade 
aparente da matéria seca (DapMS) e da matéria orgânica (DapMO) (Freitas et al., 2014). 


Fonte: Freitas et al., 2014. 


Como pode-se observar, não há valor nutricional relativamente alto com o aumento 
desse subproduto na ração. Concluindo nesse experimento, que deve-se restringir a 
inclusão na dieta não excedendo 20 % de ingestão de matéria seca (IMS), pois acarretar 
na perca de digestibilidade. 

SOUSA et al. (2021) realizou um experimento semelhante, cujo o objetivo principal 
foi avaliar a fermentação in vitro (IVRF) e a digestibilidade in vitro aparente (AIVD) de 
dietas formuladas com quatro níveis de inclusão de farelo de mesocarpo de babaçu (FMB) 
em ruminantes. Nesse experimento foram colocados 0%, 7,5%, 15%, e 22,5% de FMB 
na ração balanceada dos animais. Na Tabela 2 a seguir, pode-se observar a composição 
química do FMB reportado por Sousa et al., 2020. 








tem FMB 
Matéria seca (DM %) 87,83 
Cinza (% de MS) 5,36 
Proteína bruta (% de DM 6,09 
Extrato de éter (% de DM) 1,29 
Fibra de detergente neutra ( % de MS) 60,20 
Fibra de detergente ácido (% de MS) 44,02 
Hemicelulose (% de DM) 16,17 
Celulose (Y%de DM) 26,07 
Lignina (% de DM) 21,16 
Nitrogênio insolúvel em detergente neutro (N) (% do N total) 0,25 
Nitrogênio insolúvel em detergente ácido (N) (% de N total) 0,31 
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Fracionamento de carboidratos 





B2 (% de C (% de 


TOCH (% de DM) A+ B1 (% de TOCH) TOCH) TÓCH) 


88,25 36,42 6,05 57,53 





Tabela 2. Composição química do farelo do mesocarpo de babaçu (FMB) (Sousa et al., 2020). 


TOCH: Carboidratos totais; informaçõe fornecidas pelo provedor FMB. Fonte: Sousa et al., 2020 


Como apresentado na Tabela 2, o farelo é composto em fibras com altor de lignina e 
celulose, além de apresentar nível representativo de proteína bruta. Em relação ao extrato 
étereo (EE), é possivel observar que esta presente em pequena quantidade comparado ao 
teor de lignina, o que torna esse subproduto não tão viável em questões nutricionais, sendo 
que o nível elevado de EE é uma característica importante para a melhor digestibilidade 
do animal. Segundo Van Soest (1994), à lignina está inversalmente relacionada a 
digestibilidade. Além das características citadas anteriormente, o FMB é rico em amido, o 
que o torna uma ótima alternativa de fonte energética, pelo elevado teor de carboidratos 
(SOUSA et al., 2014; SOUSA et al., 2021). Na Tabela 3, Sousa et al.(2021) apresenta de 
que forma o FMB influenciou na degradabilidade ruminal. 











FMB (%) 
E 0.0 7.5 15.0 22.5 
A 265.91 A 258.1 B 257.99 B 250.66 C 
CT (horas: minutos) 01:50 B 02:40 A 02:41 A 02:29 A 
u 0.05993 0.05678 0.05631 0.05324 
DMED (5%) 70.99 66.57 64.29 65.98 
DMED (6%) 69.00 64.21 61.83 63.85 
DMED (7%) 67.05 61.92 59.46 61.76 
DMED (8%) 65.16 59.7 5717 59.74 





Tabela 3. Parâmetros da cinética de fermentação ruminal de dietas com farelo de mesocarpo de 
babaçu (FMB) determinados pelo modelo France e degradabilidade efetiva. (Sousa et al., 2021) 


A: Gás total (mL); CT: tempo de colonização (horas); u: taxa de degradação fracionária (h-1); DMED: 
Degradabilidade efetiva da matéria seca; Parâmetros estimados por France et al. (1993). Os valores 
acompanhados pelas mesmas letras nas linhas são idênticos pelo teste de identidade da curva de 5% 
(Regazzi & Silva, 2010). Fonte: Sousa et al., 2021. 


Os resultados apresentam que há uma maior degrabilidade dos alimentos quando 
está com o potencial máximo de produção de gás, que ocorre quando há uma melhor 
fermentação (SOUSA et al., 2020). Além, de comprovar que o FMB apresenta baixo valor 
nutritivo devido a alta concentração de lignina, o que o torna um alimento não tão viável 
em relação a aspectos nutritivos para alimentação de ruminantes, já que estes reduzem 
a degrabilidade. Porem, apresenta-se economicamente viável devido o fácil acesso ao 
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subproduto e por apresentar baixo custo quando comparado a outras forrageiras. Contudo, 
quando usado em pequena percentagem, o FMB se torna uma boa alternativa para 
substituição de outros substratos, desde que observados aspectos econômicos (SOUSA, 
et al., 2021). 


Dendê 


O dendezeiro é uma espécie exótica originária da Costa Ocidental da África e que 
foi introduzida no Brasil no século XVII. Atualmente é muito encontrado na Região Norte do 
Brasil, sendo utilizada e cobiçada pela indústria com finalidade para produção de biodiesel, 
justamente por ela produzir mais óleo do que qualquer outra palma. A produtividade do 
dendezeiro só na Região Amazônica, já é superior a 5 t de óleo por hectare/ano (MMA, 
2006). Segundo Rodrigues Peres et al. (2005), o dendê pode se tornar uma das principais 
fontes de biodiesel no Brasil, cujo mercado internacional é muito bem estruturado, 
permitindo um fluxo contínuo de subproduto (Sousa et al., 2010). 

Após o processo de extração do óleo para fins de biodiesel, tem-se a torta de 
dendê, também conhecida como torta de palmiste, que pode ser utilizada na alimentação 
de animais, afim de substituir outros substratos. Uma das vantagens da utilização deste 
como substrato, é por estar disponível permanentemente durante todo o ano (Pimentel et 
al., 2015). 

Pimentel et al. (2015), realizou um experimento onde foram utilizadas 12 vacas 
mestiças, com quatro tratamentos: 0% = controle (sem inclusão da torta de dendê na dieta); 
5% de inclusão de torta de dendê na matéria seca da dieta; 10% de inclusão de bolo de 
palmiste; e 15% de inclusão de bolo de palmiste. O principal objetivo do experimento era 
avaliar o efeito da inclusão da torta de palmiste na alimentação de animais. A composição 
química do bolo de palmiste reportado por Pimentel et al., 2015 é apresentada na Tabela 4. 





Componentes Bolo de palmiste 
Matéria seca (%) 92,25 
Proteína bruta (Y%DM) 14,34 
Extrato etéreo (%DM) 10,56 
Detergente neutro fibra cp'(Y%DM) 65,63 
Fibra de detergente ácido (Y%DM) 46,42 
Lignina (Y%DM) 18,31 
Matéria mineral (YDM) 3,13 


Tabela 4. Composição química do bolo de palmiste (Pimentel et al., 2015). 


Fibra em detergente neutro corrigida para cinzas e proteína. Fonte: Pimentel et al., 2015. 


A torta de palmiste se mostra uma ótima fonte de proteína e com potencial 
devido o teor consideravel de extrato etéreo. Outro ponto a ser levado em consideração 
é a composição lignocelulósica, pois como já visto anteriormente, esta pode afetar a 
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degrabilidade do alimento. 


Nível de bolo de palmiste (%DM) 


Eficiência de alimentação — >>> >> Eq! CV%? ps 
0 5 10 15 

DMI * (kg dia”) 14.71 14.89 15.74 15.52 18 7.93 0.133 
INDFcp * (kg dia”) 6.57 6.49 6.98 6.98 19 8.57 0.082 
TDNI* (kg dia?) 9.56 9.87 10.23 10.31 20 7.87 0.096 
FE “(gDM h?) 2911.59 2825.44 3051.49 3360.32 21 21.95 0.238 
FENDFc * (gNDFc h?) 1297.80 1234.07 1351.96 1509.23 22 22.02 0.153 
FETDN (gTDN H1) 1894.33 1871.55 1978.33 2231.28 28 22.77 0.214 
RE 'º (gDM H1) 1730.18 1914.05 1908.74 1882.69 24 9.40 0.047 
RENDFc " (gNDFc h 1) 771.91 833.27 846.84 846.55 25 10.11 0.108 
RETDN ?2 (gTDN h?) 1186.33 1372.62 1374.06 1480.76 26 10.95 0.028 
TCT “ (h day?) 830.00 798.75 820.42 802.92 27 7.57 0.568 
NRB '!* (número dia?) 546.85 488.30 529.47 512.10 28 11.19 0.106 
NCd '* (número dia) 30706.60 27702.47 29381.29 29178.79 29 10.94 0.178 
NCb *º (número dia?) 56.24 56.64 55.39 56.93 30 23.85 0.781 

TRB ” (sec bolus) 57.51 58.48 56.69 58.58 31 23.00 0.5 


Tabela 5. Parâmetros de eficiência alimentar e ruminação de vacas em lactação recebendo diferentes 
níveis de torta de dendê reportados por Pimentel et al., 2015. 


'Equações de regressão; 2coeficiente de variação em porcentagem; *probabilidade de erro; “DMI- 
consumo de matéria seca; * INDFcp- ingestão de fibra em detergente neutro corrigida para cinzas e 
proteínas; “TDNI- ingestão de nutrientes digestíveis totais; 7FE- eficiência de alimentação de matéria 

seca; FENDFc- eficiência de alimentação de fibra em detergente neutro corrigida; *FETDN- eficiência 
alimentar de nutrientes digestíveis totais; ''RE- eficiência de ruminação da matéria seca; "RENDFc- 
eficiência de ruminação da fibra em detergente neutro corrigida; '“RETDN- eficiência de ruminação dos 
nutrientes digestíveis totais; !ºTCT- tempo total de mastigação; '*NRB - número de bolos ruminados por 
dia; ''NCd- número de mastigações por dia; '“NCb- número de mastigações por bolus e !''TRB- tempo 
gasto por gasto do bolus ruminado. 18Y = 15,21; ºY = 6,75; ?Y = 9,99; 2Y = 3037,21; 2Y = 1348,26; “Y 
= 1993,87; *Y = 45,222 + 1745,86x, R2 = 0,45; SY = 824,64; *Y = 1132,14+ 88,47x, R2 = 0,49; 2Y = 
813,02; *Y = 519,18; ?Y = 29242,29; Y = 56,30; 'Y = 57,80. Fonte: Pimentel et al., 2015. 


Avaliando os resultados obtidos por Pimentel et al. (2015), pode-se observar que 
o tempo total de mastigação (TCT) diminui de acordo com o aumento da quantidade 
de substrato utilizado na ração. Outro fator relevante, é que, tanto o número de bolus 
ruminados por dia (NBR ), quanto o número de mastigações por dia (NCd) foram menores 
quando colocado 5% da torta de palmiste. Avaliando a eficiência de ruminação da matéria 
seca (RE) foi observado que ele diminiu consideravelmente no tratamento com 15% da 
torta. Conclui-se então que a inclusão da torta de dendê não acarreta em um aumento 
significativo no aumento do tempo despendido com alimentação, não alterando a eficácia 
da alimentação e ruminação quando utilizado até 15% da torta na ração (PIMENTEL et al., 
2015). 
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Subproduto da munguba (Pachira aquatica aubl.) com potencial para a 
aplicação na nutrição animal 

A munguba é uma espécie arbórea da familia das Malvaceae (Bombaceae), nativa 
da Floresta Amazônica e com ampla distribuição nas Américas (LORENZI, 2000; SANTOS 
et al., 2007; LIMA et al., 2012). Apresenta grande potencial para a restauração de áreas 
úmidas, solos degradados e florestas (Hernandez-Montero e Sosa, 2016; Camacho et al., 
2018), isso se dá devido suas características e capacidade de se adaptar em condições 
edáficas (Peixoto e Escudeiro, 2002; Silva et al., 2015). O fruto (Fig.2) dessa espécie 
contêm amêndoas organolépticas, ou seja, com alto teor lipídico (LORENZI, 2008; SILVA 
et al., 2010; JORGE; LUIZA, 2012; SILVA et al., 2012; CORREIA, 2019). Essa espécie tem 
ganhado mais espaço tanto na indústria de biocombustíveis e nutrição não convencional. 





Figura 2. Espécie jovem, fruto e flor da Pachira aquatica. 


Fonte: Anjos, 2020. 


Em estudo das propriedades funcionais das proteínas das amêndoas da Pachira 
aquatica Aubl. constatou-se que as suas sementes contêm aproximadamente 13,75% 
de teor proteico quando in natura, enquanto na forma de torta observou-se 28,27% de 
proteínas (SILVA et al., 2015; CORREIA, 2019). Este teor de proteína é expressivo quando 
avaliado a quantidade necessária para alimentação de ovinos, inclusive as ovelhas que 
estão em desmame que necessitam de uma proporção maior de proteína (SILVA et al., 
2008). Silva et al. (2008) também afirma que “ Da parição até o desmame, as ovelhas 
devem ser alimentadas com volumosos de boa qualidade e ração com 14-16% de proteína 
(400-800 g/dia), dependendo do tamanho da ovelha, número de crias e estado corporal.” 

A amêndoa da munguba apresenta uma composição riquissima antes e após 
extração do óleo, com alto teor de lipídeos, proteínas e carboidratos (Azevedo et al., 2015). 
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Essa é outra característica que torna esse subproduto potencial para a nutrição animal. Na 
Tabela 6 é apresentada a composição centesimal das amêndoas da munguba reportada 
por Azevedo et al., 2015. 


Componentes Teor (%) 
Umidade 5,58 + 0,09 
Cinza 4,89 + 0,09 
Lipídeos 46,62 + 0,85 
Proteínas (N x 6,25) 13,75 +0,46 
Carboidratos * 29,20 


Tabela 6. Composição das amêndoas de munguba (Pachira aquatica aubl.) (Azevedo et al., 2015) 
Carboidratos por diferença (100 menos o somatório dos demais componentes). 
Resultados das análises com média e desvio-padrão de três repetições. 


Fonte: Azevedo et al., 2015. 


A munguba apresenta altos teores nutricionais, e quando suas amêndoas são 
comparadas com outras oleaginosas, indicaram um percentual de proteínas inferiore a 
algumas espécies como a colza, soja (Barcelos et al., 2002; Azevedo et al., 2010). Os 
resultados mostram-se semelhantes a fava de morcego (13,3%) (Queiroga Neto V., 
2005; Azevedo et al., 2010), e superiores ás amêndoas da macadâmia (9,2%) (Ribeiro 
D., 2003; Azevedo et al., 2010). Já quando levado em consideração as espécies citadas 
nesse artigo, o teor de proteína é semelhante ao do bolo de palmiste, e superior ao de 
babaçu. A comparação com teores nutricionais de outras espécies é importante para a 
seleção da Pachira aquatica aubl. como alternativa na alimentação de ruminantes, devido 
as semelhanças entre elas e qual as procedências a serem tomadas para o uso de uma 
biomassa sem respaldo científico. 

Mesmo após o processo de extração do óleo, essa biomassa apresenta um alto teor 
protéico. Avaliando algumas características no geral, temos que o conteúdo protéico destas 
tortas é relativamente alto (35%), com variação de 14 a 60%, sugerindo a utilização como 
fonte de proteína para os animais (Jardim, 1976). As frações protéicas que constituem 
a proteína é um dos fatores que influenciam diretamente o resultado final, sendo ele o 
segundo nutriente limitante (SILVA, 2008). Na tabela 7 é possível conferir os percentuais 
de proteínas obtidas a partir do isolamento proteico em um estudo realizado por Silva et 
al. (2015). 
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ISOLADO PROTEICO 











IP2,0 IP 10,0 
Massa (9) (%) Massa (9) (%) 
Proteína total na farinha ' 

desengordurada (N x 6, 25) 28,27+0,25 100 28,27+0,25a 100 
Proteína extraída 10,89+0,60º 38,52 23,20+0,55b 82,06 
Precipitação Isoelétrica 6,60+0,20º 23,35 20,05+0,93b 70,94 
Proteína não Precipitada 4,30+0,69º 15,17 3,14+0,20b 11,12 
Proteína não extraida 17,38+0,66º 61,58 5,07+0,58b 17,94 


Tabela 7. Extração e recuperação de proteínas na obtenção dos isolados proteicos de 100 g da farinha 
desengordurada de amêndoas da munguba (Pachira aquatica Aubl.) (Silva et al. (2015)). 


Resultados das análises com média e desvio-padrão de cinco repetições. Letras diferentes na mesma 
linha diferem significativamente, de acordo com o teste t-student, ao nível de 5% de probabilidade. 


Legenda: IP 2,0- isolado proteico em pH 2,0; IP 10,0- isolado proteico em pH 10,0. 
Fonte: Silva et al., 2015. 


Silva et al. (2015) concluiu que o melhor rendimento da extração dos isolados 
proteicos foi para o IP 10,0. Além, das sementes de munguba apresentar resultados 
promissores em relação ao teor lipídico e proteico. A Tabela 8 apresenta detalhadamente 
a distribuição das frações proteicas da farinha desengordurada reportado por Azevedo et 





al., 2010. 
Frações Proteicas Índices (%) 
Albumina 22,86 + 0,02 
Globulina 56,24 + 0,26 
Prolamina 1,43 +0,01 
Gluetinas 14,22 +0,14 
Resíduo 5,36 


Tabela 8. Índices de proteínas das frações proteicas das amêndoas da munguba a partir da farinha 
desengordurada, de acordo com sua solubilidade (Azevedo et al., 2010). 


Resultados das análises com média e desvio padrão de três repetições. 


Fonte: Azevedo et al., 2010. 
O conteúdo de proteínas na farinha desengordurada, extraídas e recuperadas 


caracterizam a amêndoa da munguba como boa fonte de proteína (Azevedo, 2008). A 
globulina é uma das várias globulares que são pouco hidrossolúveis e é a fração principal 
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em maior destaque com o índice de 56,24 %, seguido pela albumina com 22,86%. 

O teor de gordura (EE) também varia consideravelmente (3 a 24%) o que pode ser 
outro benefício para os ruminantes, considerando que a inclusão de óleo na dieta pode 
auxiliar na mitigação de metano entérico (Grainger, 2008). 


41 CONCLUSÃO 


É notório que os subprodutos de espécies oleaginosas apresentam características 
nutricionais essenciais na alimentação de ruminantes. Tornando-se uma alternativa para a 
substituição de outros alimentos, e tendo uma aplicabilidade para o subproduto que antes 
era descartado sem finalidades nenhuma, agregando valor para essa biomassa residual. 

As pesquisas apontam os subprodutos da munguba como promissoras tanto 
economicamente quanto nutricionalmente na dieta animal. Mas ainda há poucas 
informações na literatura sobre o uso desses subprodutos nos sistemas alimentares. 
Sugere-se o desenvolvimento de trabalhos a respeito da caracterização, se há toxinas 
prejudiciais, e a aplicação direta desses subprodutos na alimentação de ruminantes. 
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